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FOREWORD 


This  report  was  prepared  by  Bolt  Beranek  and  Newman  Inc. , 
Cambridge,  Massachusetts,  for  the  Structural  Integrity  Branch, 
Structural  Mechanics  Division,  Air  Force  Flight  Dynamics  Labora- 
tory, Wright-Patterson  Air  Force  Base,  Ohio  under  Contract 
F33615-75-C-3017 . The  work  described  herein  is  a part  of  the 
Air  Force  System  Command  continuing  program  to  establish  methods 
of  predicting  and  controlling  the  aero-acoustic  environment  of 
flight  vehicles.  The  work  was  directed  under  Project  1471, 
"Aero-Acoustic  Problems  in  Flight  Vehicles,  Task  147102  Aero- 
Acoustics."  Messrs.  Robert  Gordon  and  Davey  L,  Smith  of 
AFFDL/FBE  served  as  Air  Force  Project  Engineers,  Their  cooper- 
ation and  assistance  are  gratefully  acknowledged. 

The  work  presented  herein  was  performed  by  several  members 
of  the  staff  of  Bolt  Beranek  and  Newman  Inc. , under  the  techni- 
cal direction  of  Dr.  Eric  E,  Ungar , who  also  was  responsible 
for  integration  and  presentation  of  the  results.  Dr.  John  Wilby 
and  Mr.  B,  Pinkel  contributed  the  portions  dealing  with  jet 
noise,  and  Dr.  A,  Galaitsis  provided  the  part  addressing  fan/ 
compressor  noise.  Mr.  Richard  Hayden  prepared  the  part 
pertaining  to  powered  lift  devices,  Mr.  Joseph  Smullin  contrib- 
uted the  section  dealing  with  propeller  noise,  and  Mr.  Robert  White 
developed  the  portion  concerned  with  armament.  Dr.  Donald  Bliss 
contributed  the  sections  dealing  with  surface  flows  and  cavity 
noise. 


Background  information  on  which  this  report  is  based  is 
presented  in  a companion  volume,  "A  Review  of  Methods  for 
Estimation  of  Aeroacoustic  Loads  on  Flight  Vehicle  Surfaces,” 
AFFDL-TR-76-91,  Vol.  II. 


This  report  concludes  the  work  on  Contract  F33615-75-C-3017 , 
which  covered  a period  from  January  1975  to  October  1976.  The 
manuscript  was  released  by  the  authors  in  October  1976. 
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SECTION  I 


INTRODUCTION 

The  fluctuating  pressures  that  act  on  the  surfaces  of  flight 
vehicles  due  to  noise  from  propulsion  or  powered  left  systems,  due 
to  surface  flows,  and  due  to  armament  constitute  potentially 
Important  loads  that  need  to  be  considered  in  the  design  process. 

In  order  to  compile  a guide  for  predicting  the  significant 
aeroacoustic  loads  on  flight  vehicles  that  represents  the  best 
currently  available  practice,  an  extensive  review  of  the  litera- 
ture was  undertaken  and  alternative  prediction  techniques  were 
evaluated.  The  results  of  these  reviews  and  evaluations  are  pre- 
sented in  the  present  report,  together  with  some  extensions, 
clarifications,  simplifications,  and  corrections  of  previous 
methods . 

Although  considerable  effort  has  been  made  to  make  this  re- 
port self-contained  and  meaningful  to  the  nonspecialist  in  aero- 
acoustics,  no  attempt  has  been  made  to  include  the  fundamentals 
of  acoustics  and  aerodynamics. 

In  order  to  conserve  space,  computer  programs  have  not  been 
reproduced;  since  such  programs  are  useful  only  with  extensive 
documentation,  the  reader  is  referred  to  the  original  sources. 
Similarly,  the  data  presentations  included  here  are  limited  to 
those  that  illustrate  important  trends  or  that  serve  as  the  bases 
for  empirical  prediction  methods;  for  more  extensive  or  detailed 
data  collections,  the  reader  is  again  referred  to  the  original 
literature. 

The  present  volume  does  not  include  discussions  of  two  topics 
appearing  in  the  estimation  guide  volume,  namely:  fsm/compressor 

noise  and  noise  of  Internal  combustion  engines.  Discussion  of 


fan/compressor  noise  has  been  omitted  here,  becav*’''  a aum..j.ent 
treatment  appears  in  the  reoert  I’te-^ature  and  because  a dis- 
cussion of  the^  background  and  physical  mechanisms  was  Included 
In  the  estimation  guide  volume,  where  this  discussion  also  pro- 
vides some  guidance  for  the  user  of  the  prediction  methods 
presented  there.  Discussion  of  Internal  combustion  engine  noise 
has  been  omitted  here,  because  no  relevant  information  has  become 
available  since  appearance  of  the  publication  from  which  the 
suggested  prediction  technique  has  been  taken. 


\ 
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SECTION  II 
^IFT  NO!S: 


2..1  Introductlori 

This  section  deals  with  the  prediction  of  noise  from  the 
exhausts  of  turbojet  and  turbofan  engines,  including  shock- 
free  noise  and  shock  noise  from  nozzles  of  any  practical  shape 
(circular,  slot,  annular,  coaxial,  etc*),  but  excluding  the 
noise  associated  with  suppressor  nozzles,  such  as  multitube 
and  multilobe  designs.  Afterburner  and  thrust  reverser  noise 
is  also  included,  but  such  other  sources  of  combustion  noise 
as  core  noise  are  not  considered  here.  Discussion  of  noise 
radiation  from  jet/airframe  interactions  is  here  confined  to 
the  case  of  thrust  reversers. 

Most  of  the  published  methods  for  the  prediction  of  exhaust 
noise  have  been  concerned  with  farfield,  rather  than  nearfleld, 
sound  pressure  levels.  However,  for  the  case  of  flow  from  cir- 
cular exhaust  nozzles,  several  prediction  methods  have  been 
developed  during  recent  years.  Since  these  prediction  methods, 
of  necessity,  need  to  servo  as  the  basis  for  nearfield  predic- 
tion methods  of  other  exhaust  flow  conditions,  special  emphasis 
is  given  to  these  methods  in  the  discussion. 

Although  none  of  the  available  prediction  methods  Is  ac- 
ceptable without  reservations,  the  procedure  of  Plumblee  et  at 
has  been  selected  to  serve  as  the  basis  for  nearfleld  noise 
predictions.  Modifications  to  this  method  are  developed  to  ac- 
count for  nozzle  shape,  afterburner  effects,  coaxial  exhausts, 
shock  noise  and  forward  motion.  However,  for  thrust  reversers, 
a new  model  is  constructed,  since  this  noise  source  mechnism 
differs  from  the  others. 
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Reflection  and  scattering  effects  associated  with  the  pres~ 
ence  of  the  airframe  structure  are  also  discussed,  pressure 
correlation  data  on  the  airframe  structure  are  reviewed,  and  a 
related  prediction  model  Is  proposed. 

2.?  Subsonic  and  Shock-Free  Supersonic  Exhausts 
2.2.1  Circular  nozzles 

The  methods  for  predicting  nearfleld  noise  fall  into  three 
categories: 

1.  Analytical  methods.  In  which  nearfleld  noise  estimates 
are  derived  from  basic  flow  parameters. 

2.  Curve-fltt Ing  methods,  in  which  equations  are  fitted 
to  test  data. 

3.  Scaling  methods,  in  which  a set  of  test  data  is  taken 
as  the  reference  base  and  rules  are  developed  for  adjusting  the 
reference  values  to  account  for  changes  in  the  important  param- 
eters. 

Analytical  Methods 

Chen,  Benzakein  and  Knott  [l,  2,  3]*  have  studied 
the  applicability  of  equations  derived  by  Pranz  [^].  These 
equations  describe  the  acoustic  pressure  field  of  an  axlsymmetrlc 
Jet,  as  a function  of  the  distance  from  the  Jet  and  the  sound 
power  radiated  by  it  per  unit  volume.  Franz's  equations  are  de- 
rived for  three  different  assumptions  regarding  the  structure  of 
the  Jet  flow  field;  (1)  laterally  oriented  quadrupoles,  (2). longi- 
tudinal quadrupolos,  and  (3)  isotropic  turbulence.  The  equation 
for  isotropic  turbulence  is  the  simplest  of  the  three  and  fits 
the  test  data  nearly  as  well  as  the  other  two  cases.  This 

^Numbers  in  brackets  refer  to  list  of  references. 
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equation » reproduced  here  for  purposes  of  Illustration,  Is 


p c W 2c® 

-2_2_  1 + __o_.  + 12 

l<irr®  r®(*)® 


r-w** 


(1) 


The  quantity  W denotes  the  radiated  sound  power  per  unit  volume 
and  Is  given  by 


W 


B 

1 


p®u^‘*V‘* 


(2) 


Where  u'  is  the  fluctuating  component  of  the  longitudinal  velo- 
city, V is  the  steady  velocity  component,  and  t is  the  local 
scale  of  tho  turbulerice  in  the  Jet.  Chen  et  at  determine,  for 
each  1/3  octave  band,  the  noise  at  a point  In  the  nearfleld  by 
integrating  the  contributions  from  all  positions  in  the  Jet, 
making  use  of  data  on  the  distribution  of  u^",  V,  and  i,  with  posi 
ticn  in  the  Jet. 

Comparisons  [1]  between  values  computed  by  this  method  and 
experimental  data  reveal  differences  In  some  areas  of  as  much  as 
10  dB.  Because  of  the  possibility  of  large  errors,  and  because 
of  its  complexity,  application  of  this  method  for  prediction  pur 
poses  is  not  recommended. 
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Maestrello’8  analytical  method  [5]  la  based  on  first  prin- 
ciples, but  needs  considerable  additional  development  to  be 
applicable  to  practical  nearfleld  acoustic  pressure  calculations • 

Curvs-Fitting  Methods 

Plumblee  et  at.  [6]  start  with  the  Llghthlll  equation  and 
Introduce  Into  It:  (a)  directionality  terms  suggested  by 

Williams  [7],  Rlbner  [8],  and  Lllley  C91  and  additional  direction- 
ality terms  needed  to  improve  agreement  with  test  data,  and  (b) 
an  expression  for  the  effect  of  distance  r,  the  form  of  which  Is 
suggested  by  Equation  (1)  obtained  by  Pranz  [4].  The  relationship 
synthesized  In  this  manner  Is 


?*(T,M,r,G) 


m 

1 

-C  9 

C e * \ 

#1  _ M cos  9 ’ 

\ +a*M* 

ll  + - 

- 

1 

V * 

-C  r/k 

l«e  ’ 

/ J 

1 \ 1+C  e " / 

* 

(3) 


with  the  parameters  a*,  and  given  by  empirically  deter- 

mined functions  of  M and  T,  by  means  of  which  Plumblee  et  al. 
attempt  to  fit  this  equation  to  an  extensive  set  of  nearfield 
acoustic  pressure  data,  accounting  for  the  effects  of  Jet  velo- 
city, temperature  and  Mach  number  on  a 3.5  inch  diameter  Jet. 
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Scaling  Methode 

Two  approaches  have  been  presented:  one  by  Thomson  [10], 

•and  one  by  cockburn  and  Jolly  [11]  using  the  test  data  and 
methods  presented  by  Hermes  and  Smith  [12],  Sutherland  and  Browr 
[13]  discuss  methods  of  correcting  the  nearfield  data  for  gi*ound 
reflections.  Rudder  and  Plumblee  [l4]  present  resumes  of  the 
Thomson  and  the  Cockburn  and  Jolly  methods,  as  well  as  of  the 
Plumblee  empirical  equation. 

In  the  empirical  methods  [10,  11],  a set  of  nearfleld 
noise  contours  is  established  as  a reference  set,  and  rules  are 
provided  for  computing  adjustments  to  these  values  for  changes  Ir 
jet  velocity  and  for  determining  noise  frequency  spectra.  The 
reference  sets  of  noise  contours  for  the  two  scaling  methods 
correspond  to  the  following  conditions: 

Vj  (ft/sec)  T^  static  (°R) 

Thomson  [10]  2000  1252 

Cockburn  [11]  1920  1130 

The  reference  contours  for  the  Cockburn  method  consist  of  the 
test  data  obtained  by  Hermes  and  Smith  on  a J57-P21  engine  at 
100/S  military  thrust. 

Accuracy  of  Empirical  Methods 

It  is  of  considerable  Interest  to  explore  the  accuracy  with 
which  the  several  empirical  methods  mentioned  above  predict  ex- 
perimental results.  Published  test  data  references- and  the  range 
of  exhaust  jet  conditions  covered  by  them  are  given  in  Table  1. 


7 


! 


\ 

t 

f 


1 

m 

o\ 

0) 

0) 

bOW 

m 

m 

(0 

c 

rH 

B bO 

• 

CM 

X) 

o 

iH 

P 

•H  C 

fV 

« 

rn 

Ctt 

bO<H 

O 

o 

p 

B 

P 

U bO 

o 

P 

o 

(d 

CQ 

iH 

0) 

CQ 

O B 

p 

43 

• 

0) 

1 1 

3 

> 0) 

o 

o 

IT\ 

0)  0) 

(0 

> 

p 

1-j  0) 

E > 

OJ 

ir» 

Oi 

f\J 

x:  x: 

p 

(d 

o 

o 

U 

O 

• 

oo 

C\J 

• 

p p 

0) 

0) 

K 

< 

m 

O XI 

CM 

iH 

rH 

C\J 

t. 

<■4 

•p 

o 


o o 
^-1  <-( 


t> 


Xi 

p 


•r^ 

VO 

O 

VO 

CM 

>»  w 

0 

0 

rl 

O 

1 

CM 

o\ 

a\ 

00 

B 4) 

CO 

r| 

IS 

X 

4> 

bOai 

• 

CM 

rH 

cd  3 

B 

P 

u 

rH 

C bO 

CM 

CM 

rH 

rH 

CO  rH 

4) 

rl 

■o 

0 

3 

X 

iH  c 

10  (d 

a 

T) 

4) 

P 

cd 

b0*H 

O 

o 

o 

o 

4)  > 

3 

C 

(0 

P r-i 

p 

B to 

p 

P 

4^ 

43 

0 

CO 

0 

3 

<tt  VO 

OV 

m 

4)  B 

4>  » 

CJ 

iH 

C- 

3 (d 

> 4) 

rH 

lA 

o 

rH 

m 

c 0 

B 

4) 

X ' 

1 

T3  0) 

E > 

OO 

OV 

H 

0 

B 

B 

O 

^-s 

t3  B 

O -H 

• 

lA 

CM 

. 

P 

Vi 

6) 

4) 

<c 

o 

bH 

< <d 

O TJ 

rH 

rH 

rH 

1 — 1 

0 T3 

B 

h- 

B C 

bO 

B 

< 

cd 

B 

3 

CO 

O 

4> 

rl 

X 

4; 

Th  X 

to 

B 

rH 

f- 

CO 

o 

c^.  bO 

B 

4) 

N 

t/) 

1 

CM 

VO 

oo 

in 

rl 

(U 

P 

X 

UJ 

a 

CO 

4) 

to« 

• 

C\J 

o\ 

B P 

> 

V, 

0 

h- 

1 

rH 

C bO 

CO 

CM 

• 

OJ 

4> 

rl 

Cd 

E 

p 

U 

P 

•H  C 

CO 

rH 

P 4) 

TJ 

Um 

Cj  t — 1 

o 

(d 

bO<H 

o 

o 

4>  > 

1 

B 

to 

O 

ir\ 

1 

P 

B bO 

p 

o 

O 

4<> 

B mH 

to 

0 

B 

X <~i 

CO 

w 

4)  B 

P 

4> 

cd  P 

B 

Vi 

rl 

</J 

p 1— 1 

ov 

3 <d 

> (U 

m 

m 

B 0 

H 

to 

Ui 

•H 

1 

r-i  q} 

C > 

oo 

•sj" 

o^ 

GO 

cd  4> 

bO 

B 

B 

o 

E 

•‘J 

•O  B 

O H 

• 

■ir 

o 

a a 

B 

4> 

4> 

ck; 

CO 

X 

< <d 

O T3 

CM 

•=r 

t 

rH 

10 

4) 

> 

f» 

4)  0> 

> 

H 

rl 

o 

B B 

B 

bO 

X) 

40 

0 

0 

1 

4> 

0 

P 

to 

VO 

>,X 

0 

B 

• 

5>> 

CM 

o\ 

m 

C P 

•O 

B 

■H 

iH 

CO 

in 

c^ 

cd 

B 

to 

4) 

a 

■sr 

• 

m 

(0 

Cd 

0 

B 

UJ 

4) 

a 

rH 

B cd 

rl 

4> 

rH 

3 

m 

o 

o 

0 

CQ 

P 

> 

CO 

■9 

CO 

r-N 

CM 

p 

o 

43 

Vi  10 

P 

Cd 

B 

< 

E VO 

<d 

CJV 

p 

0 

B 

4) 

B 

0 

h- 

3 1—1 

> — ' 

CM 

VO 

in 

CO  0 

r> 

0 

rH 

•H 

• 

CO 

oo 

tw 

4>  *H 

4> 

o< 

< 

o 

CTl 

in 

3 P 

0 

B 

•1 

rH  rl 

•H 

3 

CO 

Cd  TJ 

CO 

s 

r-> 

> B 

0 

CO 

0 

P 

0 . 

CO 

4) 

rH 

Vi 

> 0 CO 

X 

B 

4h 

V ' 

0 B 

3 

a 

H U)  4) 

CO 

0 

to 

E P 

n 

•H 

4) 

r— t 

0) 

•0  rl  4> 

B 

E 

B 

P 

o 

U 

B P 1 

0 

cd 

0 

4) 

4) 

3 

Cd  Cd  CO 

Vi 

CO 

e 

(Q 

P 

4> 

B 

(d 

\ 

<d 

CQ 

- 

X3  4)  cd 

bO 

B 

• 

4) 

tH 

P 

u 

0 

to  a a 

B 

3 

B 

a 

TJ 

Vi 

•H 

rl  0 

rl 

P 

4> 

3 

> — i 

Q.K 

P 

X B 

to 

«d 

> 

CO 

4) 

p 

8 » 

(d 

0)  41 

B 

B 

mH 

O 

rH 

IQ 

4> 

4)  B X5 

4> 

4) 

(0 

rH  ' 

0) 

K 

4) 

X 

P 

x:  01  p 

> 

a 

f— 1 ' 

O 

3 

4) 

rH 

E 

O 

4> 

Eh  W 0 

B 

B 

p 

(d  ' 

B 

O 

P 

3 

P rl 

0 

4) 

0 

4) 

to 

4) 

4) 

rl 

X 

« P 

3 

0 

Eh 

B 

rH 

rH 

V 

3 cd 

n 

• • 

0 

0) 

» 

N 

H 

O 

X 

Id  p 

n 

4) 

Pi 

4h 

O 

N 

H 

rH 

o 

x:  « 

4) 

P 

r"\ 

r-x 

«) 

rH 

O 

O 

4> 

M 

0 

Cd 

X 

0 

b 

SB 

S 

M 

Pm 

as 

m 

8 


providing  shock- free  flow. 


The  Plumblee  data  cover  a much  wider  range  of  conditions  than 
are  of  Interest  for  practical  turbojet  and  turbofan  engines  at 
takeoff  and  cruise  thrust.  To  provide  some  Insight  into  the 
range  of  conditions  of  practical  Interest,  Table  2 lists 
take-off  conditions  for  a turbojet  and  various  turbofans  based 
on  a level  of  technology  corresponding  to  the  J-93  engine. 

The  Jet  conditions  for  other  levels  of  engine  technology  will 
differ  from  those  listed  in  this  Table,  but  may  lie  expected 
largely  to  fall  within  the  range  of  conditions  indicated. 

Table  3 compares  overall  sound  pressure  levels  computed  by 
the  methods  of  [10],  [11],  and  [6]  with  each  other  and  with 
measured  data,  for  several  different  measurement  locations,  for 
a Jet  with  Vj  = 2221  ft/sec,  T^  static  “ 1360°R.  The  calcula- 
tions here  involve  only  small  corrections  for  Vj  and  Tj , because 
the  selected  engine  conditions  correspond  nearly  to  those  of 
the  reference  contours.  The  method  proposed  by  Cockburn  involves 
a correction  for  Jet  velocity  only;  the  Thomson  method  involves 
corrections  for  both  Jet  velocity  and  Jet  temperature. 

Plumblee 's  measured  and  computed  values  are  in  good  agreement 
with  each  other.  The  Cockburn  method  predicts  the  highest 
values.  The  Thomson  method  predicts  noise  values  for  the  most 
part  about  ^ dB  below  the  Cockburn  method,  and  the  Plumblee 
method  predicts  values  about  8 dB  below  the  Cockburn  method. 

The  Cockbui'n  reference  noise  contours  are  based  on  measurements 
made  on  a Jet  engine.  These  values  contain  some  machinery  noise 
and  would  be  expected  to  be  higher  than  values  determined  from 
the  measurements  by  Plumblee  on  a model  Jet  attached  to  an  air 
supply. 
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TABLE  3.  COMPARISON  OF  COMPUTED  AND  MEASURED  SOUND  PRESSURF  >^OR* 
Vj  - 2221  ft/sec  and  Tj  " 13®0«R 


Measurement 

Coordinates 

Overall 

Sound  Pressure 

Level  (dB,  re  2 

X 10“'  N/m*) 

X 

y 

K 

Computed 

Measured 

U 

u 

Thomson 

[10] 

Cockburn 

[11] 

Plumb lee 
[6] 

Plumblee 

[6] 

0 

111 

130.8 

133.8 

125.5 

127 

5 

19 

130.9 

134.7 

126 

126 

10 

23.5 

131.0 

134.7 

- 

- 

15 

28 

131.1 

134.7 

- 

0 

3 

- 

139.5 

- 

- 

5 

7.5 

IHO.9 

143.7 

137.5 

137.5 

10 

11.7 

141.3 

144.1 

136.0 

135.5 

15 

15 

141.4 

145.5 

135.2 

134.0 

20 

17.7 

14i,7 

144.7 

134 

133 

23 

19 

141.7 

144.4 

133 

5 

3 

151.1 

154.4 

15c 

148 

10 

5.2 

151.9 

155.0 

144 

144 

15 

5 

151.3 

155.0 

I4l 

140 

•The  engine  operating  condition  corresponds  to  that  of  Figure  30 
of  [6j. 
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Table  3 corresponds  to  non-afterburning  operation  of  a 
turbojet  engine.  Table  ^ similarly  corresponds  to  the  upper 
end  of  the  Jet  velocity  range  of  Interest,  that  is,  approximately 
to  the  Max  A/B  conditions  of  Table  2,  for  which  large  corrections 
to  the  reference  values  In  the  Thomson  and  Cockburn  methods  are 
required. 

Again,  Plumblee's  computed  and  measured  values  are  In  good 
agreement  with  each  other.  On  the  other  hand,  the  Thomson 
method  underestimates  the  noise,  whereas  the  Cockburn  method 
overestimates  it  by  5 to  10  dB.  These  discrepancies  are  due  to 
inadequate  consideration  of  temperature  effects.  These  are  not 
Included  at  all  in  the  Cockburn  method,  and  are  overestimated  in 
the  Thomson  method.* 

It  is  not  surprising  that  the  values  computed  by  the  Plumblee 
equation  agree  well  with  most  of  the  noise  values  measured  by 
Plumblee,  because  the  equation  represents  an  empirical  fit  to 
these  data.  However,  Tables  3 and  M represent  only  a 
small  part  of  the  range  of  data  of  interest,  and  one  needs  to 
consider  whether  there  exist  regimes  where  the  Plumblee  equation 
represents  the  data  inadequately.  Comparison  of  values  computed 
by  the  Plumblee  equation  to  test  data  reveals  discrepancies  of 
less  than  5 dB  over  most  of  the  very  extensive  test  regime  cov- 
ered. The  largest  discrepancies  occur  for  supersonic 


•In  the  Thomson  method  the  correction  for  temperature  here  is 
-7.8  dB,  based  on  the  assumption  that  the  noise  varies  as  p?, 
and  that  pj  varies  inversely  as  Tj . The  Thomson  method  may*^ 
overestimate  the  effect  of  temperature.  As  pointed  out  by  Tanna 
[18],  the  temperature  that  should  be  used  is  that  near  the 
end  of  the  potential  core,  at  a position  where  the  turbulence 
is  a maximum.  The  use  of  this  temperature  to  the  second  power 
results  in  a variation  of  the  noise  intensity  with  Tj*,  where 
Tj  is  the  temperature  at  the  nozzle  exit.  This  suggests  that 
tne  correction  for  temperature  should  be  half  of  that  given  by 
the  Thomson  procedure.  If  the  values  given  in  Table  ^ under 
the  Thomson  method  correspondingly  are  Increased  by  3-9  dB,  they 
would  be  brought  above  the  Plumblee  values  by  an  amount  compar- 
able to  that  in  Table  3. 
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TABLE  4.  COMPARISON  OF  COMPUTED  AND  MEASURED  SOUND  PRESSURES  FOR* 
Vj  « 3267  ft/sec  and  Tj  “ 3069°R 

Measurement  I _5 

Coordinates  1 Overall  Sound  Pressure  Level  (dB,  re  2 x 10“  N/m^) 


14 
19 

23.5 
28 

3 

7.5 

11.7 

15 

17.5 
19 

3 

5.2 

5 


Thomson 

[10] 

128.8 

129.1 

129.2 

129.4 

139.3 
140.7 

141.6 

143.3 
143.3 
150 
153.9 

151.6 


Computed 


Cockburn 

[11] 

143.9 

144.4 

144.4 

144.6 

146.2 

153.6 
155.0 

157.4 

154.6 

156.2 
166.2 

168.5 
168.5 


Plumblee 

[6] 


Measured 


Plumblee 

[6] 


143 

145 

143 

145 

142.5 

142 

141.5 

l40 

I4l 

- 

155 

- 

151 

149 

149 

144.5 

•The  engine  operating  condition  corresponds  to  that  of  Figure  31 
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Jets*  for  values  of  X/D  leas  than  7,  and  for  the  highest  octave 
band  investigated,  as  Illustrated  In  Figures  1,2,  and  3. 


^Although  in  each  test  case  the  nozzle  had  a converging-diverging 
contour,  there  may  have  been  some  shocks  in  the  jet,  which  may 
be  responsible  for  these  discrepancies.  Figure  8 of  [6] 

Is  a photograph  of  the  f^ach  1.5  jet,  revealing  a characteristic 
diamond-shaped  shock-pattern. 


The  frequency  of  the  noise  generated  by  shock  waves  radl- 
at*  'g  at  right  angles  to  the  jet  is  given  [19] 


3D»^S  -1789 

Wj  're  ci  denotes  the  velocity  of  sound  in  air  (ft/sec),  PR  is 
the  pressure  ratio  across  the  nozzle,  and  D represents  the 
nozzle  diameter  (ft). 

For  a Mach  number  of  1.25  (i.e.,  PR  « 2.6),  the  value  of  fD 
is  -^oOO  Hz-inch.  This  fa-ls  within  the  octave  band  8^00-16800 
Hz-inch,  where  Figures  1 and  2 reveal  the  upswing  in  measured 
noise  at  low  values  of  X/D.  For  a jet  Mach  number  of  l.:5  (i.e., 
PR  = 3*7)»  the  value  of  fD  given  by  the  above  equation  is  67OO 
Hz-ir  , h,  which  falls  in  the  ^1200-8^100  Hz-inch  octave  band,  where 
Figure  3 reveals  the  upswing  in  measured  noise.  This  upswing  is 
even  more  marked  in  the  next  higher  band. 


li| 


FIG.  1.  COMPARISON  OF  MODEL  JET  EXPERIMENTAL  AND  CALCULATED  CONTOURS 
MACH  NO.  * 1.25,  TEMP  * 1360”R,  VELOCITY  = 2221  FPS 


FIG.  2.  COMPARISON  OF  MODEL  JET  EXPERIMENTAL  AND  CALCULATED  CONTOURS 
MACH  NO.  » 1.25,  TEMP  » 3069"R,  VELOCITY  - 3267  FPS 


FI6.  3..  COMPARISON  OF  MODEL  JET  EXPERIMENTAL  AND  CALCULATED  CONTOURS 
MACH  NO.  « 1.50,  TEMP  « 1360*R,  VELOCITY  • 2677  FPS 
(FIG.  34,  REF.  6). 


Figures  ^ and  5 show  a comparison  between  the  near fie Id 
noise  contours  measured  for  a YJ-93  engine  and  values  computed 
by  the  Plumblee  equation.  The  conditions  to  which  these  figures 
pertain  are  as  follows : 


Weight 

Total 

Jet 

Exit 

Exit 

Pig.  No. 

Condition 

Thrust 

flow 

Temp 

Velocity 

Area 

Mach 

(lb) 

(Ib/sec)  (°R) 

(ft/sec) 

(ft*) 

No. 

4 

Military 

19405 

283 

2066 

2466 

6.4 

1.28 

5 

Max  A/B 

26367 

288 

3651 

3?46 

8.2 

1.26 

\ 

N, 


The  agreement  between  the  computed  and  measured  v^ues  is  very 
good  for  most  of  the  range  shown.  \ 

Ihe  YJ-93  engine  j s equipped  with  a converging-diverging  nozzle 
If  a shock  wave  contribution  to  the  noise  occurs,  one  wouXd  ex- 
pect it  to  be  observed  in  the  highest  octave  band  (fD  » 8*10b  to 
168C0  Hz-lncn) . A small  upswing  of  the  measured  noise  contours 
at  low  values  of  X/D  is  indeed  noted.  Again,  the  Plumblee  equa- ' 
tlon  fails  to  reflect  this  upswing. 

Poeoiblg  Modifioations  of  Plumblee  Method 

As  is  evident  from  Equation  (3)  and  the  equations  for  the 
eleve-i  associated  constants , the  Plumblee  method  is  somewhat 
cumbersome  to  apply.  A useful  alternative  may  be  based  on  the 
use  of  two  sets  of  reference  contours  (one  for  Mach  1.25  and 
the  other  for  subsonic  Jets)  and  on  an  equation  for  adjusting 
the  noise  levels  for  changes  in  Jet  velocity  V and  jet  static 
temperature  T.  For  any  given  case,  one  would  select  the  set 
of  contours  whose  Mach  number,  velocity,  emd  temperature  are 
closest  to  the  parameters  of  the  jet  to  be  considered  and  the 
level  would  be  adjusted  by 

A dB  • 70  log  ^ - 10  log  ^ , (5) 

r r 

where  the  subscript  r designates  the  reference  parameters. 
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FIG.  4a.  COMPARISON  OF  PREDICTED  (INCLUDING  REFLECTION)  AND 
MEASURED  SOUND  PRESSURE  LEVEL  (dB)  CONTOURS  FOR 
.SINGLE  YJ-93  ENGINE  OPERATING  AT  MILITARY  POWER. 
(a)OVEBALL  BAND,  f D - 1050  to  26880  Hz  In. 

(FROM  FIG.  15,  REF.  15). 


Predicted  by  Plumblee  Equation 
Measured 


JET  AXIS 


FIG.  4b.  COMPARISON  OF  PREDICTED  (INCLUDING  REFLECTION)  AND 
MEASURED  SOUND  PRESSURE  LEVEL  (dD)  CONTOURS  FOR 
SINGLE  YJ«93  ENGINE  OPERATING  AT  MILITARY  POWER, 
(b)  OCTAVE  BAND,  fD  » 2100  to  4200  Hz  In. 

(FROM  FIG.  15.  REF.  15). 


2! 


Pndlcted  by  Plumbtee  Equotton 
Maosured 


FIG.  4c.  COMPARISON  OF  PREDICTED  (INCLUDING  REFLECTION)  AND 
MEASURED  SOUND  PRESSURE  LEVEL  (dB)  CONTOURS  FOR 
SINGLE  YJ-93  ENGINE  OPERATING  AT  MILITARY  POWER, 
(c)  OCTAVE  BAND,  fO  • 4200  to  8400  Hz  In, 

(FROM  FIG.  15.  REF.  15). 


Predicted  by  Plumblee  Equation 
Measured 


FIG.  4d.  COMPARISON  OF  PREDICTED  (INCLUDING  REFLECTION) 

AND  MEASURED  SOUND  PRESSURE  LEVEL  (dB)  CONTOURS 
FOR  SINGLE  YJ-93  ENGINE  OPERATING  AT  MILITARY  POWER, 
(d)  OCTAVE  band.  fD  * 8400  to  16800  Hz  In. 

(from  fig.  15.  REF.  15). 


Prodlctod  by  Plumblee  Equotlon 
Measured 


FIG.  5a.  COMPARISON  OF  PREDICTED  (INCLUDING  REFLECTION)  AND 
MEASURED  SOUND  PRESSURE  LEVEL  CONTOURS  FOR  SINGLE 
YJ-93  ENGINE  OPERATING  AT  MAXIMUM  AFTERBURNER. 

(a)  OVERALL  BAND.  fD  - 1050  to  26880  Hz  in. 

(FROM  FIG.  14,  REF.  15). 


m 


- V*  # ^ « 

r >.  *> 


~—  Predicted  by  Plumblee  Equation 
Measured 


,/  A 


M § 


7 «> 

/ ^ 


135- 


X/D 

JET  AXIS 


COMPARISON  OF  PREDICTED  (INCLUDING  REFLECTION)  AND 
MEASURED  SOUND  PRESSURE  LEVEL  CONTOURS  FOR  SINGLE 
YJ-93  ENGINE  OPERATING  AT  MAXIMUM  AFTERBURNER. 


(b)  fD  » 2100  to  4200  Hz  In. 
(FROM  FIG.  14,  REF,  15). 


predicted  by  Plumblee  Equation 

■ — ■ Measured 


FIG.  Sd.  COMPARISON  OF  PREDICTED  (INCLUDING  REFLECTION)  AND 
MEASURED  SOUND  PRESSURE  LEVEL  CONTOURS  FOR  SINGLE 
YJ-93  ENGINE  OPERATING  AT  MAXIMUM  AFTERBURNER. 

(d)  fD  •>  8400  to  16800  Hz  In. 

(FROM  FIG.  14.  REF.  15). 


It  Is  beyond  the  scope  of  the  present  report  to  establish 
the  reference  sets  of  parameters.  (This  would  require  careful 
consideration  of  all  of  the  available  test  data  and  their  respec- 
tive accuracies.)  To  provide  some  Insight  Into  the  possibilities 
of  the  suggested  method,  one  may  choose  the  curves  of  Figures  ^ 
(for  the  YJ-93  engine  at  military  thrust)  as  the  reference  con- 
tours for  the  regime  around  M » 1.25.  .Figure  6 shows  the  test 
data  from  Figure  together  with  those  from  Figures  5 (Max  A/B), 
adjusted  according  to  Equation  (5).  (The  calculated  adjustment 
Is  -7  dB. ) The  agreement  here  between  the  curves  for  military 
thrust  and  the  reduced  max  A/B  curves  Is  better  than  the  agree- 
ment between  these  curves  and  the  values  computed  by  the 
Plumblee  method  (Figures  ^ and  5). 

Figure  1,7  shows  a similar  comparison  of  the  test  data  on 
the  model  Jet  from  Figures  1 and  2 with  the  levels  of  the  curves 
of  Figure  2 reduced  by  8 dB  (as  computed  from  Equation  5),  to 
adjust  these  curves  to  the  conditions  corresponding  to  Figure  1. 
Again,  the  agreement  Is  seen  to  be  very  good,  and  somewhat  better 
than  that  between  the  test  data  and  the  values  computed  by  the 
Plumblee  equation  (Figures  1 and  2). 

Figure  8 compares  two  sets  of  test  data  from  [6]  with  the 
levels  of  the  curves  of  one  set  increase  7 dB  in  accordance  with 
Equation  (5)  to  account  for  differences  in  the  parameters.  The 
agreement  Is  seen  to  be  excellent. 

A second  similar  comparison  appears  in  Figure  9.  The  agree- 
ment here  again  Is  reasonably  good,  with  the  greatest  discrepancy 
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Mil.  Thrust  (Fig.  4) 
Max.A/B  (Fig.  5) 
Uss  7dB 


FIG.  6.  COMPARISON  OF  NOISE  CONTOURS  FOR  YJ-93  ENGINE  FOR 
MILITARY  POWER  AND  MAXIMUM  AFTERBURNER. 


Fie.  7.  COMPARISON  OF  NOISE  CONTOURS  OF  MODEL  JET  AT  MACH  NO.  1.25 


FIG.  9.  COMPARISON  OF  REFERENCE  SETS  OF  CONTOURS,  BASED  ON  FIGS.  27  AND 
30  Of  [6]. 


occurring  in  the  range  previously  noted,  l.e.,  for  X/D  values 
leas  than  7 in  the  fD  »*  8400  to  I68OO  Hz-inch  octave  band,  where 
shock  waves  in  the  Jet  may  affect  the  observed  noise.  Figure 
1.1.9  Indicates  that  two  sets  of  reference  contours,  rather  than 
a single  one,  would  provide  enhanced  accuracy. 

In  the  method  of  Plumblee  et  al  C6],  frequency  information 
is  restricted  to  the  three  octave  bands  centered  on  frequencies 
given  by  fD/c^  * 0.221,  0.442,  0.884.  Although  these  bands  may 
Include  the  apectrpl  peak  for  a given  condition  and  location, 
more  spectrum  information  may  be  of  Interest.  Thus,  a method  is 
needed  for  extrapolating  the  calculated  sound  levels  to  higher 
and  lower  frequencies. 

Perhaps  the  most  appropriate  method  for  this  purpose  is  that  of 
Thomson  [10],  which  specifies  four  non-dimensional  spectra 
for  the  nearfleld  region  0^X/D^30,  0 < Y/D  ^ 30.  The  spectrum 
shapes  proposed  by  Thccuson  are  selected  here  in  preference  to 
those  of  Cockburn  and  Jolly  [11],  because  the  former  repre- 
sent a larger  sample  of  e/igine  sizes  and  operating  conditions. 
Although  Cockburn  and  Jolly  provide  more  detail,  in  that  they 
provide  spectra  for  24  locations,  including  some  upstream  of  the 
nozzle  exit  plane,  their  spectra  are  based  on  only  one  set  of 
measurements  for  the  J57-P21  turbojet  engine  [12]. 

Although  the  nearfield  riolse  prediction  method  of  Thomson  is  more 
accurate  than  that  of  Cockburn  and  Jolly,  it  is  incomplete.  The 
most  useful  available  method  involves  calculation  of  the  octave 
band  sound  pressure  levels  at  three  center  frequencies,  follow- 
ing the  method  of  Plumblee  et  al^  and  then  using  the  appropriate 
spectrum  shape  from  Figure  10  to  extrapolate  the  calculated 
levels  to  higher  and  lower  frequencies.  (For  consistency,  the 
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spectrum  levels  of  Thomson  have  been  modified  to  provide  levels 
In  constant  percentage  bandwidths.) 

Itwould  bemore  satisfactory  to  extend  the  method  of  Plumblee 
0t  al  to  other  frequencies  by  extending  the  table  of  frequency- 
dependent  parameters;  however,  this  requires  considerably  more 
time  and  effort  than  available  for  the  present  study. 

2.2.2  Non-circular  nozzles 

Of  the  several  types  of  non-circular  exhaust  nozzles  treated  In 
the  literature,  only  two  have  any  present  practical  importance. 

If  the  various  noise  suppressor  designs  are  excluded.  These  are 
plug  and  slot  nozzles.  Conical  plugs  (Figure  11)  are  some- 
times used  .with  circular  nozzles,  and  slot  nozzles  are  being 
considered  for  powered-lift  systems,  such  as  the  augmentor  wing. 

The  available  acoustic  data  on  circular  plug  and  slot  nozzles 
are  limited  and  refer  to  the  acoustic  farfield.  This  data  has 

•V» 

been  reviewed  by  Stone  [20]  for  application  to  acoustic  fir- 
field  prediction  methods.  Since  no  better  information  is  avail- 
able, it  is  proposed  that  these  farfield  methods  be  applied  also 
for  nearfleld  predictions. 

Vlug  Nozzlee 

Measurements  of  the  farfield  noise  from  plug,  nozzles  indicate 
that  the  levels  are  slightly  lower  than  those  for  a circular 
nozzle  of  the  same  cross-sectional  area.  On  the  basis  of  limited 
experimental  evidence.  Stone  [20]  introduce  an  empirical 
correction  term 

- 3 log[0.10  + 2 (h/D)]  (6) 
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EQUIVALENT  DIAMETER,  2\/Dh-h 

HYDRAULIC  DIAMETER,  2h 

DIAMETER  RATIO, 

H 


FIG.  11.  PLUG  NOZZLE  GEOMETRY. 


where  h denotes  the  gap  height  and  D the  circular  nozzle 
diameter,  as  defined  In  Figure  11.  This  correction  term 
Is  to  be  added  to  the  overall  sound  pressure  level  calculated 
for  a circular  nozzle  of  the  same  cross-sectional  area. 


Ihe  available  published  data  do  not  Indicate  any  change  in  the 
directivity  patterns  in  presence  of  a plug,  but  do  show  some 
change  In  the  frequency  spectrum.  Stone  [20]  accounts  for 
this  frequency  change  by  Introducing  into  the  Strouhal  number 
a factor  (Dj^/D^)  ® , where  and  denote  the  hydraulic  and 

equivalent  diameters,  as  defined  in  Figure  11.  The  Strouhal 
number  for  a circular  plug  nozzle  then  becomes 


S = 


0 . If 


(7) 


In  the  absence  of  other  Information,  it  is  proposed  that  the  above 
modifications,  when  applied  to  the  prediction  procedur:^  of  Sec- 
tion 1.1.1  for  circular  nozzles,  be  adopted  for  nearfleld  noise 
predictions . 


Slot  nozzles 

There  exists  little  reliable  information  that  can  be  used  to  es- 
tablish quantitative  prediction  procedures,  and  such  data  as  is 
available  pertains  to  the  farfield,  where  some  of  the  non- 
axisymmetric  effects  may  be  less  pronounced  than  In  the  nearfield. 

The  data  analysis  of  Stone  [20]  for  nozzles  with  aspect  ra- 
tios of  ^1.8  and  69  indicates  that  in  the  acoustic ' farfield  the 
asymmetry  effects  are  small.  Thus,  Stone  neglects  these  effects 
and  proposes  a farfield  prediction  method  that  uses  the  relation- 
ships for  circular  nozzles. 


However,  Stoners  proposal  must  be  treated  with  caution  since, 
as  he  indicates,  all  measurements  made  at  = 90®  (l.e..  In 
the  plane  of^the  long  dimension  of  the  nozzle)  may  be  subject  to 
ground  reflection  effects.  However,  his  data  also  lead  to  other 
observations  which  need  to  be  considered: 

(a)  The  overall  sound  pressure  levels  at  ()>  = 90®  may  be 
1 to  2 dB  lower  than  the  corresponding  values  at 

^ “ 0®  (in  the  plane  of  the  short  dimension  of  the 
nozzle) . 

(b)  The  one-third  octave  band  spectrum  at  4)  *=  90®  peaks 
at  a frequency  which  is  about  a one-third  octave  be- 
low the  peak  frequency  at  <()  = 0° . 

(c)  Non-dimensional  spectrum  levels  at  high  and  low  fre- 
quencies are  lower  at  4>  = 90°  than  at  4>  = 0° . 

Similar  trends  may  also  be  observed  in  other  data  for  rectangu- 
lar nozzles.  In  [2]  and  [22],  the  peak  frequencies  at 
4»  = 90®  are  lower  than  those  at  (j)  = 0®  by  between  a one-third 
octave  and  a whole  octave.  The  frequency  spectrum,  when 
non-dimensionallzed  with  respect  to  the  overall  level,  at  low 
frequency  is  higher  at  <|)  = 90°  than  at  4»  = 0° . 

These  observations  support  a qualitative  descriptlc  i of  the 

noise  field,  which  associates  a decrease  in  high  frequency 
sound  pressure  level  at  90®  with  shielding  of  the  noise 
generated  in  the  initial  mixing  region  close  to  the  nozzle, 
in  a manner  similar  to  that  found  in  multitube  suppressor 
nozzles.  For  such  a phenomenon,  one  would  expect  to  find 
greater  noise  reduction  at  ♦ » 90®  in  the  nearfield  than  in  the 
farfield.  However,  this  i?  speculation,  and  for  the  present 
time  it  is  proposed  that  nearfield  noise  predictions  be  obtained 
following  the  recommendation  of  Stone  [20],  using  the  ap- 
proach for  axlsymmetric  circular  jets.  Thl^  approach  is  expected 
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to  give  conservative  results,  i.e.,  to  yield  predicted  levels 
that  will  overestimate  the  actual  levels  at  frequencies  greater 
than  the  peak  frequency  of  the  spectrum. 

2.2.3  Effective  location  of  jet  noise  sources 

In  many  prediction  methods,  it  is  assumed  that  the  noise 
sources  are  located  near  the  nozzle  exit  plane.  This  assumption 
is  of  little  consequence  In  the  prediction  of  farfield  sound  levels, 
but  is  important  in  the  estimation  of  nearfield  noise.  Thus,  it 
is  useful  here  to  revlev/  the  available  information  regarding  the 
effective  location  of  noise  sources  in  a Jet. 

The  experimental  techniques  that  have  been  used  to  measure  the 
location  of  noise  sources  in  a Jet  efflux  include  directional 
microphones,  reflectors.  Jets  discharging  through  a hole  in  a 
wall,  correlation  methods,  and  analysis  of  sound  pressure  level 
contours.  The  results  Indicate  that,  for  any  given  frequency 
band,  sources  are  distriuuted  over  a fairly  large  axial  distance. 
Some  of  this  observed  spatial  distribution  may  be  due  to  low 
resolution  of  the  measurement  techniques,  but  there  is  no  physi- 
cal reason  that  acoustic  energy  at  a given  frequency  be  generated 
at  a unique  location  along  the  Jet  axis. 

On  the  other  hand,  it  is  convenient  for  predicting  nearfleld 
noise  levels  to  assume  that  sound  at  a given  frequency  is  generated 
at  a unique  axial  distance  along  the  Jet,  the  distance  being  that 
associated  with  the  maximum  in  the  source  intensity  distribution 
curve  for  the  appropriate  frequency  band.  This  simplification 
is  used,  for  example,  by  Plumblee  et  al.  [6].  As  another  example, 
the  data  of  MacGregor  and  Slmcox  [233  show  that  the  source  dis- 
tribution half-pov/er  points  <i.e.,  locations  at  which  the  source 
power  is  3 dB  below  the  maximum)  occur  over  axial  distances  of 
to  12D,  centered  approximately  at  the  location  of  maximum 
sound  generation. 
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In  order  co  apply  the  experimental  results  for  source  location 
In  a noise  prediction  procedure,  one  needs  an  empirical  rela- 
tion between^ source  location,  nozzle  diameter,  Mach  number, 
temperature  and/or  other  parameters.  Investigators  have  non- 
dim:  '?ialized  source  location  with  respect  to  nozzle  diameter 

(F:./ue  12),  but  this  does  not  provide  a very  satisfactory 
collapse  of  the  data  for  different  flow  conditions.  The  data 
collapse  Is  Improved  if  nozzle  pressure  ratio  Is  introduced  as 
an  additional  scaling  parameter,  as  in  Figure  13,  but  even 
then  there  still  exists  considerable  variation  from  one  test 
condition  to  another.  Also,  the  method  does  not  take  into  ac- 
count variations  of  exhaust  temperature. 

It  might  be  argued  from  physical  reasons  that  jet  core  length 
is  a better  parameter  for  non-dlmenslonalizing  the  source  loca- 
tion X^,  since  the  core  length  includes  the  effects  of  pressure 
ratio  and  temperature  and  is  a more  meaningful  dimension  than 
jet  nozzle  diameter.  Potter  and  Crocker  [2^]  have  used  such 
non-dimenslonallzation  and  have  shown  that  this  reduced  the  data 
spread  considerably,  as  compared  with  an  X^/D  scaling;  however, 
no  comparison  is  available  with  scaling  in  relation  to  both 
nozzle  diameter  and  pressure  ratio. 

One  now  needs  to  consider  amethod  for  estimating  the  length  of 
the  potential  core.  Eldred  et  al,  [25]  have  presented  a re- 
lationship for  calculating  the  core  length  of  subsonic  and  super- 
sonic Jets.  They  contend  that  the  data  Indicate  no  influence 
of  temperature  between  300®F  and  2000°F,  and  therefore  include  no 
temperature  effect  In  their  empirical  relation, 

ij.  - 3.‘»5D[1  + 0.38M^3^  (8' 

where  M4  is  the  exit  Mach  number  of  the  jet. 
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D PR 


FOR  LEGEND  SEE  FIG  12. 


FIG.  13.  JET  SOURCE  LOCATION,  NONDIMENSIONALIZED  WITH 
RESPECT  TO  NOZZLE  DIAMETER  AND  PRESSURE  RATIO. 
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In  contrast,  Abramovich  [26]  sho»fs  a strong  dependence  of 
potential  core  length  on  Jet  exhaust  temperature.  Figure  1^1 
contains  a family  of  curves  relating  core  length  to  flow  Mach 
number  for  several  total  temperatures.  This  figure  also  In- 
cludes a curve  corresponding  to  Equation  (8),  from  which  one  may 
see  that  the  two  methods  predict  similar  core  lengths  for  ambient 
temperature  Jets,  but  not  for  hot  Jets.  The  curves  of  Figure  1^1 
also  disagree  with  the  data  of  MacGregor  and  Slmcox  [23]  for  Jets 
with  a temperature  ratio  Tj/T^  » 3.1,  which  data,  as  evident 
from  the  figure,  fall  near  the  values  predicted  for  ambient 
temperature  Jets. 


Figure  15  shows  the  source  location  data  of  Figure  13  plotted 

In  terms  of  where  Is  the  core  length  calculated  from 

o c c 

Abramovich’s  method  (Figure  1^).  It  Is  seen  that  there  Is  little 
or  no  Improvement  In  the  data  collapse,  as  compared  to  Figure  15. 
If  core  lengths  are  calculated  on  the  basis  of  ambient  temperature 
Jets  (by  the  method  of  Eldred  et  al,),  the  data  collapse  Is 
slightly  better  than  that  shown  In  Figure  15.  Thus,  the  question 
of  optimum  scaling  Is  unresolved. 


The  relationship  proposed  by  Plumblee  et  al.  [6]  does  not 
lend  itself  readily  to  plotting  In  the  forms  of  Figure  13  or  15. 
However,  source  distributions  can  be  calculated  if  specific 
assumptions  are  made.  The  source  location  is  given  by 


> 1 1» 


VITO)/ 


(9) 


where  coefficients  ,,,  and  are  given  in 

[6]  for  the  three  octave  bands  centered  at  Strouhal  numbers 
fp/c©  ■ 0,221,  O.^^I2,  and  0.88^.  Under  the  assumption  that 
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FI6.  14.  jet  core  length  as  FUNCTION  OF  JET  TEMPERATURE 
AND  MACH  NUMBER. 
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FOR  LEGEND  SEE  FIG.  12. 

FIG.  15.  JET  NOISE  SOURCE  LOCATION*  NONDIMENSIONALI ZED 
MITH  RESPECT  TO  CORE  LENGTH. 
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the  nozzle  pressure  ratio  Is  2.0  and  the  Jet  temperature  is 
either  519°  (ambient  Jet)  or  1600°R  (typical  engine  condition), 
one  obtains  ^;he  distributions  plotted  In  Figure  l6.  These 
values  are  consistent  with  the  experimental  data  of  Figure 
13. 


Apparently  there  Is  considerable  uncertainty  in  the  Interpreta- 
tion of  available  Information  regarding  the  location  of  noise 
sources  In  a Jet.  The  equation  derived  by  Plumblee  et  al  pre- 
dicts source  locations  which  are  within  the  range  of  measured 
values,  but  the  equation  Is  available  for  only  three  frequen- 
cies. A method  which  can  be  used  over  a wider  range  of  frequen- 
cies Is  required.  Since  the  scatter  associated  with  each  of 
three  non-dlmenslonallzlng  procedures,  - that  Is  In  terms  of 

^0.1  ^o 

IT  m * r(hot  fibwr  ocotttowT* 

c c 

Is  approximately  the  same,  no  method  can  be  selected  on  the 

X 

basis  of  superior  performance.  Thus,  the  scaling  factor  _jo.  ^ 

D PR 

may  be  selected  on  the  basis  of  convenience.  With  this  non- 
dimensional  relationship,  a representative  curve  has  been  con- 
structed that  lies  approximately  In  the  center  of  the  experimen- 
tal range  of  data  and  also  Is  consistent  with  the  predictions  of 
Plumblee  et  at.  This  proposed  model  distribution  Is  shown  In 
Figure  16. 
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STROUHAL  NUMBSl  fDAj 


POINTS  PREDICTED  BY  V»ETHOD  OF  PLUMBLEE  ET  AL  FIG.  6 . 

o T|--1600°R  A T^=520®R 


FIG.  16.  PROPOSED  MODEL  FOR  NOISE  SOURCE  LOCATION  IN  JET 
FLOW. 
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2.3  Afterburner  Noise 

Combustion  noise  In  a Jet  engine  Is  generated  In  the  com- 
bustion chambers  upstream  of  the  turbines  and  In  the  after- 
burner, If  such  a device  Is  In  operation.  Combustion  chamber 
noise  Is  Important  only  if  the  Jet  exhaust  velocities  are  low. 
Afterburner  noise  is  generated  when  fuel  is  burned  in  the  Jet 
tailpipe,  downstream  of  the  turbine,  resulting  In  high-speed  Jet 
flow.  However,  the  associated  Jet  Mach  number  may  be  no  higher 
than  that  for  the  non-afterburning  case. 

2.3.1  Available  data  and  early  prediction  methods 

Several  measurements  of  afterburner  noise  have  been  made, 
some  in  the  far  field  [32,  33]  and  others  in  the  near  field  [12, 
15].  Although  no  models  have  been  proposed  to  describe  the 
afterburner  noise  characteristics,  one  nearfleld  study  [15] 
compared  the  measured  sound  pressure  levels  with  those  predicted 
by  the  method  of  Plumblee  et  al.  [6]. 

An  early  prediction  method  for  the  nearfield  noise  of  a Jet 
engine  with  afterburner,  proposed  by  Pranken  et  al , [3^1],  con- 
sists of  three  steps; 

(1)  Calculation  of  overall  sound  pressure  level  for 
military  power; 

(2)  Addition  of  a correction  term  for  afterburner  noise; 
and 

v" 

(3)  Shifting  of  the  sound  pressure  level  spectrum,  based 
on  non-afterburner  operation,  down  in  frequency  by 
half  an  octave. 

The  method  makes  no  allowance  for  different  levels  of  afterburner 
operation,  and  does  not  account  for  directivity. 


2.3.2  Proposed  nearfield  noise  prediction  method 

In  order  to  develop  an  alternative  approach,  one  may 
assume  the  main  component  of  afterburner  noise  to  be  due  to  the 
Increased  flow  velocity  and  nozzle  diameter,  so  that  one  may  apply  the 
method  of  Plumblee.  In  [13]>  application  of  maximum  after- 
burner thrust  was  reported  to  produce  the  following  changes  in 
Jet  characteristics: 


Military 

Power 

Maximum 

Afterburner 

Change 

Total  Temperature  (°R) 

2066 

3651 

^77^r 

Exit  Velocity  (ft/sec) 

2466 

3246 

+ 32% 

Exit  Area  (ft^) 

6.4 

8.2 

+28? 

Exit  Mach  Number 

1.28 

1.26 

-2% 

The  corresponding  comparison  between  theory  and  experiment  shows 
agreement  here  is  as  good  as  that  for  the  non-afterburner  case. 

For  example.  Figure  17  shows  the  difference  between  measured 
and  predicted  o\  -^rall  sound  pressure  levels  at  all  microphone 
locations  when  ground  reflection  is  taken  into  account.  The  mean 
difference  is  1.8  dB  for  the  maximum  afterburner  stage  and  2.6  dB 
for  military  power.  The  corresponding  standard  deviations  are 
2.4  d3  and  2.1  dB. 

Comparisons  in  terms  of  octave  band  sound  pressure  level 
contours  are  shown  in  Figures  I8  and  19,  for  center  fre- 
quencies of  approximately  80  and  I60  Hz,  respectively.  These 
are  the  lower  two  frequency  bands  available  from  the  Plumblee  et  al 
method.  Just  as  for  non-afterburner  noise,  one  finds  here  that 
the  agreement  between  measured  and  predicted  levels  Is  quite 
good  at  the  higher  frequency,  but  that  the  agreement  is  poor  at 
the  lower  frequency. 


The  method  of  Plumblee  et  al,  [6]  has  the  disadvantage 
that  noise  due  to  the  combustion  process  Is  not  Included  in  the 
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FIG.  17.  COMPARISON  OF  MEASURED  AND  PREDICTED  OVERALL  SOUND 
PRESSURE  LEVELS  FOR  AFTERBURNER  OPERATION  [15], 
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prediction  procedure  (unless  some  combustion  noise  was  present. 
Inadvertently,  in  the  basic  data  used  to  construct  the  analyt- 
ical model).  At  the  present  time,  there  is  available  no  method 
for  predicting  the  combustion  noise,  or  even  for  determining 
whether  it  Is  a significant  contributor  to  afterburner  noise. 
Thus,  further  studies  are  required,  but  at  present  the  method 
of  Plumblee  may  be  used  as  an  adequate  approximate  tool. 
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2.4  Coaxial  Jets 


Turbofan  engines  with  bypass  flow  have  coaxial  nozzles 
thi'ough  which  pass  flows  that  usually  have  different  tempera- 
tures and  velocities.  Depending  on  the  particular  engine  de- 
sign, the  nozzles  may  be  coplanar,  the  primary  nozzle  may  be 
retracted  (as  in  the  JT8D)  or  the  fan  duct  or  secondary  nozzle 
may  be  shortened  (as  in  the  JT9D).  A schematic  representation 
of  a coaxial  nozzle,  showing  the  various  noise  source  regions, 
is  presented  In  Figure  20. 

2.4.1  Methods  for  prediction  of  farfield  noise 

Several  experimental  investigations  of  the  farfield  noise 
of  coaxial  Jets  have  been  conducted  on  model  and  full-scale 
engines.  The  conclusions  reached  by  the  investigators  differ, 
and  several  empirical  relationships  have  been  proposed  to  des- 
cribe the  radiated  sound  field.  The  main  point  of  contention 
is  whether  the  primary  Jet  velocity  or  the  secondary  velo- 
city Vj2  should  be  used  as  the  principal  parameter  of  the  noise 
field.  No  attention  appears  to  have  been  given  to  the  nearfield 

Prom  a physical  standpoint,  it  is  logical  to  consider  the 
noise  generation  in  two  parts:  (1)  the  interaction  between  the 

primary  and  secondary  flows,  and  (2)  the  Interaction  between  the 
secondary  flow  and  the  surrounding  air.  With  this  approach  the 
effect  of  forward  motion  is  assumed  to  be  restricted  to  the 
secondary  flow  region,  as  discussed  in  Section  1.5. 

Both  of  the  available  prediction  methods  for  the  acoustical 
farfield  of  coaxial  Jets  are  empirical  and  based  mainly  on 
model-scale,  cold-jet  experiments,  and  both  are  based  on  predic 
tlon  methods  for  single  Jets,  with  suitable  modifications  added 
to  account  for  the  dual  flow. 
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FIG.  20,  SCHEMATIC  SKETCH  OF  NOISE  SOURCE  REGIONS  FOR 
CO-ANNULAR  JETS. 


55 


stone  [20]  calculated  the  farfleld  OASPL^^  at 
e ■ yo*  for  the  primary  Jet  alona,  c.nd  then  modifies  the  result 
to  account  for  the  secondary  flow,  according:  to  the  following 
equation: 

OASPLgg  - OASPL^q^^  + 5 log  CT^j/Tj23 


(1C) 


The  suffixes  1,  2 here  refer  to  the  primary  and  secondary  flows, 
respectively.  The  effect  of  Jet  temperature  is  Included  rather 
arbitrarily,  since  it  Is  not  directly  associated  with  the  appro- 
priate flow  area  and  Jet  velocity.  (An  improvement  is  obviously 
required  here.)  Adjustments  can  be  made  to  Equation  (10)  to 
give  the  overall  sound  pressure  level  at  any  required  angle  0, 
and  one-third  octave  band  spectra  can  be  obtained  by  using  typi- 
cal shapes  for  single  Jets. 

Dunn  and  Peart  [35],  on  the  other  hand,  calculate  the 
sound  pressure  levels  associated  with  the  primary  and  secondary 
flows  separately.  (The  calculation  for  the  secondary  flor  can 
be  applied  also  to  the  mixed  flow  region,  when  the  primary  and 
secondary  flows  are  mixed.)  For  the  primary  Jet,  they  calculate 
the  one-third  octave  band  spectrum  on  the  basis  of  the  area 
and  Jet  velocity  of  the  primary  Jet.  They  then  account  for 
the  presence  of  the  secondary  flow  by  introducing  a relative 
velocity  correction  term  (f),  such  that 

» 5Fi:(f)  + AjgCf)  (11) 
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where  Sj?L(f ) Is  the  space-averaged  one-third  octave  band  sound 
pressure  level  calculated  directly  from  single  Jet  analysis, 
(SPL(f))j^  is  the  corrected  value  for  the  primary  flow  and 
^ijCf)  is  the  correction  term,  where 


10  log 


T - V 

hi 


'ji 


m 


(12) 


The  index  m = m (A^/A^^,  f/f^)  is  given  by  Figure  21. 


For  prediction  of  the  noise  field  from  the  secondary  flow,  they 
use  the  prediction  method  for  a single  jet  directly,  but  with 
the  density,  area  and  Jet  velocity  parameters  defined  as  follows: 


A.  »= 


V,  = 


= flow  density  of  secondary  flow 
A,  + Ar, 


^l^n  ^ ^2^J2 


1/2 


(13) 


It  should  be  noted  that  Equation  (13)  was  taken  by  Dunn  and 
Peart  from  the  work  of  Bielak  [36].  His  data  were  measured 
on  model-scale  Jets,  where  the  flow  from  the  primary  and  second- 
ary nozzles  was  at  ambient  temperature,  and  he  developed  the 
equation  for  the  effective  Jet  velocity  on  the  basis  of  con- 
tinuity of  momentum  for  ambient  temperature  flow.  Thus,  it 
appears  appropriate  to  introduce  new  equations  which  take  into 
account  hot  primary  flows,  such  as  occur  in  a typical  turbofan 
engine.  If  one  again  uses  continuity  of  momentum,  one  finds 
that  the  effective  Jet  velocity  obeys 


^l^l^Jl  Pp^2^J2 
^l^l^Jl  * ^2*2^J2 


(m) 


A similar  approach  could  be  followed  for  equivalent  density,  if 
required. 
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FIG.  21.  INDEX  m AS  FUNCTION  OF  STROUHAL  NUMBER  f/fj  [36]. 


2.4.2  Suggested  method  for  prediction  of  nearfield  noise 

Dunn  and  Peart  C35]  have  developed  their  own  method  for 
predicting  farfield  noise  of  single  jets,  but  have  not  dealt 
with  nearfield  sound.  However,  since  their  general  approach  is 
more  readily  adaptable  to  the  prediction  procedure  developed  by 
Plumbiee  et  al  (Section  1.1)  than  is  the  procedure  of  Stone,  it 
is  proposed  that  the  procedure  of  Dunn  and  Peart  be  used  in 
setting  up  a preliminary  method  for  predicting  nearfield  sound 
levels  of  a coaxial  Jet.  Of  course,  comparisons  between  theory 
and  experiment  are  required  to  Justify  this  approach. 


The  nearfield  Jet  noise  procedure  of  Plumbiee  et  al.  [6], 
which  was  extended  to  additional  frequencies  in  the  present  re- 
port, requires  as  input  parameters  the  Jet  diameter  D,  exit  Mach 
number  M,  Jet  total  temperature  T,  and  Jet  noise  source  location 
x^.  These  parameters  here  need  to  be  established  for  the  pri- 
mary and  secondary  flows.  For  the  primary  flow,  the  choice  of 
parameters  is  readily  apparent;  viz.,  D = D^,  M » T = T^, 
and  the  source  location  x^  may  be  calculated  as  for  the  single 
Jet.  Calculation  of  the  correction  term  •’’s  also  straight- 

forward . 


For  the  secondary  flow,  the  effective  diameter  can  be 
obtained  from  the  effective  area  A^,  Equation  (13) > and  the 
effective  Jet  velocity  Vj  can  be  calculated  from  Equation 
(1^).  An  approximation  to  the  effective  total  temperature 
T^  is  given  by 


T^  + BTg 


‘e 


(15) 


where  $ is  the  engine  bypass  ratio.  (Although  Equation  (15) 

Is  an  engineering  approximation,  rather  than  an  exact  relation. 


Its  accuracy  should  be  adequate  for  nearfleld  noise  estimation.*) 

The  effective  Jet  Mach  number  M can  be  estimated  on  the  basis 

e 

of  the  above  -values  of  effective  Jet  velocity  and  total  tem- 
perature . 

The  total  sound  pressure  level  SPL  in  a given  frequency 
band  then  may  be  determined  from 

SPL  = 10  log  (16) 

where  (SPL)^  and  (SPL)2  are  the  sound  pressure  levels  calculated 
for  the  primary  and  secondary  flows,  respectively,  and  where  the 
two  noise  sources  are  assumed  to  be  incoherent. 

The  method  of  Plumblee  et  al  may  be  less  accurate  for  the 
prediction  of  coaxial  flow  noise,  because  of  the  low  velocities 
associated  with  secondary  flow  in  turbofan  engines.  Table  2 
lists  typical  fan  Jet  velocities  as  being  in  the  range  of  700 
to  1200  ft/sec,  whereas  the  Plumblee  method  is  based  on  measure- 
ments in  the  range  of  1200  to  3000  ft/sec. 


• The  definitions  of  A , V.  and  T are  not  mutually  consistent 
in  terms  of  mass  floS.  ’'(This  Is  also  true  of  Equation  (13)). 
However,  the  proposed  definitions  provided  reasonable  approxi- 
mations for  nearfleld  noise  estimation,  because  much  of  the 
Jet  noise  is  generated  downstream  of  the  Jet  core,  where  the 
differences  between  total  and  static  temperatures  are  less 
than  at  the  noezle. 
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2.5  Shock  Noise 


Most  investigations  of  noise  from  supersonic  Jet  exhausts 
have  dealt  with  shock-free  flow.  However,  additional  noise 
sources  arise  from  the  presence  of  shock  waves  in  flow  which  is 
not  fully  expanded.  This  shock  noise  can  be  of  broadband  or 
discrete-frequency  nature;  broadband  noise  is  caused  by  inter- 
action between  turbulence  and  the  shockwave,  and  discrete  fre- 
quency noise  is  generated  by  a feedback  mechanism  between  the 
shock  fronts  and  Jet  nozzle. 

Shocks  occur  in  any  supersonic  flow  from  a convergent 
nozzle,  and  in  supersonic  flow  from  a convergent-divergent 
nozzle  operating  away  from  its  design  point.  Shocked  flow  can 
be  classified  as  either  overexpanded  or  underexpanded.  In  the 
overexpanded  flow  case,  the  nozzle  static  pressure  p^  is  below 
ambient  (p^<p„),  whereas  in  underexpanded  flow,  p^>p«.  For  fully 
expanded  flow,  p^  » p^.  Typical  flow  regimes  are  illustrated 
in  Figure  22. 

2.5.1  Shock  noise  vs  shock-free  noise  ^ 

Although  broadband  and  discrete- frequency  shock  noise  has 
been  measured  by  several  investigators,  very  little  effort  has 
been  directed  towards  establishing  prediction  procedures,  even 
for  the  farfield  case.  Indeed,  most  authors,  including  Chenefa 
[2]  specifically  exclude  shock  noise  from  their  prediction 
methods  for  supersonic  Jet  noise,  because  of  the  complexity  of 
the  problem.  Stone  [20]  proposes  a method  for  calculating 
farfield,  broadband  shock  noise,  but  excludes  discrete  frequency 
noise  because  *'such  tones  can  often  be  eliminated  by  slight  de- 
sign moQlflcatlons''. 

Prom  the  description  of  the  development  of  the  predltion 
method  of  Plumblee  at  [6]  it  Ic  not  clear  whether  shock 


61 


(a)  Ov«rex|xind«cl.  Nozzle  Static  Pressure  P^<Pg 


(b)  Fully  dixpanded.  P^«  P^ 

rf"~ 

J 

(c)  Underexpanded . P^  > P^ 


Expansion  Rafllection 


FIG.  22.  TYPICAL  SHOCK  PATTERNS  IN  OVER-  AND 
UNDER-EXPANDED  JETS. 
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noise  was  Included.  However,  It  Is  possible  that  some  of  the 
dlsci-epancles  between  predicted  and  measured  levels  are  due  to 
the  presence  of  shock  noise  In  the  experimental  data.  It  Is 
assumed  for  ohe  present  that  the  Plumblee  model  does  not  take 
shock  noise  Into  account,  and  that  a method  for  predicting  this 
noise  needs  to  be  established. 

The  contribution  made  by  shock  noise  can  be  determined  by 
comparing  noise  data  for  a convergent-divergent  nozzle  operating 
at  or-design  conditions  with  data  on  the  same  nozzle  at  off- 
design  conditions,  or  by  comparing  data  for  a convergent- 
divergent  nozzle  at  on-deslgn  conditions  with  similar  data  for 
a convergent  noz^^le  o].  rfi.ing  at  the  same  nozzle  pressure  ratio. 
Such  a comparison  for  i i Leld  noise  has  been  made  by  Slncox 
[37];  sample  data  <ire  wn  In  Figure  23.  This  figure 
compares  data  for  a con  orgent-divergent  nozzle  designed  for 
M » 1.^1  operation  with  msasui*ements  for  a convergent  nozzle. 

At  a flow  temperature  of  520°R,  Figure  23a  Indicates  a de- 
crease In  radiated  sound  in  the  neighborhood  of  M “ 1.^  for  the 
convergent-divergent  nozzle,  but  no  corresponding  drop  for  the 
convergent  nozzle.  At  the  higher  flow  temperature  of  2^60®R, 
no  decrease  In  sound  level  is  observed  at  the  on-design  condition, 
and  the  noise  levels  for  convergent  and  convergent-divergent 
nozzles  are  nearly  the  same  (Figure  23b).  Simcox  claims 
that  this  phenomenon  occurs  because  at  a given  nozzle  pressure 
ratio  the  shcck-free  noise  increases  much  more  rapidly  with  Jet 
temperature  than  does  the  shock  noise.  (In  fact,  the  data  sug- 
gest that  shock  noise  Is  independent  of  temperature.) 
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2.5.2  Broadband  Shock  Noise 

An  Interim  prediction  method  for  f.irfleld  noise  generated 
by  shock  wave  effects  In  supersonic  flow  has  been  proposed  by 
Stone  [20].  The  method  assumes  that  overall  and  spectral 
sound  levels  are  non-dlrectlonal  and  Independent  of  Jet  temper- 
ature. According  to  this  method,  the  overall  sound  pressure 
level  Is  given  by 

OASPL,  - H7  + 10  log  [(Aj/r®)  (Po'o^PiSA'iSa’^] 


+ 10 

log 

12.5 

(Mj  - 1)3 

- 10 

log  j 

[o.04 

+ (M^  - 1) 

- 10 

log 

[- 

2 (Mj  - 1)^ 

] 

- 10 

log  1 

* 

(0.05/Mj^) 

- 10 

log 

(0.0086/Mj^ 

where  M^>  1,  One-third  octave  band  spectral  levels  are  obtained 
from  the  non-dlmenslonal  spectrum  shown  In  Pigure*24.  The 
shock-noise  spectrum  levels  are  then  added  to  corresponding  lev- 
els for  shick-free  noise. 

This  approach  agrees  with  the  experimental  finding  of 
Slmcox  f i7 1 that  t he  shook  nolr.c  Is  Independent  of  tempera- 
ture, but  falls  to  predict  the  nearfleld  shock  noise  levels 
measured  for  a model-scale  heated  Jet  by  Knott  et  al  [38] 

A typical  exai^>le,  for  a nossle  pressure  ratio  of  3*067  end  Jet 
total  temperature  of  2220*R,  is  snown  In  Figure  ?'i.  Tf  the 
Irregularities  associated  with  discrete  frequency  shock  noise 
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are  excluded,  the  estimated  sound  levels  turn  out  to  be  about 
16  dB  higher  than  the  measured  values.  The  predicted  spectrum 
shapes, on  the  other  hand,  agree  reasonably  well  with  the  ob- 
served ones. 

Two  alternatives  present  themselves.  One  is  the  modifica- 
tion of  Equation  17  to  bring  the  predicted  and  measured  levels 
into  agreement.  This  can  be  done  for  the  data  in  Figure  24 
by  changing  the  constant  in  Equation  17  from  147  to  131. 

However,  there  still  remains  the  question  as  to  whether  the 
(l/r  ) distance  law  Is  adequate  for  the  near  field;  If  it  is 
not,  the  constant  131  may  not  be  applicable  at  other  locations 
in  the  near  field. 

The  second  alternative  is  co  use  available  experimental 
data  to  determine  Incremental  changes  in  spectrum  levels  due  to 
the  presence  of  shock  noise.  These  incremental  changes  could 
then  be  added  to  predicted  spectra  for  shock-free  noise.  That 
Is,  If  the  Incremental  change  due  to  broadband  shock  noise  is 
AgB  total  spectrum  level  is  SPL^  « SPL  (shock  free) 

^ ^SB 

Experimental  results  which  can  be  used  to  determine  values 
of  Agg  are  found  In  [38  and  373.  If  one  takes  for  each 
test  condition  the  Increases  In  sound  pressure  level  due  to  the 
presence  of  shock  noise.  If  one  averages  these  Increases  over 
all  frequencies  above  the  shock  fundamental  discrete  frequency  f^, 
and  If  one  plots  the  resulting  values  as  a function  of  Jet  total 
tMq>erature,  one  obtains  the  plot  shown  In  Figure  25.  One  may 
then  draw  a curve  to  represent  the  typical  relation  between  the 
sound  level  Increment  and  jet  total  teRq>erature. 


JET  TOTAL  TEMPERATURE  Ty  ®R 


FIG.  25.  SOUND  LEVEL  INCREMENT  DUE  TO  BROADBAND  SHOCK  NOISE, 
AS  A FUNCTION  OF  JET  TEMPERATURE. 
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On  the  basis  of  the  experimental  evidence  given  In 
[37]  and  [38],  one  may  Introduce  the  following  reasonable 
approximations : 

(1)  Shock  noise  makes  a negligible  contribution  to  Jet 
nearfleld  noise  at  frequencies  below  the  shock  fundamental 
discrete  frequency  f^^  (this  frequency  is  discussed  in  Section 
1.^1. 3). 

(2)  At  frequencies  above  f^,  the  presence  of  broadband 
shock  noise  increases  the  Jet  noise  spectrum  levels  by  an  in- 
crement that  is  independent  of  frequency,  but  dependent  on 
Jet  temperature. 

(3)  The  incremental  increase  in  sound  level  is  independent 

of  Jet  Mach  number  when  Below  M = l.^»  the  increment  de- 

creases with  Mach  number  to  a value  of  zero  at  M = 1.0. 

(i|)  An  exception  to  (3)  occurs  in  the  neighborhood  of  the 
design  pressure  ratio  for  a convergent-divergent  nozzle.  In 
this  case,  the  increment  is  zero  if  the  pressure  ratio  is  be- 
tween 0.9  and  1.1  times  the  design  pressure  ratio. 

(5)  The  increment  is  applied  at  all  measurement  locations 
where  ^>*15*^.  For  the  increment  is  zr^o. 

By  means  of  these  simplifying  approximations,  one  may  transform 
the  empirical  curve  in  Figures  25  into  the  family  of  predic- 
tion curves  shown  in  Figure  26. 

2.5.3  Di screte-frequency  shock  noise 

Discrete-frequency  shock  noise  has  been  studied  extensively 
in  the  laboratory,  by  use  of  carefully  constructed  model-scale 
nozxles.  Under  these  conditions,  high  levels  of  discrete- 
frequency  shock  noise  were  usually  observed  [37,  39,  ^0,  ^11], 
but, In  some  cases  the  tones  were  low  In  level  or  were  not 
observed  at  all  [38,  »Dn  full-scale  Jet  engines,  discrete- 

frequency  shock  noise  has  been  encountered  Infrequently  [42,^33. 
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FIG.  26.  PROPOSED  INTERIM  MODEL  FOR  PREDICTION 
OF  BROADBAND  SHOCK  . '^ISE  INCREMENT. 
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Two  explanations  have  been  suggested  for  the  variability  in 
generation  of  discrete-frequency  shock  noise.  One  explanation 
relates  this  effect  to  the  jurbulence  in  the  flow  upstream  of  the 
nozzle;  Hewlett  [42]  ooserved  that  the  level  c discrete- 
frequency  shock  noise  decreased  as  the  turbulence  level  increased. 
The  second  explanation  ascribes  thif;  effect  to  Irregularities  in 
the  construction  of  the  nozzle;  Powell  [44]  demonstrated  that 
a small  notch  at  the  nozzle  exit  of  a model-scale  Jet  at  ambient 
temperature  would  reduce  the  intensity  of  the  discrete-frequency 
shock  noise  or  eliminate  it  completely.  Both  of  these  explana- 
tions could  account  for  the  apparent  absence  of  this  form  of 
shock  noise  in  the  acoustic  fields  of  many  Jet  engines,  even 
though  they  may  operate  at  high  nozzle  pressure  ratios.  It  is 
also  reported  [42]  that  the  discrete-frequency  noise  of  the 
R-R  Spey  engines  on  the  HS  121  airplane  disappeared  when  a 6- 
chute  suppressor  nozzle  was  fitted  to  the  Jet  pipe. 

There  is  available  no  criterion  for  the  onset  of  discrete- 
frequency  shock  noise;  this  shortcoming  is  a major  obstacle  to 
establisJjment  of  a method  for  predicting  nearf  1 eld  sound  levels. 
Since  experimental  data  necessary  to  construct  such  a criterion 
are  not  available,  an  alternative  approach  is  suggested  in  the 
following  discussion.  This  approach  uses  sound  field  information 
from  model-scale  tests  which  exhibited  strong  discrete-frequency 
shock  noise  and  makes  the  assumption  that  these  sound  levels  pro- 
vide an  upper  bound  to  the  levels  likely  to  be  encountered  on 
full  scale  engines.  (Turbulence  and  nozzle  irregularities  would 
normally  reduce  the  full-scale  engine  discrete-frequency  shock 
noise  levels.) 
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The  original  work  of  Powell  [19]  on  ambient  temperature 
Jets  established  an  equation  for  the  fundamental  frequency  of 
dlscrete-freqpency  shock  noise.  This  equation  is 


f « 
1 


K DXR-R^]^ 
1 j c 


(18) 


where 


and 


Dj  “nozzle  diameter 

R “ nozzle  pressure  ratio 

= critical  nozzle  pressure  ratio  = 1.893 

K “ constant. 

1 


Powell  obtained  a value  K = 3.0.  but  subsequent  experimental 

1 

studies  by  other  Investigators  resulted  in  lower  values  of  K^, 


These  values  are  shown 

in  the  following  table: 

Reference 

Test  Condition 

Rowlett  [^*2] 

Cold  model  Jet 

2.i\ 

Slmcox  [373 

Cold  model  Jet 

2.7 

Westley  [^1] 
and  Woolley 

Cold  model  Jet 

1.0-1. 7 

Knott  et  al  [38] 

Heated  model  Jet 

1. 7-2.2 

Hay  [il3] 

Pull  scale  engine 

1.1 

Since  the  experimental  values  of  are  fairly  uniformly  distribu- 
ted in  the  range  between  1.0  and  3«0,  a value  Kj  » 2.0  will  be 
used  for  the  present  model.  The  frequency  fa  of  the  second  har- 
monic of  the  shock  noise  is  given  by 


f - 2f 

t 1 


(19) 
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In  practice,  the  situation  Is  rather  more  complicated  than  Is 
Implied  by  the  above  two  frequency  equations.  For  example, 
WestJ.ey  and  Woolley  [^0]  have  shown  that  several  families  of 
shock  tones  can  be  present  and  the  noise  field  may  change  from 
one  dominant  family  to  another.  Howev  the  simplified  ap- 
proach is  adequate  for  present  purposes,  if  it  is  acknowledged 
that  the  error  in  predicted  frequency  may  be  as  large  as  +50%, 

The  detailed  data  on  nearfleld  acoustic  pressures  for  the 
fundamental  of  the  discrete-frequency  shock  noise,  which  have 
been  amassed  by  Westley  and  Woolley  C^O,  41]^  will  be 
used  here  to  locate  the  effective  source,  construct  simplified 
noise  contours  and  determine  upper  bounds  for  sound  levels. 
Because  the  data  include  very  little  information  for  the  second 
harmonic,  the  proposed  model  for  this  component  will  rely  on 
engineering  Judgement.  Obviously,  experimental  verification  of 
the  proposed  model  is  required. 

Westley  and  Woolley  [^1]  observe  that  the  maximum  sound 
levels  for  the  dominant  fundamentaD  frequency  occur  at  a dis- 
tance from  the  nozzle  which  increases  with  nozzle  pressure  ratio. 
E.g.,  for  a nozzle  pressure  ratio  of  2.67>  the  maximum  levels 
occur  in  the  neighborhood  of  the  third  and  fourth  shock  cells, 
whereao  at  a pressure  ratio  of  5* 67,  the  maxima  occur  at  the 
fifth  and  sixth  cells.  For  present  engines  typical  nozzle  pres- 
sure ratios  lie  below  about  3*0  (see  Table  2).  Hence,  the 
proposed  model  assumes  that  the  effective  location  of  the  sources 
of  discrete-frequency  shock  noise  (fundamental  and  all  harmonics) 
Is  approximately  3«3  shock  cell  lengths  downstream  of  txhe  nozzle 
exit  plane,  and  lies  on  the  Jet  centerline. 
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The  characteristic  length  s of  a shock  cell  can  be  taken 
appro ximatelj'  to  obey 

' § - [R-R,]‘‘  , (20 > 

and  thus  the  effective  source  location  x may  be  approximated 

s 

by 


= 3.5  CR-R^]’’  (21) 

J 

Theoretical  directivity  patterns  constructed  by  Powell 
[9]  for  specific  operating  conditions  show  a series  of  lobes 
(Figure  27a);  for  the  fundamental,  the  main  lobe  is  in  the 
upstream  direction;  for  the  second  harmonic,  the  main  lobe  is  at 
right  angles  to  the  Jet  axis.  The  nearfield  measurements  of 
Westley  and  Woolley  [AO,  i*l]  show  an  extremely  complicated 
pattern  of  sound  pressure  level  contours  for  the  fundamental 
frequency  component,  as  Illustrated  in  Figure  27b.  Because  the 
actual  contours  are  too  complicated  for  use  in  an  engineering  , 
prediction  model,  series  of  smoothed  contours  have  been  con- 
structed to  approximate  the  experimental  data,  as  also  shown  in 
the  aforementioned  figure. 

No  nearfield  data  are  available  for  the  second  harmonic, 
but  circular  contours  as  shown  in  Figure  27c  (with  oontorr. 
on  the  normal  to  the  Jet  axis  at  the  effective  source  location, 
and  tangent  to  each  other  at  the  effective  source  location)  are 
proposed  as  a simplified  representation  of  the  contours  of 
Powell. 

In  order  tr  predict  upper  bounds  on  the  sound  pressure 
level,  one  may  refer  to  data  taken  along  the  normal  to  the  Jet 
axis  at  the  effective  source  location,  which  data  f'lO,  'll] 
Indicates  a sound  pressure  dependence  on  the  fourth  power  of 
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distance.  Although  from  a physical  standpoint  It  might  be  more 
satisfactory  to  have  a model  With  an  origin  located  on  the  Jet 
boundary,  such  a model  results  In  a more  complex  dependence  on 
distance,  and'the  simpler  method  is  adopted  at  present  for  the 
sake  of  simplicity. 

If  one  assumes  that  the  acoustic  power  generated  by  the 
dlscrete~frequency  shock  noise  is  directly  proportional  to  the 
nozzle  area  Aj , one  finds  that  one  may  write  a reduced  sound 
pressure  level  as 

♦ - SPL^^  -10  log  (Aj/Apgf)  + 40  log  (y/4D)  (22) 

where  the  reference  area  « 0.028  ft*,  corresponding  to  the 

model  tests  of  Westley  and  Woolley.  Application  of  this  equa- 
tion to  the  data  of  Westley  and  Woolley  I’esults  in  the  non- 
dimensional  curves  of  Figure  2fl. 

A modification  to  Equation  (22)  can  be  made  to  peimit 
its  application  to  airplane  cruise  conditions,  under  which  con- 
ditions discrete-frequency  shock  noise  has  been  found  to  be  a 
problem  on  several  airplanes  [42,  43].  'This  modification 
involves  correcting  for  the  characteristic  impedance 
of  the  ambient  air;  with  it  Equation  (22)  becomes 


where  Cgj^  are  sea  level  values  of  ambient  density  and  sound- 
speed,  and  where  and  are  corresponding  values  at  the  cruis- 
ing altitude. 

Although  the  preceding  discussion  has  been  directed  primarily 
at  convergent  nozzles,  the  method  can  be  applied  also  to  convergent- 
divergent  nozzles  which  are  operating  away  from  the  design  pressure 
ratio. 
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28  NON-DIMENSIONAL  SOUND  PRESSURE  LEVELS 
FOR  DISCRETE  FREQUENCY  SHOCK  NOISE 
(Based  on  [40]). 


As  for  broadband  shock  noise  (Section  2.^.1),  it  can  be 
assumed  here  that  no  discrete-frequency  shock  noise  occurs  for 
a convergent-divergent  nozzle  operating  at  pressure  ratios  be- 
tween 0.9  and  1.1  times  the  design  pressure  ratio.  (This  is  a 
gross  simplification  of  the  problem,  requiring  considerable  ex- 
perimental study.) 

It  is  instructive  to  compare  results  calculated  by  means 
of  the  proposed  procedure  for  prediction  of  the  upper  bounds  of 
the  sound  pressure  levels  with  the  limited  amount  of  published 
data.  Figure  IV-A-3-16  of  Knott  et  al,  [38]  shows  that  a 
model  Jet  with  total  temperature  of  2200°R  and  nozzle  pressure 
ratio  of  3.066  produced  at  x « 0 y = 3*630,  a sound  pressure 
le/el  of  12^  dB  at  the  I68O  Hz  shock  noise  fundamental  frequency. 
The  proposed  prediction  procedure  yields  a corresponding  upper- 
bound  sound  pressure  level  of  1^^  dB  at  a fundamental  frequency 
of  l^i|0  Hz.  The  flight  test  data  reported  by  Hay  [43]  for 
an  engine  with  a nozzle  diameter  of  22  inches  and  pressure  ratio 
of  3-3»  operating  at  an  altitude  of  35 » 000  ft,  show  a maximum 
pressure  level  of  140  dB,  measured  on  the  horizontal  stabilizer. 
The  proposed  prediction  procedure  yields  a corresponding  upper 
bound  of  148  dB  i^r  the  second  harmonic,  if  a 6 dB  increase  is 
Included  to  account  for  pressure  doubling  due  to  surface  reflec- 
tion. In  this  case  the  frequency  of  the  second  harmonic  is  in- 
fluenced by  the  airplane  speed. 

These  comparisons  appear  to  justify  the  use  of  the  upper 
bound  estimation  procedure.  However,  structures  designed  to 
withstand  the  estimated  acoustic  pressures  may  be  excessively 
heavy,  so  that  development  of  a more  precise  method  is  desirable. 


V 2.6  Forward  Motion  Effects 

Although  there  exists  much  empirical  evidence  that  forward 
flight  reduces  the  acoustic  power  generated  by  a Jet  flow,  the 
mechanism  Is  not  fully  understood.  Consequently,  no  widely 
accepted  method  Is  available  for  predicting  the  effect  of  for- 
ward motion  even  on  the  acoustic  farfleld,  which  Is  the  only 
region  that  has  been  Investigated.  Nevertheless,  the  available 
farfleld  results  here  will  be  used  as  a basis  for  the  construc- 
tion of  an  Interim  nearfleld  prediction  procedure,  and  the 
available  data  (which  correspond  mainly  to  turbojet  exhausts 
which  are  shock- free)  will  be  extrapolated  to  include  shock 
noise  and  coaxial  jet  effects. 

2,6.1  Shock-free  jet  noise 

Forward  motion  affects  the  jet  noise  field  in  terms  of  its 
acoustic  power  generation,  a Doppler  frequency  shift,  and  the 
propagation  distance  and  angle.  The  change  In  acoustic  power 
generation  has  been  taken  into  account  in  several  investigations 
by  Introduction  of  the  relative  velocity  (V.-V  ),  where  V. 

J ft  J 

denotes  the  jet  velocity  relative  to  the  noszle  and  the 

airplane  (or  forward)  v velocity . In  early  studies,  Vj  was 

simply  replaced  by  (V.-V  ).  Subsequent  analyses  reduced  the 

J ft 

ro.le  played  by  the  relative  velocity:  Pfowes  Williams  [45] 

7 /8 

suggested  that  be  replaced  by  Vj  (1-V^/Vj)"  , Kobrynski 
[46]  replaced  by  Vj (1-V^/Vj )^/^  and,  on  the  basis  of 

recent  experimental  data  [53],  Stone  [20]  proposed  that  the 

term  V,  (1-V„A,  be  used  in  place  of  V.. 

j a j j 

Analytical  studies  by  Plnkel  [4?]  showed  that  the  influence 
of  forward  flight  is  more  complicated  than  Is  suggested  by  the 
above  relative  velocity  terras.  However,  If  Plnkel’ s results 
are  simplified  co  the  form  Yj (1-V^/Vj )”,  the  exponent  n is  found 
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to  be  between  1/2  and  3/^*  Cocking  and  Bryce  [48]  also  ob- 
tained exponent  values  In  the  same  range.  Examination  of  the 
aforementioned  experimental  and  analytical  results  indicates 
that  an  exponent  of  3/4  represents  a reasonable  average  value; 
this  value  therefore  la  proposed  for  the  present  prediction  model. 

Although  one  might  at  first  glance  expect  no  Doppler  fre- 
quency shift,  because  the  airframe  and  engine  move  together  as 
a rigid  body,  Pranken  et  al,  [49]  found  no  experimental  evidence 
to  support  this  expectation  with  respect  to  jet  noise.  Pranken 
et  al.  adopted  Powell's  hypothesis  [50]  that  the  airplane  moves 
away  from  the  jet  noise  sources,  with  the  result  that  Doppler 
frequency  shifts  do  occur,  and  also  that  there  result  changes  in 
the  distance  travelled  by  a sound  wave  from  a source  to  the  air- 
frame, and  in  the  angle  of  propagation  associated  with  the  given 
receiver  location.  The  method  of  Pranken  et  al. , which  covers 
these  factors  has  been  adopted  by  Rudder  and  Plumblee  [51]  for 
nearfleld  noise  and  is  reproduced  here,  in  the  absence  of  any 
better  available  method.  The  notation  used  here,  however,  is  con- 
sistent with  that  of  other  parts  of  this  section. 

Pranken  et  al.  [49]  assume  that  forward  motion  of  the  air- 
frame moves  the  receiver  from  its  actual  position  r relative  to 

s 

the  source  to  an  apparent  position  r*  on  a line  parallel  to  the 

s 

direction  of  motion.  (See  Figure  29.)  Correspondingly,  the 
angle  between  the  Jet  axis  and  the  source/receiver  line  is  changed 

from  6!  to  4* . Pranken  et  al  assumed  that  all  sources  are  close 

8 

to  the  Jet  exit  plane,  but  this  assumption  is  neither  necessary 
nor  desirable;  Figure  29  shows  a typical  noise  source  located 
at  distance  downstream  of  the  exit.  (Primed  symbols  are  used 
here  to  Indicate  that  the  noise  sources  are  assumed  to  be  located 
on  the  Jet  centerline,  whereas  the  Jet  no^se  prediction  procedure 
of  Plumblee  et  at,  [6]  assumes  a source  distribution  along  the 
Jet  boundary.) 
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In  order  to  determine  the  relationship  between  b*  and  <j>*, 

S 8 

one  may  note  that  the  sound  wave  travels  a distance  r*  at  velo- 
city In  the  same  time  In  which  the  receiver  moves  through  a 
distance  Ax  at  velocity  V^,  so  that 


« 


c 

o 


(2^0 


Prom  Figure  29, 


(r*)2  = y2  + (x-x^  - Ax)^ 

x-x 

since,  by  definition,  cot  0'  = 

s y 

x-x  -Ax 
cot  . -f— 

and 

a a o' 


(25) 


one  finds  that 
1 


cot 


(1-M„  ) 


Lot  e; 


- M_  V l-M.*+cot*e: 


a 


a 


(26) 


This  relation  between  0'  and  Is  plotted  in  Figure  30, 

from  which  one  may  note  that  for  supersonic  flight  speeds  no  Jet 

noise  propagates  upstream  of  the  source,  as  one  would  expect. 

If  one  assumes  a stationary  source  and  a moving  receiver, 
one  finds  that  there  occurs  a Doppler  shift,  such  that  the  ob- 
served frequency  f*  on  the  airframe  Is  related  to  the  source  fre- 
quency f as 

f*  • (1+M^cos  ♦')  f 
& 0 
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HG.  30.  EFFECTS  OF  FORWARD  MOTION  ON  ANGLE 
BETWEEN  JET  AXIS  AND  LINE  FROM 
SOURCE  TO  RECEIVER  [49]. 


This  shifted  frequency  needs  to  be  taken  Into  account  when  one 
calculates  the  Strouhal  number  associated  with  a given  fre- 
quency; for  example,  if  f*  is  the  observed  frequency  on  the 
airframe,  the  Strouhal  number  should  be  calculated  correspond- 
ing to  the  frequency  f » f*/(l+M_cos 

cl  S 

For  illustrative  purposes,  Pranken  et  at  . [49]  applied 
their  method  (which  assumes  that  the  c>ound  sources  are  located 
near  the  Jet  exhaust  nozzle)  to  the  nearfield  contours  of  a 
turbojet  engine  for  a flight  Mach  number  M « 0.8.  Figure  31 
shows  the  corresponding  nearfleld  sound  pressure  level  contours 
for  static  operation  and  forward  motion. 

It  must  be  noted  that  the  hypothesis  that  the  receiver 
moves  away  from  the  source  was  originated  by  Powell  [50], 
passed  through  the  literature  by  Pranken  at  at.  [49],  and  ac- 
cepted by  Rudder  and  Plumblee  [51]  without  any  direct  support- 
ing evidence.  The  hypothesis  was  based  on  the  observation  that 
the  significant  Increase  in  sound  pressure  level  that  is  expected 
upstream  of  rigidly  connected  sources  moving  near  Mach  1 was  not 
observed  on  rocket  structures.  No  experimental  da^;a  have  been 
published  that  show  that  there  exists  a Doppler  shift,  and  no 
data  are  available  that  show  that  there  exists  no  Doppler  shift. 
Under  these  circumstances,  the  analysis  incorporating  a Doppler 
shift  is  presented  here  with  strong  reservations. 

2.6.2  Coaxial  Jets 

The  effects  of  flight  velocity  on  the  noise  produced  by  co- 
axial jets  have  been  studied  solely  with  respect  to  the  acoustic 
farfield.  The  corresponding  results  are  reviewed  here  in  terms 
of  acoustic  power  generated,  since  the  propagation  effects  here 
will  be  essentially  the  same  as  those  for  a single  Jet  (Section 

2.5.1). 


FIG.  31.  EFFECT  OF  FORWARD  MOTION  ON  NEAR-FIELD  JET  NOISE [49]. 


The  results  of  vai^lous  published  studies  show  some  diver- 
gence. Reed  [52]  end  von  Gl&hn  *»t  al.  [53]  consider  a dual- 
flow Jet  as  a single  mixed  flow>  and  relate  relative  velocity 
to  the  primary  Jet  velocity.  In  contrast,  Cocking  and 
Bryce  [48]  and  Dunn  and  Peart  [35]  consider  the  primary 
and  secondary  flows  separately.  Most  of  the  available  measure- 
ments were  obtained  from  ground  tests,  using  a wind  tunnel  or 
tertiary  Jet  to  simulate  the  forward  velocity  effect.  (Such 
tests  have  the  advantage  that  they  exclude  the  Doppler  shift 
effect,  but  they  have  the  disadvantage  that  the  microphone  either 
Is  Immersed  In  the  moving  alrstream  or  outside  the  flow,  In  which 
case  there  may  be  refraction  effects.)  To  minimize  the  experi- 
mental uncertainties,  the  various  reported  effects  of  forward 
motion  are  compared  In  Table  5 for  an  angle  of  90°  to  the 
Jet  axis.  This  table  lists  the  relative  velocity  terms  for  the 
primary,  secondary  or  combined  flows,  which  terms  replace  the 
appropriate  Jet  velocity  term  in  the  prediction  of  the  acov.stic 
farfield. 

The  relative  velocity  exponents  indicated  in  Table  5 
show  considerable  variation,  as  was  the  case  for  single  Jet  flows. 
The  results  of  Dunn  and  Peart  show  unity  exponent,  but  It  Is 
probable  that  little  data  on  relative  velocity  effects  was  avail- 
able at  the  time  of  the  study:  all  later  studies  show  exponents 
of  lens  than  unity.  However,  the  work  of  Dunn  and  Peart  does 
take  into  account  the  area  ratio  A^/A^,  a factor  that  should  be 
Included  In  any  general  prediction  method. 

Based  on  the  small  amount  of  available  Information,  it  ap- 
pears reasonable  to  assume  that  forward  motion  of  the  airplane 
will  not  Influence  the  noise  contribution  from  the  primary  * jw, 
and  that  the  relative  velocity  effect  for  the  secondary  flow 
may  be  taken  to  be  similar  to  that  for  single  Jets.  Thus,  one 
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TABLE  5. 


RELATIVE  VELOCITY  TERMS  FOR  ACOUSTIC  FAR-FIELO 
OR  COAXIAL  JETS. 


REFERENCE 

TERM  REPLACING  JET  VELOCITY,  TO  YIELD 
CORRECT  SOUND  PRESSURE  VELELS  AT  90*> 
FROM  JET  AXIS 

Reed  [52] 

Vjj  j ; n - 0.63  to  0.79 

Von  Glahn  et  al.  [531 

1 ^12 

Cocking  and  Bryce  [48] 

Primary:  / 

/ \ \0.63 

secondary:  Vjj  |l  - y-  j 

Dunn  and  Peart  [35] 

Primary: 

Secondary* : ^ ^ j 

« 


.11  * ^2^.12 
Aj  Aj 


U 


TIm  •xprttcalon  given  here  eceumee  the  angle  between  the  gross  thrust 
veotor  and  the  direction  of  notion  to  be  small. 
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may  replace  the  effective  flow  velocity  in  the  equation  for 

nearfleld  sound  presure  levels  by  the  term  V.(l-V 

j a j 

where 


P A V, 


+ p A V 


r ” V 4 

V.  i-  1. Jj 

J p A v”  + p A 


I. 


i t 


u 


Justification  of  this  choice  of  parameters  must  wait  until 
corresponding  nearfleld  experimental  data  are  available. 


2.6.3  Shock  noise 

There  is  little  evidence  available  concerning  the  effect  of 
forward  flight  on  shock  noise.  Thus,  it  appears  appropriate  for 
the  present  to  assume  that  forward  motion  has  no  influence  on 
the  acoustic  power  generated  by  both  discrete-frequency  and 
broadband  shock  noise.  However,  one  may  expect  a spatial  redis- 
tribution of  acoustic  energy,  as  well  as  changes  in  the  frequency. 

In  the  case  of  broadband  shock  noise,  it  is  reasonable  to 
make  the  same  assumptions  regarding  Doppler  shift  and  acoustic 
propagation  as  those  made  for  shock-free  jet  noise.  These  are 
discussed  in  Sections  2.5.1  and  2.5.2,  and  will  not  be  repeated 

here. 


The  influence  of  forward  motion  on  nearfield  discrete- 
frequency  shock  noise  has  been  investigated  by  Hay  [i<3]  and 
Hewlett  [^<2]  under  actual  flight  conditions.  However,  these 
investigations  refer  only  to  frequency,  and  do  not  show  the 
round  pressure  level.  Both  studies  are  based  on  the  result  of 

Powell  [19],  which  gives  the  fundamental  frequency  in  still  air 

as 

• _ °o  1 , 

^1  • Kp  • (28) 

with  « 3. 
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Both  studies  develop  seml-emplrlcal  r-iatlons  for  the  shock 
noise  frequency  by  adjusting  the  time  taken  for  the  sound  wave 
to  travel  in  the  upstream  portion  of  the  feedback  loop  (i.e., 
in  the  external  flow). 


The  equation  developed  by  Hay  for  the  fundamental  fre- 
quency f^  Is 


[M  + K,  (M,  - M )]  (1  - M ) (1  + K,H,  ) 


(29) 


where 


R ■ Jet  pressure  ratio 

R ■ critical  Jet  pressure  ratio  « 1.89 
c 

M.  ■ Jet  Mach  number  relative  to  ambient 
^ speed  of  sound 

■ Airplane  Mach  number 

D ■ Jet  diameter 
Cq  ■ ambient  speed  of  sound 
Kg  - 0.625 
» 1.1 


There  Is  some  divergence  here  between  the  « 3 of  Powell  and 
■ 1.1  of  Hay,  but  the  results  of  several  investigations  at 
zero  forward  speed  also  show  considerable  variation  In  the 
value  of  K^,  as  was  discussed  In  Section  1.4.  Thus,  for  a 
general  prediction  procedure.  It  appears  appropriate  to  use 
Equation  (29)  with  “ 2.0,  so  that  the  prediction  methods 
for  static  and  moving  aircraft  would  be  consistent.  The  corre 
latlon  between  Equation  (2?)  and  related  experimental  data 
on  BACl-11  and  VClp  alrple  .es  [54]  Is  shown  In  Figure  32. 
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(l-M.) 

[Mo  ♦ 0.625  (Mj-AA„) 

] (U0.625Mj) 

(R-R^)''‘  [k 0.625 (M|-Mj] 

0.625  AAj 

FIG.  32.  SHOCK  NOISE  FUNDAMENTAL  FREQUENCY  AS  A 
FUNCTION  OF  JET  AND  AIRPLANE  MACH 
NUMBER  [54]. 
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Th6  fx*6(|U6ncy  fg  of  t»hG  sGcond  harmonic  Is  given  by 

^2  " (30) 
as  Is  the  case  for  zero  forward  motion. 

The  spatial  distribution  of  the  sound  pressures  associated 
with  discrete-frequency  shock  noise  will  be  affected  by  forward 
motion,  as  for  broadband  shock  noise.  However,  no  Doppler 
shift  may  be  expected,  because  the  mechanism  of  the  discrete- 
frequency  noise  is  fixed  in  relation  to  the  Jet  nozzle. 


2,7  Thrust  Reversers 

Thrust  reverser  noise  has  been  investigated  only  recently  and 
such  little  information  as  is  available  refers  to  far-field  sound. 
M »t  of  the  available  data  is  concerned  with  target-type  thrust 
reversers.  Stone  [20]  has  summarized  the  data  into  prediction 
procedures  for  far-fleld  noise,  and  his  summary  forms  the  basis 
of  the  present  discussion. 

2.7.1  Target-type  reversers  ^ 

Target-type  thrust  reversers  (Figure  33)  have  been  inves- 
tigated by  several  authors  [55*  56,  57]  and  their  results 
have  been  anaA-yzed  by  Stone  [20]  for  the  purpose  of  developing 
generalized  relationships  for  overall  sound  pressure  levels  and 
one-third  octave  band  spectra.  This  analysis  indicates  that  the 
radiated  sound  is  dominated  by  noise  due  to  interaction  between 
the  Jet  and  the  airplane  structure;  there  is  no  evidence  that  the 
exhaust  Jet  noise  itself  plays  a significant  role. 

Stone  proposes  the  following  empirical  relation  for  the  far- 
field  overall  sound  pressure  level  at  an  angle  of  90°  to  the  in- 
let axis: 


OASPLgo  - ■.  10  log  (Aj/r*)  f 10  log  |pjP,cVp^3^ 


- 10(cos*»j/)  logCO^/Dg]  + 55  log 

- 10  log  [1  + O.OKV./c  )*-®]  (33 

J 0 

where  t)^  is  the  angle  between  the  source  to  observer  line  and  the 
y axis 

( 1^9  for  semi-cylindrical  reversers 


and  K. 


15^  for  V-gutter  reversers 
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FIG.  33.  DiAGRAMHATIC  REPRESENTATION  OF  TARGET-TYPE  THRUST  REVERSERS. 


stone  configured  the  velocity  terms  arbitrarily  to  give  a 
relationship  at  high  velocities",  to  match  his  formulation  for  Jet 

noise  at  high  supersonic  speeds.  However  for  V./c  the  inten- 

5 R J ° 6 

slty  varies  approximately  as  * which  is  similar  to  the  Vj 

relationship  for  dipole  noise.  Dipole-like  noise  sources  would 

be  expected  to  occur  from  sources  associated  with  flow  Impact  on 

relatively  small-scale  structures. 

As  an  Interim  procedure,  one  may  take  the  directivity  in 
the  6-plane  to  be  given  by 

OASPLg  - OASPLgj  + 20  log  (32) 

This  differs  from  the  directivity  corresponding  to  a single 
dipole,  but  this  Influence  is  not  surprising,  because  in  a thrust 
reverser  one  encounters  dipole  sources  arrayed  along  several 
axes. 


It  is  likely  that  in  the  nearfield  the  sound  does  not  spread 
according  to  the  farfield  inverse-square  law.  For  jet  noise, 
Plumbblee  et  al  [6]  and  Chen  [1]  used  the  formulation  of  Pink 
[57]  for  the  nearfleld  of  quadrupole  noise  sources.  Franz's 
work  also  Includes  a formulation  for  the  nearfielci  of  a dipole 
source;  this  result  will  be  adopted  here  for  thrust  reversers. 

In  order  to  represent  this  distance-dependence,  one  may  add  to 
the  right-hand  side  of  equation  (3D  the  term 

AOASPL  • 10  log  1^1  + (cQ/rw^)^j  (33) 

where  is  a characteristic  radian  frequency,  discussed  below. 


For  Jet  nolee,  Chen  [1]  Introduced  a characteristic 
frequency  that  changes  with  axial  position  In  the  Jet.  In  con- 
trast, Plumblee  et  at  [6]  obtained  an  empirical  relationship 
that  Is  Independent  of  frequency. 


It  is  proposed  that  for  target-type  thrust  reversers  the 
characteristic  frequency  be  taken  as  the  peak  frequency  in  the 
noise  spectrum,  which  frequency  is  given  by 

2ir 

% “ U (TVT 

6 n G J a 

1/2 

where  D » equivalent  circular  nozzle  diameter  = (4A/it) 

* hydraulic  diameter  = (llA/perimeter) 

6"  *=  effective  angle  = 0(Vj/c^)®'^ 

Tj  **  Jet  total  temperature 

It  is  assumed  here  that  the  noise  sources  associated  with 
thrust  reverser  operation  move  with  the  airframe.  Consequently, 
no  Doppler  shift  corrections  have  been  applied. 


A directivity  correction  term  to  account  for  the  effect  of 
forward  motion  of  the  aircraft  should  be  applied  to  the  overall 
sound  pressure  level.  This  term  is 


AOASPL  » 10  log 


(l  . 


cos8 


/ 


(35) 


However,  since  thrust  reverser  operation  occurs  only  at  low 
airplane  Mach  numbers,  the  correction  usually  is  likely  to  be 
negligible. 


The  frequency  distribution  of  target  thrust  reverser  noise 
has  been  represented  by  Stone  [20]  in  terms  of  a single  reduced 
one-third  octave  band  spectrum.  The  collapse  cf  the  data  is  shown 


95 


t 


by  Figure  3^»  together  with  the  curve  proposed  by  Stone  for  far- 

field  noise  prediction.  This  proposed  curve  has  its  maximum  at  a 
% 

Strouhal  number  of  unity,  and  this  value  was  used  in  deriving 
equation  (39).  The  Strouhal  number  S is  defined  by 


— 


(36) 


J 


where  it  has  been  assumed  that  no  Doppler  shift  occurs  for  receiv- 
ers located  on  the  airplane. 

2.7.2  Cascade-type  reversers 

Information  on  the  noise  fields  of  cascade-type  thrust  re- 
versers (Figure  35)  is  extremely  sparse,  so  that  any  prediction 
method  must  be  regarded  as  very  approximate.  The  small  amount  of 
data  available  [38]  suggests  that  cascade-type  reversers  are 
somewhat  quieter  than  those  of  the  target-type;  consequently, 
jet  noise  may  make  a significant  contribution  to  the  farfield  (and 
presumably  the  nearfield)  noise.  The  method  proposed  by  Stone 
accounts  for  the  Jet  noise  contribution  in  addition  to  that  a- 
rlsing  from  interaction  of  the  jet  with  the  reverser  structure. 

In  a manner  similar  to  that  used  for  target-type  reversers.  Stone 
derives  an  empirical  equation  for  the  farfield  noise  at  an  angle 
of  90®  to  the  inlet  axis  of  a cascade-type  reverser: 


OASPL 


90.cr  ■ >'cr  * " 1“  ) 


♦ 50  log  (v./c  ) - 10  log  [1  + 0.0l(V./c  )*] 
J • J • 


(37) 
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FIG.  34.  NOfl-DIMENSIOHAL  ONE-THIRD  OCTAVE  BAND  SPECTRUM 
FOR  TARGET-TYPE  THRUST  REVERSERS  [60] 


Here  the  empirical  constant  Is  given  by 


■ "-sr.2  * <V*t> 


(38) 


with 

0 If  there  Is  an  Internal  flow  deflector  to  guide 

K , ■ the  flow  Into  the  cascades 
cr 

5 If  there  Is  no  deflector 


!0  If  alrfoll-shaped  vanes  are  used 
6 if  con&tant-thlckness  vanes  are  used 

Ag/A^  • ratio  of  cascade  exit  area  to  tailpipe  area 
(generally  greater  than  unity) 


Equation  (36)  can  be  adapted  to  nearfield  conditions  by  the 
same  reasoning  as  for  target-type  reversers:  an  additional  term, 

given  by  equation  (33)  is  added  to  the  right-hand  side  of 
equation  (37)  and  the  characteristic  frequency  is  here 
given  by 

2.8ir  Vj 

“o  * D^'  (D^/D^)®-'*  ♦ ■=<>= 

Where  and  are  defined  In  Figure  35,  and  3”  has  been 
defined  earlier.  As  before, aperiraental  Justification  of  this 
nearfleld  modification  is  recuired. 

Following  the  proposed  procedure  of  Stone  [20],  the  over- 
all sound  pressure  level  at  90®,  as  estimated  by  equations  (37) 

and  (39)  is  combined  with  the  predicted  Jet  noise  contribution 
to  give  the  total  overall  sound  pressure  level: 

r (OASPLqo  „-)/10  (OASPLon  < )/lol 

OASPL  • 10  log  jlO  90, cr  10  J (AO) 
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’’•asurements  of  cascade-rype  thru'^t  reverser  noise  show 

directivity  patterns  which  are  complicated  In  shape,  althougn  the 

variations  are  relatively  small  In  magnitude.  Stone  proposes  a 

simplified  relationship  for  the  directivity  of  the  farfleld  noise: 
% 

OASPLq  - OASPL^q  + 20  log  (1  + H sin  26) 

It  is  proposed  that  this  also  be  used  for  the  r^earfleld. 

The  effect  of  forward  motion  Is  more  difficult  to  Include  for 
cascade-type  type  than  for  target-type  reversers.  For  the  Jet- 
structure  Interaction  noise,  the  same  approach  can  be  followed  as 
In  the  case  of  target  reversers,  with  the  correction  term  being 
given  by  equation  (35)*  For  the  Jet  exhaust  noise,  the  situa- 
tion is  different  from  that  of  a typical  Jet  engine,  because  with 
reversers  the  exhaust  axis  is  not  parallel  to  the  direction  of  the 
airplane  motion,  so  that  the  relative  velocity  factors  discussed 
in  Section  1.5  do  not  apply  here.  Extensive  analysis  will  be 
required  to  establish  a detailed  procedure.  For  the  present,  it 
Is  proposed  that  the  correction  term  given  by  equation  (35)  be 
used  for  the  Jet  noise  also.  Thus,  the  correction  can  be  applied 
directly  to  equation  (^<1). 

The  one-third  octave  band  spectrum  levels  measured  for  cas- 
cade-type thrust  reverser  have  been  reduce i by  Stone  [35]  as 
shown  in  Figure  36.  This  figure  also  shows  a curve  proposed 
by  Stone  for  a farfleld  prediction  procedure.  It  Is  suggested 
that  the  same  curve  be  used  for  nearfield  noise  prediction,  with 
the  overall  sound  pressure  level  given  by  equation  (ill),  the 
Strouhal  number  given  by  equation  (36),  and  assuming  that 
there  is  no  Doppler  shift. 
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2.8  Reflection  of  Acoustic  Pressures  At  The  Fuselage 

The  preceding  prediction  methods  pertain  to  acoustic  pres- 
sures that  occur  under  free-fiold  conditions,  where  there  are 
no  reflecting  surfaces.  In  a practical  environment,  however, 
there  occur  reflections  due  to  the  presence  of  the  ground  plane 
and  of  the  airframe  structure.  Only  the  effect  of  the  latter 
Is  discussed  here. 

It  is  well  known  that  for  the  simplest  case  of  normal  in- 
cidence of  an  acoustic  v;ave  on  a large  plane  surface,  reflection 
at  the  surface  produces  a 6 dB  increase  In  the  sound  pressure 
level  at  the  surface.  However,  for  grazing  incidence,  there  oc- 
curs no  Increase,  and  for  oblique  incidence  on  a curved  surface, 
the  Increase  lies  between  0 and  6 dB.  Increases  greater  than  6 
dB  can  occur  at  the  Intersection  of  two  planes,  such  as  the 
Junction  of  the  horizontal  stabilizer  and  the  rear  fuselage 
structure. 

In  considering  the  influence  of  only  the  cylindrical  fuse- 
lage, one  can  analyze  the  two  limiting  cases:  (1)  that  whex'e 

an  engine  is  mounted  near  the  fuselage,  and  (2)  that  where  the 
engine  is  at  a considerable  distance  from  the  fuselage.  These 
two  cases  were  analyzed  by  Cockhurn  and  Jolly  [ll]  and  the 
results  were  applied  by  Rudder  and  Plumblee  [1^4]  to  Jet  en- 
gine noise. 

2.8.1  Engines  Close  To  The  Fuselage 

For  an  engine  mounted  within  15  nozzle  diameters  of  the 
fuselage,  as  illustrated  in  Figure  37  for  a rear-mounted  en- 
gine, an  analysis  based  on  scattering  is  not  appropriate. 

Rudder  and  Plumbloe  calculate  the  pressure  p at  a point  on  the 


Ff«o  Fi«l«i  Sound  Prosturo) 

FIG.  37.  GEOMETRY  AND  PRESSURE  DISTRIBUTION  FOR 
ENGINES  MOUNTED  CLOSE  TO  THE  FUSELAGE 
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fuselage  surface  from 


3 + cos  2a 
2 


(42) 


where  a is  the  angle  between  the  wavefront  and  the  tangent  to 
the  fuselage  surface  (-TT/2<a<7r/2)  and  is  the  free-field  pres- 
sure. A typical  pressure  loading  on  the  fuselage  is  shown  at 
the  bottom  of  Figure  37  • 

Equation  (^2),  when  rewritten  in  terms  of  the  sound 
pressure  level  (SPL)^,  becomes 


(SPL)„  = (SPL)^  + 20  log 
s o 


[ 


3-tcos2a 
2 


(43) 


2.8.2  Engines  Far  From  The  Fuselage 


For  an  engine  mounted  more  than  15  Jet  nozzle  diameters 
from  the  fuselage,  one  may  use  the  acoustical  scattering  analy- 
sis of  Potter  [593.  A typical  geometric  configuration  for 
this  case  is  shown  in  Figure  38,  where  a source  location  is 
identified  in  the  Jet  flow.  The  analysis  Involves  determination 
of  the  sound  field  for  a single  source  S at  a point  Q on  the  fuse- 
lage surface,  for  single  frequency  w,  and  later  accounting  for 
all  sources  and  frequencies. 


The  analysis  assumes  that  the  acoustic  wave  arriving  at  the 
fuselage  is  a plane  wave  of  strength  p^  and  with  the  wave  direc- 
tion such  that  the  normal  to  the  wave  front  makes  an  angle  3 
with  the  normal  to  the  axis  of  the  fuselage  (see  Figure  39). 

The  point  on  the  surface  of  the  fuselage  is  specified  by  coordi- 
nate (a,  X,  (|i),  where  a Is  the  radius  of  the  fuselage,  x the 
dimension  along  the  axis  and  the  angle  around  the  fuselage 
measured  from  the  plane  containing  the  source  and  the  fuselage 
center-line.  («|>  is  zero  when  the  point  Is  on  the  far  side  of 


FIG.  38.  GEOMETRY  OF  FUSELAGE  AND  ENGINE  MOUNTED  FAR  FROM 
THE  FUSELAGE. 


the  fuselage^  completely  hidden  from  the  source.) 

The  surface  pressure  has  been  shown  [59]  to  be  given 
by  the  real  part  of 


P 


ka  CCS0  ^ 
m-0 


where 


SP.g. Jit  exp[l(-Y^  + •.rm/2)]  expClkx  slnP^]  exp[-2irlft] 


C ' ■ 2 f J*  (ka  COS0)  + N*  (ka  cos8)  ] ^ 

oil  I J 
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Y t ■ 
'O 


ivi 

^ I 

tan 


(ka  cosP)  - Nju^j^(ka  cos0)  | ^ 

[-J^(ka  cos3) 

Klj(ka  cosdl 


y’ 

'm 


tan 


Jj^^l(ka  co9$)  - cos8) 

Njjj_j^(ka  cosP)  - Vi  (ka  co8$) 


(45) 


and  where  J and  N are  Bessel  functions  of  the  first  and  second 
kind,  respectively. 

It  should  be  noted  that  this  result  Is  applicable  only  to 
wlng^mounted  engines » where  the  sound  field  for  a given  source 
can  be  approximated  as  a plane  wave  at  the  point  where  It  strikes 
the  fuselage;  the  effects  of  apherloel  radiation  have  not  been 
included.  INirtherdore,  the  surface  of  the  fuselage  was  assumed 
to  he,  perfectly  rigid,  which  is  not  the  once  in  practice.  The 

' ■ . 
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complete  solution  would  be  extremely  complicated  and  the  present 
method  is  considered  to  be  sufficiently  accurate  for  practical 
situations. 

Rudder  and  Plumblee  present  equation  (l|4)  In  a more  con- 
venient form  for  computatlonj  they  express  the  pressure  at  a 
point  on  the  fuselage  surface  due  to  an  obliquely  Incident  plane 
wave  at  frequency  f as 

P 
A* 


where 


* Ap_cos(2irft  + #)  (46) 

o 

- C"  + n* 


t ! 


ii 


is  streni^^ 


♦ ■ tan“*  (n/4) 
C » Dq  (EG-PH) 
n - (PG+EH) 

D - 5- 

o irka  cos3 


0 « cos  (kx  sin?) 

H -sin  (kx  sin?) 

of  the  free-fleld  sound  pressure. 


2.9  Pressure  Correlations 


Although  pressure  correlation  coefficients  in  the  near- 
field  of  model  and  fullscale  Jets  have  been  measured  by 
several  investigators  [16,  30,  60-65],  very  little  effort  has 
been  directed  to  providing  a corresponding  generalized  analyt- 
ical model,  and  no  attempt  appears  to  have  been  made  at  study- 
ing the  effects  of  parameters,  such  as  nozzle  pressure  ratio 
and  Jet  temperature.  The  published  information  on  nearfield 
pressure  correlations  is  reviewed  here,  so  that  a provision- 
analytical  model  can  be  proposed  for  use  in  a prediction 
procedure. 

2.9.1  Subregions 

Clarkson  C65,  66]  postulated  that,  from  the  point 
of  view  of  pressure  correlations,  the  region  around  a Jet  can 
be  divided  into  two  subregions: 

In  Subregion  1 (near  the  Jet  boundary)  the  pressure  fluctu'a- 
tlons  are  strongly  Influenced  by  the  convection  of  the  aero- 
dynamic turbulent  field.  This  subregion  extends  up  to  about 
2 nozzle  diameters  out  from  the  Jet  boundary,  and  downstream 
to  a plane  about  12  to  15  diameters  from  the  nozzle  exit  plane. 

In  Subregion  2 (nearfield)  the  spatial  distribution  of  the 
noise  sources  is  on  a scale  comparable  to  the  dimensions  considered 
in  correlation  measurements.  This  subregion  extends  to  about 
10  diameters  out  from  the  Jet  boundary,  20  diameters  downstream 
and  10  diameters  upstream  of  the  nozzle  exit  plane.  The  posi- 
tion of  the  outer  surface  of  this  subregion  Is  not  clearly  de- 
fined and  is  probably  frequency -dependent . 


Clarkson  does  not  Identify  the  criteria  which  he  used 
to  define  these  two  subregions. 


Eldred  et  al  [25]  observed  a conflict  between  differ- 
ent sets  of  published  experimental  data:  some  indicated  tnat 

the  close-in  correlations  were  dominated  by  the  aerodynamic 
field,  while  other  results  showed  domination  by  acoustic  radi- 
ation. Eldred  et  al  argued  that  the  radiated  sound  field 

' g 

varies  on  the  average  as  V.  , whereas  the  aerodynamic  pressure 

ii  J 


field  increases  as  V 


J 


so  that  the  field  of  influence  of  the 


aerodynamic  pressures  would  decrease  as  V.  increases. 

J 


2.9.2  Freefield  Data 

Figure  ^0  summarizes  the  locations  at  which  various 
investigators  have  carried  out  correlation  measurements.  (In 
most  cai  .0  the  correlation  coefficients  were  measured  along 
axes  parallel  to  the  Jet  centerline.)  At  locations  identified 
by  a angular  symbol  only  broadband  correlation  measurements 
were  me.  ej  at  those  indicated  by  a circular  symbol,  narrowband 
(one-thj.*'d  or  one-half  octave  band)  measurements  were  made. 

Open  symbols  indicate  that  the  convection  velocity  was  super- 
sonic for  11  frequencies  measured;  fllled-in  symbols  indicate 
that  all  \ch  velocities  were  subsonic;  half  fllled-ln  symbols 
correspond  to  the  case  where  the  convection  velocities  were  sub- 
sonic at  some  frequencies  and  supersonic  at  others. 


Figure  ^0  also  indicates  the  extent  of  the  two  subre- 
gions postulated  by  Clarkson  [65].  It  is  seen  that  it  is 
difficult  to  reconcile  the  subregions  with  the  data  on  the 
basis  ot  convection  velocity. 
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Broadband  Corretatione 

Two  Investigations  have  resulted  In  analytical  models 
for  the  broadband  pressure  correlation  coefficient.  Unfor- 
tunately, In  both  cases  the  model  referred  to  only  one  loca- 
tion In  the  acoustic  field,  at  which  location  the  Investiga- 
tors also  made  response  measurements  on  a model  panel. 
Maestrello  et  at  [61l]  found  that  the  local  pressure  field 
was  dominated  by  aerodynamic  effects  and  proposed  an  analyti- 
cal model  for  the  pressure  correlation  coefficient  of  the 
form 

P(5  ,t)  - e“®>  Uos  o).(T-e, /U^) 

IS  u * u 

(n) 

where  U denotes  the  convect Ion  velocity  of  the  pressure  field, 
c 

T the  time  delay  and  and  the  separation  distances  in  the 
axial  and  lateral  directions,  respectively.*  The  coefficients 
in  the  exponents  of  this  equation  were  found  to  have  the  values 

a ■ 300  sec 
1 

W 0.H9  ff“^ 

o - 1.M2  ft’^ 

s 

and  the  convection  velocity  was  found  to  be  O.65  times  the  jot 
velocity  Vj . The  characteristic  frequency  Wq  corresponds  to 
the  peak  in  the  pressure  spectrum;  Maestrello  at  al  propose  for 


It  .is  Implicit  in  this  type  of  representation  that  the  pressure 
field  Is  homogeneous.  This  assumption  Is  not  correct  and  the 
representation  is  thus  only  an  approximate  form. 


the  frequency  at  any  axial  distance  x the  relation 


(t) 


0 


l.it 

[1  + 0.2x/D]^ 


where  D denotes  the  Jet  nozzle  diameter. 


(^<8) 


Mollo-Chrlstensen  [6l]  did  not  express  his  measure- 
ment results  for  the  aerodynamic  region  In  analytical  form, 
but  his  data  Indicate  an  average  convection  velocity  of  0.68 
Vj , which  is  close  to  the  above  indicated  value  obtained  by 
Maestrello  et  al. 


Olson  and  Llnderg  [633  made  measurements  in  the  acous- 
tic region  and  proposed  a model  of  the  form 


(‘I9> 


Where  c^  represents  the  ambient  speed  of  sound.  The  character- 
istic frequency  In  this  equation  Is  • 8000  radians/sec  (1273 
Hz)  and  the  coefficients  in  the  exponents  have  the  values 


a 

- 4500 

sec" 

a 

■ 0.60 

ft”^ 

s 

a 

- 0.44 

ft"^ 
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Navvowband  Correlations 

In  principle,  the  correlation  coefficients  given  by  equa- 
tions (47)  and  (49)  can  be  transformed  into  the  frequency 
domain  by  application  of  the  Fourier  transform.  However,  an 
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alternative  approach  that  is  widely  uaed  In  practice  consists 
of  fllterl/.g  the  pressure  signals  prior  to  correlation.  The 
correlation  coefficients  obtained  In  this  manner  are  functions 
of  frequency  and  also  depend  to  some  extent  on  bandwidth. 

Important  frequency-dependent  correlation  measurements 
were  performed  by  Howes  at  al  [30],  Cox  et  al  [16']  and 
Clarkson  [62,  65],  all  on  fullscale  Jet  engines.  The 
frequency  bandwldths  used  ranged  from  a one-third  octave  to 
slightly  over  a one-half  octave. 

Narrowband  correlation  coefficients  have  the  general  form 
of  a decaying  cosine.  The  experimental  results  refer  mainly  to 
the  spatial  correlation  coefficient  (for  t®0),  which  can  be 
represented  Uy  a function  of  the  type  cos  (wC/U^,),  where  w is  the 
center  frequency  of  the  band.  In  the  present  context,  the  two 
important  parameters  of  the  correlation  coefficient  are:  (1) 

the  first  zero  crossing,  which  provides  a measure  of  the  length 
scale  of  the  correlation  coefficient,  and  (2)  the  spatial  rate 
of  decay. 

For  correlations  along  an  axis  parallel  to  the  Jet  center- 
line,  the  first  zero  crossing  of  the  correlation  coefficient 
is  given  by 


«...  - V'"'  • (50) 

where  U is  the  convection  velocity  of  the  pressure  field  along 
c 

the  axis  of  correlation.  For  the  aerodynamic  region,  U is 

c 

equal  to  the  flow  velocity,  and  for  the  acoustic  region 
is  given  by 

U_  • e^/cos8 
o o 


(51) 


where  0 Is  the  angle  between  the  direction  of  propagation  of 
the  Incident  acoustic  wave  and  the  axis  of  correlation.  Be- 
cause of  the  implicit  assumption  of  an  incident  plane  wave, 
this  Is  essentially  a farfleld  approximation. 

In  the  region  dominated  by  the  aerodynamic  pressure  field, 
Clarkson  [65]  assumes  that  U is  proportional  to  f^,  so  that 

C 

the  first  zero  crossing  is  proportional  to  f . He  obtained 
values  of  n between  0.37  and  0.50  for  fullscale  engines  and 
between  0.6I  and  0.71  for  model  scale  jets.  For  the  acoustic 
field,  Clarkson  takes  n=0.  However,  this  latter  value  is  not 
strictly  correct,  since  0 In  equation  (51)  is  frequency- 
dependent  , 

In  order  to  demonstrate  the  contribution  of  0,  the  results 
of  Interpreting  the  data  from  [30]  in  terms  of  equation 
(51)  are  Indicated  by  hatched  regions  in  Figure  ^0,  which 
regions  indicate  the  corresponding  range  of  angles  of  propaga- 
tion. In  all  cases  the  higher-frequency  components  come  from 
areas  in  the  flow  that  are  closer  to  the  nozzle  exit  plane. 
(This  is  typical  of  the  noise  source  distribution  within  a jet 
exhaust,  but  the  locations  of  the  noise  sources  estimated  on 
the  basis  of  correlation  data  are  not  in  close  agreement  with 
more  directly  measured  source  locations,  probably  because  of 
the  simplifications  Incorporated  In  the  correlation  model.) 

The  results  of  an  analysis  by  Clarkson  [65],  consisting  of 
application  of  equations  (50)  and  (5I)  to  the  measured 
first  zero  crossings,  are  shown  in  Figure  ^1.  Figure  *12 
shows  the  results  of  another  analysis,  where  the  angle  0 Is 
determined  by  source  locations  obtained  from  the  proposed  model 
curve  of  Figure  16.  The  theoretical  value  of  0.25  for 
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I 1 NEAR  FIELD  MEASUREMENTS 


FIG.  41.  EFFECT  OF  ANGLE  9 ON  FIRST  ZERO  CROSSING  OF 

ACOUSTIC  PRESSURE  CORRELATION  COEFFICIENT  [65]. 


X16 


Ks«; 


FIG.  42.  COMPARISON  OF  NON-DIMENSIONAL  VALUES  FOR  PRESSURE  CORRELATION 
FIRST  ZERO-CROSSING  LENGTH. 


fC  (cos  9)/c  shown  In  this  figure  follows  from  equation 
1 > 0 o 

(50).  It  is  seen  that  the  experimental  values  here  tend  to 
be  higher  than  the  predicted  ones,  particularly  for  small 
values  of  x/D  and  y/D). 

In  relation  to  the  correlation  decay  in  the  longitudinal 
direction,  if  one  assumes  that  the  correlation  coefficient  can 
be  written  as 

(-m  oj5,cos  e \ ro)5,cos  O'] 

c“ ) I c J (52) 

then  the  experimental  results  show  that  m^  lies  between  O.O6 
and  0.73,  with  an  average  value  of  0.23.  The  value  of  m^  de- 
pends to  some  extent  on  the  filter  bandwidth,  but  the  varia- 
tion is  small  and  can  be  neglected  for  practical  purposes. 

The  highest  value  of  m^  was  obtained  close  to  the  nozzle  exit; 
however,  the  data  show  no  consistent  trend  with  either  dis- 
tance from  the  nozzle  exit  or  angle  from  the  Jet  axis.  Thus, 
it  is  proposed  that  the  average  value  of  0.23  be  used  in  the 
preliminary  model. 

Correlation  data  for  the  lateral  or  circumferential  di- 
rections are  very  sparse.  Some  measurements  were  made  by 
Howes  at  al  [30]  in  a direction  normal  to  the  plane  pass- 
ing through  the  Jet  centerline.  The  resulting  correlation 
coefficients  showed  characteristics  similar  to  those  for  the 
longitudinal  direction.  However,  it  is  more  difficult  here  ' o 
relate  the  first  zero  crossing  of  the  correlation  to  acoustic 
wave  propagation  parameters. 
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If  plane  wave  acoustics  (farfleld)  are  applied  to  the 
data  obtained  by  Howes  et  at  [30]  for  x/D  = 17.3  and  y/D 
■ 3.^9,  one  finds  the  first  zero  crossing  r.  to  be  given 
approximately  by 


sin 


X 


0.06 


(53) 


corresponding  to  the  correlation  measured  along  an  axis  per- 
pendicular to  a plane  through  the  Jet  axis;  x denotes  the 
angle  between  this  plane  and  a typical  ray  from  the  effective 
acoustic  source  to  the  mean  measurement  location. 


When  one  compares  equation  (53)  with 


^ ^1,0  Q 


0.25 


(5^) 


which  is  obtained  from  equations  (50)  and  (51) > one  notes 
that  the  correlation  lengtn  in  the  lateral  direction  i's  shorter 
than  that  predicted  on  the  basis  of  farfleld  plane  wave  propa- 
gation. Part  of  this  discrepancy  is  probably  due  to  the  fact 
that  the  noise  at  a given  frequency  does  not  come  from  a single 
coherent  source,  but  from  a spatial  distribution  of  uncorrelated 
sources.  (The  circumferential  distribution  of  the  sources  will 
have  a greater  influence  on  the  lateral  correlation  than  will 
the  axial  distribution  on  the  longitudinal  correlation.) 


The  lateral  correlation  decay  index  m^,  which  is  used  In 
association  with  the  function  (i>C(sln  x)/c^,  in  analogy  to  equa- 
tion (52)  is  found  to  have  a value  of  1.21  on  the  basis  of 
the  very  limited  data  of  Howes  et  aZ  [30]. 
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2.9.3  Surface  Pressure  Data 


The  prec'edlng  discussion  referred  to  correlation  measure- 
ments under  freefield  conditions.  Several  investigators  have 
also  measured  pressure  correlation  coefficients  on  flat  or 
curved  surfaces,  with  diverse  results. 

For  the  case  of  a single  surface,  such  as  a panel,  data 
of  Howes  et  at  [30]  and  Cox  et  at  [l6]  show  that  the 
correlation  coefficients  are  similar  in  form  to  thone  obtained 
in  the  freefield,  and  that  the  distance  to  the  first  zero  cross- 
ing and  the  decay  rate  also  are  similar.  Figure  ^3  contains 
typical  data  from  [16],  comparing  pressure  correlations  mea- 
sured in  the  freefield,  on  a rigid  panel,  end  on  a flexible 
structure  of  skin-stringer  configuration,  1th  a skin  thick- 
ness of  0.0^”.  The  data  show  no  significant  change  due  to  the 
presence  of  either  panel. 

The  situation  is  more  complicated  vrhere  there  are  two  sur- 
faces, such  as  near  the  Intersection  of  the  vertical  and  hori- 
zontal stabilizers.  Clarkson  [65]  presents  data  measured 
on  a configuration  of  this  type,  for  both  model  and  fullscale 
systems.  Here  the  reflected  pressure  field  resulting  from  the 
second  surface  Influences  the  shape  of  the  spatial  correlation 
and  the  distance  to  the  first  zero  crossing  point  (see  Figure 
In  certain  frequency  ranges,  the  distance  to  the  first 
zero  crossing  is  constant  or  increases  with  frequency,  in  con- 
trast to  the  freefield  case,  where  distance  to  the  first  zero 
crossing  decreases  as  frequency  increases.  Clarkson  constructs 
a simple  model  for  noise  sources  in  the  Jet  to  describe  the  re- 
flection effect,  but  the  results  show  only  qualitative  agreement 
with  the  measurements.  The  presently  proposed  prediction  model 
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FIG.  44.  EFFECT  OF  REFLECTION  FROM  PERPENDICULAR  PLANE  ON 
ZERO  CROSSING  POINT  OF  CORRELATION  COEFFICIENT  IN 
LONGITUDINAL  DIRECTION  [65]. 


does  not  account  for  the  effects  of  reflecting  surfaces.  The 
correlation  coefficient  Is  assumed  to  have  a form  which  Is 


separable  In  the  coordinates  and 

and  is  taken 

P “ exp 

-0.23wC^cos  0 

exp 

-1.210)5 ^sln  x' 

c^ 

c^ 

0 

0 

m m 

o>5  cos  6 
1 

cos 

i|.17  0)5  sin  X 
s 

c^ 

c^ 

0 

o 

(55) 


A factor  of  ^.17  has  been  Introduced  In  the  terms  Involving 
X»  In  order  to  provide  agreement  with  the  small  amount  of  ex- 
perimental data;  It  Is  an  artifact  which  has  no  physical  inter- 
pretation at  the  present.  Obviously  there  is  a need  for  further 
understanding  and  for  improvement  of  the  pressure  correlation 
model . 


2.10  List  of  Symbols  for  Sec.  II* 


Symbol 

Definition 

Units 

t 

A 

nozzle  area 

ft* 

(m*) 

D 

nozzle  diameter 

ft 

(m) 

equivalent  diameter 

ft 

(m) 

^h 

hydraulic  diameter 

ft 

(m) 

Jn 

Bessel  function  of  first 
kind 

— 

N 

Mach  number 

— 

Mach  number  of  aircraft 

— 

^n 

exit  Mach  number  of  Jet 

Bessel  function  of  second 
kind 

— 

OASPL 

overall  sound  pressure  level 

dB 

PR 

nozzle  pressure  ratio 

— 

R 

nozzle  pressure  ratio 

«C 

critical  nozzle  pressure 
r*vtlo 

— 

S 

Strouhal  number 

— 

SPL 

sound  i-  ' assure  level 

dB 

T 

absolute  temperature 

(°K) 

total  temperature  of  Jet 

»R 

n) 

"a 

temperature  of  2unblent  air 

®R 

CK) 

Jet  temperature 

®R 

(°K) 

reference  temperature 

°R 

(‘^K) 

% 

convection  ve’  ^Ity 

ft/sec 

(m/s) 

V 

flow  o>ioclt.v 

ft/sec 

(m/s) 

^a 

speed  of  aircraft 

ft/sec 

(m/s) 

•Specially  defined  symbols;  ■"ch  as  those  for  empirically  deter- 
mined constants,  which  ar«  -ed  only  once  In  the  text^  are  not 
Included  here.  Such  symb  ..s  are  defined  In  the  text  where  they 
occur. 

tfnie  units  given  here  are  typical  ones.  SI  units  are  given  In 
parentheses  where  appropriate.  Note  that  some  empirical  predlc- 
tlhn  methods  require  the  use  of  specific  units. 
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Definition 


Units 


Jet  velocity 

ft/sec 

(m/s) 

refer  mce  velocity 

ft/sec 

(m/s) 

sound  power  per  unit 
volume 

(lb/ft*sec ) 

(N/m*s) 

axial  distance 

ft 

(m) 

source  location  along  axis 

ft 

(m) 

lateral  distance 

ft 

(m) 

sound  speed  In  Interna- 
tional standard  atmosphere 

ft/sec 

(m/s) 

sound  speed  In  air  at 
sea  level 

ft/sec 

(m/s) 

speed  of  sound  In  ambient 
air 

ft/sec 

(m/c) 

frequency 

Hs 

thickness  of  annulus 

ft 

(m) 

local  scale  of  Jet 
turbulence 

ft 

(m) 

length  of  potential  core 

ft 

(m) 

sound  pressure 

lb /ft* 

(N/m* ) 

free-fleld  sound  pressure 

Ib/f t * 

(N/m*) 

source-to-observer  distance 

ft 

(m) 

characteristic  length  of 
shock  cell 

ft 

(m) 

fluctuating  velocity 

ft/sec 

(m/s) 

coordinates 

ft 

(m) 

Impingement  angle  (Fig.  37) 

deg 

engine  bypass  ratio 

— 

angle  fron  Jet  axis 

deg 

separation  distance 

ft 

(m) 

Symbol 

Definition 

Uni  ts 

p 

density  of  Jet  flow 

slug/ft* 

(kg/m*) 

^ISA 

density  of  International 
standard  atmosphere 

slug/ft  * 

(kg/m*) 

PSL 

air  density  at  sea  level 

slug/ft  * 

(kg/m* ) 

Pe 

equivalent  flow  density 

slug/ft’ 

(kg/m* ) 

P# 

density  of  ambient 
atmosphere 

alug/ft * 

(kg/m* ) 

p(e,»c,,T) 

pressure  correlation 
coefficient 

T 

time  delay 

sec 

(s) 

angle  In  nozzle  exit 
plane 

deg 

<i) 

radian  frequency 

sec"  * 

(s"M 

Wo 

radian  frequ-^ncy  at 
spectrum  peak 

sec"  * 

(s“* ) 

reduced  soundpressure 
level 

dB 

SECTION  III 


POWERED  LIFT  DEVICES 

3.1  Introduction 

Powered  lift  devices,  which  have  been  the  subject  of 
active  study  In  only  the  past  few  years,  serve  to  Increase 
lift  coefficients  and  thus  to  reduce  landing  and  takeoff 
distances. 

Powered  lift  or  ’’propulsive  lift”  systems  are  characterized 
by  considerable  integration  of  the  propulsion  and  flap 
systems.  Therefore,  the  noise  of  powered  lift  systems  tends  to 
depend  on  more  parameters  than  propulsion  noise  by  itself. 
Because  of  this  complexity  and  the  relatively  newness  of  powered 
lift  concepts,  the  available  noise  prediction  methods  are 
largely  based  on  empirical  data  based  on  model  studies  j 
fortunately,  the  scaling  laws  are  known  well  enough  to  permit 
extrapolation  of  model  data  to  full  scale. 

3.1.1  Overview  of  powered  lift  systems 

Propulsive  lift  systems  fall  into  two  basic  categories, 
called  "externally  blown”  and  "internally  blown”  configurations. 
In  externally  blown  systems,  the  lift-producing  flov;  produced 
by  the  propulsion  system  Is  entirely  outside  of  the  flap 
structures;  In  Internally  blown  configurations  this  flow  begins 
inside  a portion  of  the  flaps. 

Externally  blown  configurations  are  of  two  types:  (1) 

under-the-wlng  externally  blown  flaps,  where  the  Jet  exhaust 
Impinges  directly  on  the  flaps  and  Is  turned  downward  by  them, 
and  (2)  upper  surface  blown  flaps,  where  the  Jet  exhaust  stream 
Is  directed  along  the  upper  surface  of  the  wing  and  the 
attached  flow  is  turned  downward  by  (the  upper  surfaces  of)  a 


series  of  flaps.  The  Internally  blown  configurations  also 
have  developed  Into  two  classes:  (1)  Internally  blown  flaps » 

which  are  basically  Jet  flaps,  where  high-pressure  engine 
air  Is  ducted  to  a slot  nozzle  blowing  over  a downward  directed 
flap,  and  (2)augmenter  wings,  which  essentially  use  a Jet 
flap  with  a second  flap  (or  shroud)  above  It,  so  as  to  form  a 
channel  that  guides  the  slot  nozzle  flow  and  entrains 
additional  flow  from  the  upper  wing  surface.  Numerous 
variations  and  hybrids  of  these  propulsive  lift  system 
classes  have  also  been  proposed  and  investigated  to  various 
extents . 

3.1.2  General  noise  mechanisms 

All  powered  lift  systems  Involve  the  Impingement  of 
high-speed  flow  on  flaps  and  other  parts  of  the  airframe 
structure.  This  Impingement  produces  noise  as  the  result  of 
the  several  separate  acoustic  mechanisms  summarized  below. 

1.  Body  Sources.  Lift  and  drag  fluctuations  that  act 
in  phase  over  a substantial  portion  of  the  surface 
of  a flap  radiate  sound  like  an  acoustic  dipole 

at  low  frequencies  (for  which  the  acoustic  wave- 
length Is  gre'iter  than  the  flap  chord). 

2.  Edge  Sources.  Pressure  fluctuations  of  small 
spatial  extenti  Interacting  with  the  leading  or 
trailing  edges  of  a flap,  radiate  sound  like  span- 
wise  distributions  of  partially  baffled  dipole 
sources.  Edge  sources  are  Important  primarily  at 
frequencies  at  which  the  acoustic  wavelength  is 
snaller  than  the  flap  chord. 
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3.  Surface  Sources.  Small-scale  pressure  fluctuations, 
such  as  present  In  a turbulent  shear  layer  on  a flap 
surface  or  In  a region  where  flow  Interacts  with  surface 
discontinuities  or  protuberances,  radiate  sound  directly. 

The  aforementioned  acoustic  sources  result  from  fluid- 
dynamic  disturbances,  which  may  be  separated  into  the  following 
three  categories  for  the  caae:s  of  interest  here; 

1.  In  flow  Turbulence » such  as  that  arising  from  the 
shear  layer  at  a Jet  edge  or  from  internal  engine 
turbulence,  may  produce  lift  and  drag  fluctuations 
over  the  entire  flap  surface  and  also  tends  to 
contain  small-scale  fluctuations  in  pressure  and 
velocity,  which  interact  with  the  leading  edge. 

2.  Boundary  Layer  and  Attached  Wall  Jet  Flow  leads  to 
edge  noise  generation  as  it  passes  over  the  trailing 
edge  and  also  produces  surface  sources. 

3.  The  Turbulent  Wake  of  an  airfoil  can  produce  whole- 
body  lift  and  drag  fluctuations  associated  with 
large-scale  eddies,  as  well'vas  aex’odynamlc  interactions 

i 

of  smaller-scale  eddies  with  the  trailing  edge.  In 
multiple-flap  systems  or  airfoil  cascades,  the  wake 
from  one  flap  also  produces  inflow  turbulence  for 
the  next  one  downstream. 

3.1.3  Acoustic  source  characterization 

The  predominant  sources  of  noise  produced  by  flow/surface 
Interaction  have  the  character  of  acoustic  dipoles.  The  acoustic 
intensity  I In  the  far  field  of  such  an  Interaction  source  obeys 
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I(«) 


I 

\ pc*  r* 


cos*e 


dS 


(56) 


where  w denotes  the  radian  frequency,  ♦p(w)  the  local  aero- 
dynamic force  spectrum,  6 the  angle  between  the  dipole  axis 
and  the  source-to-observer  line,  r the  source-to-observer 
distance,  p the  density,  and  c the  speed  of  sound  in  air. 

For  distributed  sources,  4>p,0,  and  r differ  for  each 
component  on  a surface  Sj^,and  the  net  effect  must  be  obtained 
by  integration  over  this  surface. 


For  evaluation  of  Equation  (56),  the  force  spectrum  dis- 
tribution must  be  known.  Since  this  information  is  not  available 
in  practice,  direct  evaluation  of  Equation  (56)  cannot  be  accom- 
plished. However,  the  form  of  this  equation  has  served  as 
the  basis  for  scaling  laws  that  have  been  useful  for  extra- 
polating model  data  to  full-scale  conditions.  Table  6 
summarizes  the  scaling  relationships  for  the  various  dipole- 
like  source  mechanisms  that  have  been  investigated  experimen- 
tally. 

The  directivity  of  an  isolated  flap  element  source 
depends  on  the  mechanism,  as  shown  schematically  in  Figure  ^5. 

Figure  46  illustrates  the  near-field  directivity  of  a 
trailing  edge  source.  Note  that  the  sound  is  maximum  along 
the  surface.  In  a plane  that  is  normal  to  the  surface,  the 
experimentally  observed  directivity  differs  from  that  of  an 
ideal  dipole  because  of  scattering  and  finlte-source-slze 
effects. 
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TABLE  6.  PARAMETRIC  DEPENDENCES  OF  FLOW/SURFACE  - 
INTERACTION  NOISE  MECHANISMS  [67]. 


Body  Source : I 

Trailing  or 

Leading  Edge  Dipole:  I 

« <o*F*  cos'^0 
« U®A^cos*e 
« nw*F*  sin®6  cos*(t|)/2) 

« W sln*e  cos*(V2) 

1 

I - 

1 

acoustic  intensity 

1 w * 

radian  frequency 

p2  . 

mean-square  fluctuating  force 

u - 

mean  velocity  of  flow 

A_» 

CO 

correlation  area 

n « 

effective  number  of 

sources  on  surface 

h " 

spanwlse  eddy  scale 

w = 

wetted  span 

See  Figure  ^15  for  definition  of  angles  e and 

(a)  Ip  plane  nornal  to  (b)  In  plane  normal  to  (c)  In  plane  of  the  surface 

surface  and  to  span  surface  and  parallel 

to  span 

FIG.  46.  DIRECTIVITY  OF  €DGE  SOURCES  [88]. 


3.2  Nedrfield  Noise  from  Spatially  D1 stributed Sources 


In  powered  lift  devices,  the  noise  sources  typically 
are  spread  over  a significant  portion  of  the  wing  span.  Thus, 
only  observation  points  at  considerable  distances  can  be 
expected  to  be  in  the  geometric  farfleld,  where  the  source  can 
be  approximated  by  a single  point. 

It  is  Instructive  to  compare  the  results  one  obtains  by 
assuming  farfleld  (point  source)  conditions  with  similar 
results  for  the  nearfleld  of  a line  source  made  up  of  a linear 
array  of  incoherently  radiating  point  sources.  Such  a com- 
parison is  made  in  Figure  47,  from  which  one  may  observe 
that  the  point  source  (farfleld)  and  line  source  (nearfield) 
results  are  appreciably  different  only  if  the  distance  r from 
the  center  of  the  line  source  to  the  observation  point  is 
smaller  than  the  length  w of  the  lihe  source. 

The  foregoing  observation  suggest  that  one  may  calculate 
the  sound  pressure  levels  in  the  geometric  nearfield  of 
distributed  sources  simply  by  subdividing  each  source  into 
elements  that  are  short  compared  to  the  observation  point 
distance  (i.e.,  w < r),  calculating  the  contribution  of  each 
such  element  as  if  it  were  a point  source,  and  then  combining 
the  results  to  account  for  all  elements. 

3.3  Externally  Blown  Configurations 

Both  the  under-the-wing  externally  blown  flap  (EBF) 
and  the  over-the-wlng  upper  surface  blown  (USB)  configuration 
produce  low-frequency  noise  due  to  flow  interaction  with 
turning  flaps,  and  both  may  be  expected  to  expose  the  fuselage 
to  some  Jet  flow.  As  shown  schematically  in  Figure  48,  the  EBP 
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results  In  a lateral  spreading  of  the  Jet  flow,  Increasing 
the  likelihood  of  Its  Impingement  on  the  fuselage  at  low 
forward  speeds;  In  the  USB  configuration  the  Induced  flow 
Is  directed  Inward,  leading  to  the  possibility  of  having 
highly  turbulent  shear  layers  Impinge  on  the  fuselage.  Thus, 
flow  excitation  must  be  considered  In  addition  to  acoustic 
excitation. 

3.3.1  Externally  blown  flaps  (EBP) 

General 

EBP  configurations  typically  involve  high  bypass  ratio 
engines,  arranged  under  the  wing,  blowing  on  two  or  three 
enlarged  flaps.  In  such  configurations  one  encounters  noise 
due  to  the  engine  and  Its  free  exhaust  Jet,  in  addition  to 
the  noise  due  to  sources  associated  with  interaction  of  flow 
with  the  flaps.  The  present  discussion  focuses  on  the 
latter  noise  sources;  however,  flap  noise  Is  caused  by  pressure 
fluctuations  in  the  flow  that  reaches  the  flaps,  and  this 
Inflow  is  due  to  the  Jet.  Indeed,  a variety  of  mixer  or 
"decayer”  nozzles  - e.g.,  see  [68]  - have  been  considered 
for  reducing  Inflow  turbulence  and  flap  noise. 

Flap  noiee 

Figure  ^9  Illustrates  some  typical  noise  spectra  measured 
on  a large-scale  EBP  model.  The  shape  of  the  spectrum  measured 
at  the  wing  tip  corresponds  to  what  one  may  expect  to  observe 
on  the  fuselage. 


FIG.  49,  SOUND  SPFCTRA  .lEASURED  AT  500  FT.  WINGTIP 
SIDELINE  AND  MAXIMUM  FOR  500  FT.  FLYOVER. 
THREE-FLAP  EBF  MODEL  WITH  SINGLE  TF-34 
ENGINE  (6000  LB.  MAXIMUM  THRUST)  [68] 
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A prediction  method  has  been  synthesized  from  two 
methods  previously  reported  by  Hayden  et  iil.  [69,  70]  and 
Clark  et  al.  [71].  This  suggested  method  Is  summarized  below. 


The  overall  sound  pressure  level  resulting  from  an 
individual  Jet  nozzle/flap  Interaction  is  given  by 

OASPL  a 29  + O.lil  + 10  log  (p!  U®  D*) 


+ 10  log 


’F 

1-^slntti 


cos 


(^)1  - = 


(51) 


0 log  r 


where  OASPL  = overall  Lound  pressure  level,  (dB,\re 

2x10*"  * N/m 

p_  « mean  density  of  Jet  (slugs/ft*)  \ 

“ effective  velocity  of  Jet  (ft/sec),  \ 

© 

discussed  below. 

Do  “ effective  nozzle  diameter  (ft);  discussed 
below 

r a distance  from  flow  impingement  point  to 
observer  (ft) 

4p  a flap  angle  (degrees);  see  Figure  50 
e,y»  a observation  point  coordinate  angles 
(degrees);  see  Figure  50. 


The  effective  velocity  U-  of  the  Jet  at  the  impingement 
point  depends  on  the  nozzle  type  and  on  the  axial  distance  x 
from  the  Jet  exit  plane  to  the  impingement  point.  The  variation 
of  Jet  velocity  with  axial  distance  is  illustrated  in  Figure  51 
for  several  types  of  nozzles.  Note  that  the  exhaust  gas 
velocity  from  conical  nozzles  decreases  only  slightly  for  the 
first  few  diameter's  distance,  whereas  the  velocity  from  decayer 
nozzles  decreases  markedly  within  two  diameters. 
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51.  VELOCITY  DECAY  CHARACTERISTICS  OF 
TYPICAL  EBF  HOZZLES  [68]. 
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The  first  flap  of  typical  two  engine  EBF  configurations 

Is  between  1.0  and  3.0  nozzle  diameters  from  the  exit  plane, 

% 

whereas  the  first  flap  of  four-engine  configurations  Is 
between  2.0  and  diameters  from  the  exit  plane.  Thus, 
all  but  decayer  nozzles  can  be  approximated  by  distance-inde- 
pendent Jet  velocity  expressions  for  the  EBP  application  of 
Interest  here. 

For  nozzles  other  than  external  mixers  (decayers),  one 

may  calculate  the  effective  jet  velocity  U from  the  following 

expressions  (which  apply  for  .x/D  < 5): 

0 ~ 

For  simple  round  nozzles, 

Ug  ■ Uq  " mean  exit  velocity  (58a) 


(58b) 


(58c) 

» area  of  core  exit  nozzle 
c 

A|.  » area  of  bypass  nozzle  annulus 
» core  exit  velocity 
» fan  (annulus)  exit  velocity 
BPR”  bypass  ratio 


For  (d  annular  nozzles  (without  internal  mixing), 
..  I*c  ''f  1 

• I j 

For  nozzles  with  Internally  mixed  flows, 

BPR  • 

II  _ L c 

^e  BFr^  1 

where 


I 


U2 


For  decayer  nozzles,  one  may  estimate  the  effective  Jet 
velocity  at  a given  x/D^  by  using  the  appropriate  ratio  of 
decayer  velocity  to  the  conical  Jet  velocity  derived  from 
Figure  51. 

The  effective  nozzle  diameter  is  equal  to  the  actual 
nozzle  exit  diameter  where  there  exists  a single  Jet.  For 
coannular  Jets, 

■ V?  ^*0  * (59) 


The  one-third  octave  band  sound  pressure  levels 
associated  with  the  foregoing  overall  levels  (for  nozzles 
other  than  decayers)  are  given  by 


= OASPL  - 7.5  + 10  log  SVS®'*  + 

i./  OWD  P P 


(60) 


where 

S “ strouhal  number  at  spectrum  peak 
**  0.3  for  typical  EBF  configurations 

S ■ fD^/Ug  ■ Strouhal  number  at  center  frequency 
f (Hz)  of  band 

The  relation  of  Equation  (60)  has  been  found  to  coofrespond  within 
at  most  i 3 dB  to  data  for  all  azimuth  angles.  For  decayer 
nozzles,  the  1/3-octave  band  spectrum  Is  nearly  flat  (wl+-hln 
i 2 dD;  see  Figure  52)  for  such  nozzles. 


where 
C - 


SPI^X/30B  ■ OASPL  - 11  + C , (for  0.1  < S < 4.0)  , (6l) 


j 60  log  [U 


for  x/D,  between  3 and  5 


“e.deeayer  **’«  8®“  velocity  of  a deoayer,  and 
is  the  gas  velocity  of  a conical  nozzle  at  the  x/D^  value  of 
interest » as  found  from  Figure  51. 

Unsteady  Aepodynamia  Pressures 

The  sources  of  fluctuating  pressures  on  blown  flaps  are 
summarized  in  Figure  53*  Figure  54  presents  the  results  of 
a survey  of  flap  pressures  carried  out  on  a large-scale 
EBF  model,  which  indicates  the  greater  spatial  variability 
that  may  be  expected  in  general. 

Figure  55  indicates  the  distribution  of  overall 
fluctuating  pressure  levels  measured  on  an  EBP  configuration 
in  a plane  through  jet  center  line.  Considerable  variability 
is  again  evident.  The  pressure  data  may  be  normalized  with 
respect  to  the  dynamic  pressure  of  the  jet,  and  spectra  may 
be  related  to  the  corresponding  Strouhal  number.  Figures 
56  and  57  show  some  typical  reduced  spectra.  The  latter 
figure  and  related  extensive  measurements  show  that  the  reduced 
spectrum  shape  is  essentially  the  same  under  a wide  variety 
of  conditions. 
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(a)  CONICAL  NOZZLE. 


FIG.  SZ.  NORMALIZED  SPECTRA  AS  A FUNCTION  Of  STROUHAL 
NUMBER  FOR  FLAP  DEFLECTION  OF  10®-20°. 
DIRECTIVITY  ANGLE,  100®j  U^.  53  B/sec  [72]. 
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FIG.  53.  SOURCES  OF  FLUCTUATING  PRESSURE 
ON  EXTERNALLY  iiLOWN  FLAPS  [92]. 
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SMALL-SCALE  EBF  MODEL  [73]. 


Ik7 


FIG.  55  OVERALL  LEVELS  OF  FLUCTUATING 
SLTIFACE  PRESSURES  ON  TF-34 
EBF  MODEL  (75) . 
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Correlation  length  Information  is  sparse;  only  spanwlse 
correlations  have  been  reported  [76].  A corresponding 
reduced  parameter  curve  suitable  for  estimation  purposes 
appears  in  Figure  58. 

3.3.2  Upper  surface  blown  flaps  (USB) 

General 

In  USB  systems,  flow  turning  is  produced  by  the  Coanda 
effect.  Complex  flow  fields,  as  illustrated  in  Figure  59,  are 
associated  with  such  systems.  The  flow  fields  are  characterized 
by  wall-jet  flows  along  the  Jet  exit  direction  and  a pair  of 
vortices  on  the  edges  of  the  wall  Jet.  Fluctuations  of  these 
vortices  produce  a low-frequency  peak  in  the  radiated  noise, 
as  well  as  in  the  fluctuating  pressure  spectra. 

Flap  Noise 

In  addition  to  the  aforementioned  low-frequency  peak,  the 
spectrum  of  flap  noise  contains  a peak  at  higher  frequencies 
associated  with  the  attached  flow.  As  the  nozzle  exit  velocity 
Increases,  the  high-frequency  peak  becomes  predominant,  as 
illustrated  in  Figure  60.  [67,  70,  77,  78], 

The  presence  of  two  peaks  complicates  the  prediction  of 
' USB  noise;  the  low-frequency  dipole  force  varies  as  the  fourth 
power  of  nozzle  exit  velocity  (due  to  the  finite  size  of  the 
surface),  whereas  the  dipole  forces  responsible  for  the  high- 
frequency  part  of  the  spectrum  vary  as  the  sixth  power  of  this 
velocity  [703. 
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FIG.  59.  FLOW  FIELDS  ON  UPPER  SURFACE  BLOWN 
CONFIGURATIONS  [92]. 
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FIG.  60.  EFFECT  OF  JET  VELOCITY  ON  NOISE 
OF  USB  CONFIGURATION  [78]. 


(AR4  NOZZLE,  NOZZLE  IMPINGEMENT  ANGLE  20°, 
NOZZLE  LOCATION  20%  CHORD,  30°  FLAP  ANGLE) 


Nearfietd  Noise  Data 

Figure  61  Illustrates  the  considerable  variations  In 
magnitude  and  spectrum  shape  that  have  been  observed  even  over 
small  areas  on  fuselage  sidewalls.  The  dependences  of  the 
corresponding  overall  fluctuating  pressure  levels  on  Jet 
velocity  are  shown  In  Figure  62;  the  fact  that  these  dependences 
exhibit  different  slopes  - corresponding  to  different  powers  of 
velocity  - Indicates  that  different  mechanisms  predominate. 


Figures  63  and  6k  show  some  detailed  data  obtained 
from  measurements  on  a USB  model.  Figure  63  indicates  the 
measurement  (microphone)  locations  on  a simulated  fuselage 
wall; Figure  6^  shows  the  spectra  observed  at  these  locations 
at  a particular  Jet  exit  speed.  One  may  note  that  the  high- 
frequency  portions  of  all  these  spectra  are  comparable,  but 
not  the  low-frequency  portions.  A thorough  study  of  this 
and  related  data  [80]  revealed  that  different  portions  of 
the  various  spectra  varied  differently  with  Jet  velocity, 
as  shown  in  Figures  65  and  66. 

Accurate  estimation  of  the  nearfleld  noise  requires  one 
to  understand  and  predict  all  of  the  contributing  noise  sources. 
Such  understanding  can  be  developed  for  specific  configurations 
only  on  the  basis  of  extensive  measurements;  a general  prediction 
technique  applicable  to  all  configurations  appears  to  be  beyond 
the  present  state  of  the  art. 

In  absence  of  such  a general  technique,  use  of  Figure  67 
Is  recommended  for  spectrum  prediction  purposes,  with  the 
overall  sound  pressure  level  estimated  by  application  of  Equa- 
tion ( 57 )i.  modified  by  replacing  0.14  6p  by  0.016p.  Figure  67 
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FIG.  64.  NEARFIELD  NOISE  LEVELS  ON  USB  MODEL  FUSELAGE 
SIDEWALL  AT  680  fps.  JET  VELOCITY  [SOJ." 
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NEARFIELD  NOISE  SCALING  FOR  SENSOR  LOCATION  PS#1  [80]. 


applies  for  a nominal  nozsle  aspect  ratio  of  ^ for  a nominal 
jet  Impingement  angle  of  20°;  the  effects  of  these  pc.rameters 
are  discussed  below. 

Effects  of  Nozzle  Shape  and  Impingement  Angle 

The  shape  of  the  nozzle  and  the  "klckdown”  angle  at  which 
the  jet  impinges  on  the  upper  surface  of  the  wing  have  been 
shown  to  influence  the  USB  noise  spectrum  peak  [78,  79].  Fig- 
ure 68  shows  noise  spectra  (in  terms  of  Strouhal  number, 
referred  to  L^,  the  axial  distance  from  the  nozzle  exit  plane 
to  the  flap  trailing  edge)  for  nozzles  with  several  aspect 
ratios  and  the  same  Impingement  angle;  Figure  69  shows  similar 
spectra  for  the  same  aspect  ratio  and  several  impingement 
angles . 

For  nozzle  aspect  ratios  other  than  ^.0  and  impingement 
angles  other  than  20°,  one  may  use  the  differences  between 
corresponding  spectra  shown  in  Figures  68  and  69  to  obtain 
corresponding  approximate  corrections  to  the  general 
estimated  spectrum  of  .Figure  67. 

If  a deflector  is  used,  its  effect  on  the  spectrum  may 
similarly  be  determined  by  use  of  Figure  70. 

Unsteady  Aerodynamic  Pressursa 

m 

Flap  pressures  on  USB  configurations  have  been  found  to 
vary  considerably  with  location  on  the  flap,  nozzle  shape, 
and  exit  velocity.  Figure  71  summarizes  some  data  taken 
on  the  USAP  YC-1<I  USB  transport  aircraft.  The  observed  low- 
frequency  peak  that  occurs  at  pressure  ratios  above  1.36  is  as 
yet  not  understood. 
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RE  LEVEL  (dB)  SOUND  PRESSURE  LEVEL  (dB) 


LEVEL  (dB) 


FIQ.  71. 


TYPICAL  KING  SURFACE  PRESSURE  DATA 
OBTAINED  ON  YC-14  AIRCRAFT  (31.5®/54® 
FLAP  SETTING)  FOR  SEVERAL  NOZZLE 
PRESSURE  RATIOS  [81]. 


Reduced  data  from  a full-scale  USB  static  tests  are 
shown  In  Figure  72.  This  reduced  data,  when  normalized  in 
terms  of  Strouhal  number,  may  be  used  for  general  first- 
approximation  estimates.  Figure  73  indicates  correlation 
length,  information. 

3.4  Internally  Blown  Configurations 

Both  the  augmenter  wing  and  the  Jet  flap  involve  the 
ducting  of  high-pressure  air  to  a slot  nozzle  or  to  a series 
of  small  nozzles  spanning  essentially  the  entire  wing  trailing 
edge.  Because  of  the  associated  large  spatial  distribution 
of  noise  sources,  the  geometric  nearfleid  of  Internally  blown 
configurations  is  quite  extensive. 

3.4.1  Augmenter  wing  (AW) 

Noiee 


For  an  augmenter  wing,  a typical  cross-section  of  which 
is  shown  schematically  in  Figure  7^,  noise  generation  can  be 
treated  in  two  frequency  regions  [82]: 

1.  At  high  frequencies.  Jet  noise  components  are 
generated  within  the  augmenter,  and  their  noise 
levels  typically  are  15  dB  below  those  due  to 
fluctuating  aerodynamic  pressures  on  the  flap  and 
shroud.  (However,  this  Jet  noise  may  dominate 
farfleld  noise  levels  if  the  augmenter  Is  not 
treated  acoustically.) 


1€6 


2. 


At  low  frequencies,  jet  noise  Is  generated  outside 
the  augmenter.  This  noise  and  that  from  Interaction 
of  flow  with  the  trailing  edgeo  of  the  flap  and 
shroud  both  make  s.-* ,’nlflcant  contributions. 

Some  augmenter  wings  have  been  observed  also  to  produce 
Intense  "screech"  tones  due  to  an  aeroacourjtlc  feedback  action 
Involving  the  nozzle  flow  Impinging  on  the  leading  edge  of 
the  flap  [83].  Because  of  the  relative  rarity  of  this 
phenomenon  and  a lack  of  understanding  of  it.  It  Is  not  treated 
further  In  this  report. 

The  sound  field  from  an  augmenter  has  a complex  spectrum 
and  directivity,  as  will  be  illustrated  here  for  a particular 
design  that  has  been  studied  extensively.  This  design,  v'hich 
Is  shown  in  section  in  Figure  75»  was  provided  with  a screech 
suppressor  consisting  of  flutes  (about  0.1  In.  wide  and  2.0  in. 
long,  spaced  1.75  In.  apart)  cut  into  the  nozzle  edge  near 
the  flap.  Some  observed  farfield  noise  spectra  and  the  far- 
fleld  directivity  are  shown  In  Figures  76  to  78,  Prom 
these  figures  one  may  deduce,  for  example,  that  In  the  143° 
direction  (where  there  occurs  a general  directivity  peak),  noise 
due  to  the  slit  Jet  predominates,  whereas  In  the  82°  direction, 
low-frequency  noise  from  the  flaps  and  augmenter  exhaust 
predominates. 

The  best  available  prediction  stethod  [86j  gives  a some- 
what different  dependence  of  noise  on  velocity  than  that 
observed  for  the  previously  discussed  configuration  [83].  This 
prediction  method  results  In  levels  that  are  too  high  (by  no 
more  than  10  dB),  and  thus  conservatively  overestimates  the 
aeroacoustlc  environment  for  design  purposes. 
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FIG.  72.  SURFACE  PRESSURES  HEASURED  IN  STATiC 
TEST  OF  FULL-SCALE  USB  MODEL  [93]. 
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FIG.  76.  ACOUSTIC  FAR  FIELD  SPECTRA  FOR  AUGHEHTER  WING 
OF  FIGURE  75  AT  0 FORWARD  SPEED  [85]. 
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FIG.  77.  FAR  FIELD  ACOUSTIC  SPECTRA  FOR  MODEL  NOZZLE 
KITH  AND  UITHCUT  AUGMENTOR  WINS  [85]. 


3 


Surfaae  Preeauvea 

The  only  data  on  unsteady  aerodynamic  pressures  on  the 
flap  and  shroud  that  are  available  to  date  are  those  of  [85]. 
Some  of  this  data  Is  reproduced  here  in  Figures  79  to  8l.  In 
absence  of  better  information,  this  data,  normalized  with 
respect  to  the  dynamic  pressure  of  the  Jet,  may  be  used  for 
predictive  estimates. 

Data  on  correlation  lengths  is  extremely  limited;  what 
there  exists  exhibits  such  unusual  behavior  that  the  data’s 
validity  is  questionable  [85].  Thus,  oriy  vex*y  gross 
estimates  can  be  made. 

3.4.2  Jet  flap  (JF) 

The  Jet  flap,  of  which  a schematic  sketch  appears  in 
Figure  82,  has  been  studied  much  less  extensively  from 
the  noise  standpoint  than  have  the  three  other  powered  lift 
concepts,  lha  most  definitive  available  study,  which  was 
published  in  1959  [87],  provided  evidence  that  two  sources 

of  noise  are  Important:  (1)  Jet  noise  from  the  slot  nozzle, 

radiating  into  the  hemisphere  above  the  wing  only,  and 
(2)  sources  associated  with  the  flap  (now  known  to  be 
’’trailing  edge  noise”). 

One  may  therefore  synthesize  an  estimation  approach  by 
accounting  separately  for  Jet  noise  and  for  flap  noise, 
treating  the  latter  in  terms  of  a two-dimensional  wall-jet 
[88]. 


No  data  is  available  on  surface  pressures.  For  first 
order  estimates,  similar  augmenter  wing  data  [85]  and  related 
wall  Jet  data  [893  should  be  useful. 
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FIG.  81.  NON-DIMENSIONAL  PRESSURE  SPECTRA  AT  DIFFERENT 
CHORDWISE  LOCATIONS  ON  SURFACE  OF  FLAP 
SEE  FIGURE  75  [85}. 
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3.5  Forwar(<  Speed  Effects 


Because^of  the  complex,  multi-component  nature  of  the 

noise  of  propulsive  lift  devices,  the  evaluation  of  the 

% 

effects  of  freeatream  velocity  (aircraft  motion)  on  propulsive 
lift  noise  are  extremely  difficult  to  predict  analytically. 
However,  generalized  empirical  relations  have  been  derived 
[86]  on  the  basis  of  studies  performed  by  several  investi- 
gators. 


These  empirical  relations  give  the  change  A dB  in  the 
sound  pressure  level  due  to  fon»ard  motion  at  speed  as 


A dB 


10  log 


10  log 


1 - 


1 - 


a 


K 


V cos6 

8i 


K 


for  EBP,  USB,  AW 


for  JP 


where  Vj  denotes  the  Jet  exhaust  velocity  and  5 the  flap 
angle  (see  Figure  82).  The  exponent  K is  a function  of  the 
powered  lift  configuration  and  observation  position,  as 
given  in  Table  7. 


(62) 
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3.6  List  of  Symbols 

for  Section  III* 

Symbol 

Definition 

Unlts^ 

Ac 

area  of  core  exit  nozzle 

ft*  (m*) 

Aco 

correlation  area 

ft*  (m*) 

Af 

area  of  bypass  nozzle 
annulus 

ft*  (m*) 

BPR 

bypass  ratio 

effective  nozzle 

ft  (m) 

0 

diameter 

p* 

mean-square  fluctuat- 
ing force 

lb*  (N*) 

I(w) 

acoustic  Intensity 

Ib/ft  sec  (N/ms) 

OASPL 

overall  sound  pressure 
level 

dB 

S 

Strouhal  number 

SPI. 

sound  pressure  level 

dB 

surface  area 

ft*  (m*) 

effective  velocity  of 

ft /sec  (m/s) 

jet 

• \ 

\ 

speed  of  aircraft 

ft/sec  (m.'s) 

core  exit  velocity 

ft/sec  (m/s) 

Vf 

fan  annulus  exit 
velocity 

ft/sec  (m/s) 

w 

wetted  span 

ft  (m) 

'Specially  defined  syr*iJola,  auch  as  those  for  empirically  oeter- 
mined  constants,  which  are  used  only  once  in  the  text,  are  not 
Included  here.  Such  symbols  are  defined  in  the  text  where  they 
occur. 

4* 

The  units  given  hero  are  typical.  SI  units  are  gi’-  . ' in  paren- 
theses where  appropriate.  Note  that  some  empirical  prediction 
methods  require  the  use  of  specific  units. 
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Symbol 

Definition 

Units 

c 

speed  of  sound  In 
ambient  air 

ft/bec  (m/s) 

t 

correlation  length 

ft  (m) 

ly 

n 

spanwlse  eddy  scale 
number  of  sources 

ft  (m) 

r 

source-to-observer 

distance 

ft  (m) 

w 

width  of  line  source 

ft  (m) 

X 

axial  distance 

ft  (m) 

«P 

flap  angle 

deg 

6 

angle  (see  Eq.  56, 
Pigs.  M5,1»7) 

deg 

P 

density  of  ambient  air 

slug/ft*  (kg/m 

Pc 

mean  density  of  Jet 

slug/ft*  (kg/m 

♦p(w) 

aerodynamic  force 
spectrum 

IbVft*  (N*/m* 

angle  (see  Figs.  ^5,50) 

deg 

(D 

radian  frequency 

sec“*  (s“*) 

SECTION  IV 
PROPELLER  NOISE 


4.1  Introduction 

In  propeller-driven  aircraft,  the  propellers  frequently 
constitute  the  dominant  sources  of  farfleld  and  nearfleld 
noise.  Propeller  noise  typically  Is  considered  In  terms  of 
two  categories:  (1)  rotational  noise  Includes  all  tonal 

(single-frequency)  components  at  Integer  multiples  of  the  blade 
passage  frequency,  and  (2)  broadband  or  vortex  noise  Includes 
all  of  the  non-tonal  noises. 

4.2.  Noise  Generation  Mechanisms 

4.2.1  Rotational  Noise 

Rotational  noise,  which  usually  dominates  the  overall  sound 
pressure  level  and  Is  most  pronounced  at  low  frequencies.  Is 
due  to;  (1)  rotation  of  the  steady  forces  acting  on  the  pro- 
peller blade,  as  the  blade  generates  thrust  and  torque  and  as 
the  blade  displaces  air  by  virtue  of  Its  thickness,  and  (2)  the 
unsteady  thrust  and  torque  (or  "harmonic  loading”)  acting  on 
the  blades  as  the  result  of  spatial  variation  of  the  Inflow  to 
the  propeller  disk.  The  steady  force  sources  of  rotational 
noise  have  received  the  most  attention  to  date.  The  noise  due 
to  harmonic  loading  has  been  studied  less  intensively,  and  gen- 
erally not  enough  Is  known  about  the  inflow  nonuniformities  that 
cause  this  type  of  noise. 

4.2.2  Broadband  Noise 

Broadband  noise  Is  principally  caused  by  (1)  dipole  radia- 
tion due  to  the  turbulent  boundary  layer  flowing  off  the  blade 
trailing  edge  and  (2)  turbulence  In  the  Inflow  to  the  blade 


leading  edge.  The  latter  Is  negligible  In  standard  tractor 
propeller  Installations. 

4.3  Extent  of  the  Nearfield 

The  geometric  nearfield  of  a propeller  Is  the  region  In 
which  the  propeller  disk  can  not  he  approximated  as  a point 
source.  The  acoustic  near  field  is  the  region  in  which  the 
"inverse-square”  spherical  spreading  law  does  not  hold  for*  the 
variation  of  the  pressure  with  distance  from  a source.  As  a 
consequence  of  the  large  size  of  aircraft  propellers  and  of 
the  low  frequencies  at  which  they  generate  noise,  the  nearfield 
of  aircraft  propellers  generally  is  quite  large. 

4.3.1  Rotational  Noise 

Figure  83  Illustrates  the  difference  between  sound  pres- 
sure levels  calculated  on  the  basis  of  a farfield  approximation 
and  those  obtained  by  considering  the  acoustic  and  geometric 
nearfields,  for  a particular  propeller.  One  may  observe  that 
here  the  nearfield  is  significant  at  distances  that  are  less 
than  1.75  diameters  from  the  center  of  the  propeller  disk.  For 
typical  propeller  aircraft  configurations,  the  fuselage,  nacelle, 
and  part  of  the  wing  are  likely  to  be  exposed  to  the  nearfield. 

4.3.2  Broadband  Noise 

None  of  the  available  literature  appears  to  treat  the  sub- 
ject of  the  nearfleld  region  for  broadband  noise.  Because  the 
broadband  nrise  is  significant  principally  at  middle  and  high 
frequencies,  the  Importance  of  its  nearfield  is  confined  to  lo- 
cations that  are  considerably  closer  to  the  propeller  disk  than 
those  where  rotational  nearfield  becomes  significant. 

Because  of  the  noiwcoherent  nature  of  broadband  noise,  there 
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are  no  distance-related  or  pngle-related  Interference  patterns 
to  distinguish  a nearfleld.  However,  one  expects  a geometric 
nearfleld  effect,  because  of  the  location  of  the  sound  sources 
within  the  propeller  disk. 

4.4  Rotational  Noise  Analysis 

Since  forces  that  cause  rotational  noise  are  distributed 
over  the  surface  of  the  propeller  blades,  calculation  of  the 
sound  generated  by  the  propeller  requires  Integration  of  the 
effects  of  these  forces  over  the  entire  propel) ei*  disk,  both 
In  the  radial  and  the  circumferential  directions.  At  observa- 
tion points  located  at  large  distances  from  the  propeller,  the 
distances  to  all  points  In  the  propeller  disk  may  be  considered 
as  equal.  Then  the  farfleld  circumferential  Integral  may  be 
simplified  to  yield  the  Bessel  function  directivity  derived  by 
Gutln  ClOO],  eliminating  the  requirement  to  recompute  this 
Integral  In  each  case.  A further  approximation  to  the  farfleld, 
which  has'^been  used  by  Gutln  and  by  Garrick  and  Watkins  [lOl], 
among  others,  makes  use  of  the  "effective  ring"  concept.  In 
which  all  propeller  forces  are  assumed  to  lie  In  a ring  at  a 
single  radius  - usually  taken  as  0.8  times  the  propeller  radius. 
This  approximation  eliminates  the  radial  Integration. 

In  the  nearfleld,  the  assumptions  governing  the  approxima- 
tion leading  to  elimination  of  the  circumferential  Integral  cer- 
tainly do  not  hold;  path  lengths  and  angles  to  a llstenlnc'  point 
from  different  points  on  the  propeller  disk  may  vary  by  large 
percentages.  However,  the  effective  radius  approximation  Is 
workable  for  the  nearfleld.  Thus» the  calculation  of  the  near- 
fleld noise  of  a propeller  In  each  case  requires  at  least  evalu- 
ation of  the  circumferential  Integral  over  the  disk. 
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4.5  Prediction  Methods 


ii 


4.5.1  Rotational  Noise 

The  nearfleld  noise  prediction  technique  of  Kerwin  and 
Pranken  [ 34  ] was  developed  by  combining  the  Gutin  noise  model 
[100]  with  the  static  propeller  noise  measurements  of  Hv.bbard 
[102],  This  recommended  prediction  method,  which  is  applicable 
at  distances  of  one  propeller  diameter  or  less  from  the  edge 
of  the  propeller  disk,  is  accurate  only  for  positions  outside 
the  axial  projection  of  the  propeller  disk. 

The  most  recently  developed  prediction  method  [105]  is 
not  recommended,  because  its  results  differ  drastically  from 
those  of  earlier  ones;  in  effect,  this  method  is  based  on  the 
assumption  that  the  noise  due  to  torque  is  negligible. 

At  distances  greater  than  one  propeller  diameter  from  the 
disk,  the  farfleld  prediction  method  from  Reference  [103 ] is 
recommended. 

4.5.2  Broadband  Noise 

The  most  useful  prediction  method  here  is  that  of  Refer- 
ence [104]  this  method  is  a theoretical/empirical  hybrid  and 
differs  from  earlier  ones  in  that  it  predicts  finite  levels  at 
points  in  the  plane  of  the  propeller,  where  the  earlier  methods 
predict  a null,  at  variance  with  test  data.  This  difference  is 
obviously  of  great  significance  for  predicting  noise  on  the 
fuselage  sides  adjacent  to  the  propellers.  This  prediction 
method  is  valid  for  the  farfleld,  and  tends  to  result  in  some 
underestimation  at  the  pressures  very  near  the  propeller  in  the 
propeller  plane,  and  in  overestimation  of  pressures  near  the 
propeller  axis  at  distances  that  are  within  about  0.7  of  the 
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propeller  radius  from  the  edge  of  the  propeller  disk. 

4.5.3  Reflective  Surfaces 

The  effect  of  flat  and  cylindrical  reflective  surfaces  on 
the  nearfield  pressures  has  been  studied  by  Hubbard  [102]  and 
used  in  the  nearfleld  prediction  technique  of  [103] 

4.5.4  Spectral  Coherence  and  Correlation 

Past  work  on  propeller  noise  prediction  has  generally  over- 
looked the  subject  of  spectral  coherence  and  correlation 
patterns.  The  generation  mechanisms  for  propeller  rotational 
noise,  hpwever,  allow  for  making  some  simple  observations, 
which  are  summarized  below. 

The  rotational  noise  from  the  thrust,  torque,  and  thickness 
effects  has  a phase  pattern  which  is  fixed  relative  to  the  pro- 
peller blades  and  rotates  with  them.  The  phase  speeds  in  the 
radial  and  axial  directions  are  sonic,  but  a doppler  correction 
must  be  applied  to  the  wave  length  and  to  the  phase  speed  as 
measured  along  the  aircraft  longitudinal  axis. 

The  acoustic  field  of  harmonic  loading  noise  differs  from 
that  from  thrust,  torque,  and  thickness,  because  its  pattern 
consists  of  standing  waves  that  are  determined  by  the  inflow 
distribution.  Otherwise,  all  of  the  wavelengths  and  phase 
speeds  are  the  same  as  for  the  rotating  sound  field. 
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SECTION  V 
SURFACE  FLOWS 


5.1  Pressure  Fluctuations  In  Turbulent  Boundary  Layers 

5.1.1  Introduction 

Boundary  layers  develop  as  a consequence  of  the  effect 
of  viscosity.  These  boundary  layers  become  turbulent  where 
the  Reynolds  number,  based  on  the  distance  from  the  point  of 
streamline  attachment,  becomes  sufficiently  large.  This  crit- 
ical Reynolds  number  depends  on  the  local  pressure  gradient; 
for  zero  pressure  gradient  it  typically  is  of  the  order  of 
5x10^,  but  it  can  be  smaller  or  greater  by  a factor  of  about 
100  for  non-zero  gradients. 

The  structure  of  a turbulent  boundary  layer  is  affected 
not  only  by  the  basic  fluid  properties  (e.g.,  viscosity, 
density)  and  the  Reynolds  number  and  Mach  number,  but  also  by 
the  surface  roughness,  the  pressure  gradient,  and  the  velocity 
field  outside  of  the  boundary  layer.  For  prediction  of  the 
surface  pressure  fluctuations  associated  with  a given  boundary 
layer,  one  needs  to  know  certain  of  its  local  properties,  par- 
ticularly the  appropriate  thickness.  Unfortunately,  these 
properties  are  difficult  to  predict,  because  they  depend  on 
the  upstream  history  of  the  boundary  layer. 

Current  prediction  methods  are  based  largely  on  experimen- 
tal data  obtained  on  a limited  number  of  configurations,  and 
for  a limited  range  of  parameters.  Thus,  the  available  empir- 
ical data  often  may  not  correspond  well  to  the  case  for  which 
a prediction  is  desired. 
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5.1.2  boundary  Layer  Thickness  Parameters 

Two  different  thickness  parameters  are  used  to  nondimen- 
slonalize  boundary  layer  data:  the  boundary  layer  thickness 

6 and  the  boundary  layer  displacement  thickness  6*. 

The  boundary  layer  thickness  is  defined  as  the  distance, 
measured  from  the  surface,  at  which  the  flow  velocity  reaches 
99X  of  its  local  free-stream  value. 

The  boundary  layer  displacement  thickness  is  defined  by 

r - X"  [l  - (63) 

and  represents  the  distance  for  which  pU6  equals  the  mass 
flow  deficit  produced  by  the  presence  of  the  boundary  layer 
(as  compared  to  the  mass  flow  based  on  the  free  stream  velo- 
city). Here  U denotes  the  free-stream  velocity,  p the  fluid 
density  in  the  free  stream,  and  u(y)  and  p(y)  denote  the 
corresponding  properties  in  the  boundary  layer,  at  distance  y 
from  the  solid  bounding  surfaces.  Clearly,  6 depends  on  the 
velocity  and  density  profiles  in  the  boundary  layer.  For  an 
incompressible  turbulent  boundary  layer  on  a flat  plate, 

6 « 86*. 

5.1.3  RMS  Pressure  Fluctuations 

Data  indicating  the  overall  level  of  pressure  fluctuations 
beneath  developed  turbulent  boundary  layers  is  available  for 
several  configurations  (tunnel  models,  flight  tests ; flat  plates, 
cones,  aerodynamic  surfaces)  and  for  a broad  range  of  Mach  num- 
bers. As  Indicated  in  Figure  8l|,  the  data  falls  within  a fairly 
narrow  band,  if  it  is  nondlmensionallzed  with  respect  to  the 
local  free-stream  dynamic  pressure  q.  Several  investigators 
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LOCAL  FREE-STREAM 


FIG.  84.  RMS  PRESSURE  FLUCTUATI 
TURBULENT  BOUNDARY  LAY 
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have  suggested  oeml-emplrlcal  curves  or  expressions  to  fit 
the  data.  No  single  approach  appears  to  have  any  definite 
advantage  over  the  others.  The  form  suggested  by  Lowson 
[107]  for  M < 5>  and  a constant  for  M > 5 is  as  good  os  any; 
thus,  the  ratio  of  the  root-mean-square  fluctuating  pressure 
p to  the  dynamic  pressure  q may  be  taken  to  obey 


0.006 


E . 

q 


1+O.l^lM^ 


0.00133 


M < 5 
M > 5 


(6ila) 

(6i|b) 


At  low  subsonic  speeds,  this  expression  reduces  to  p/q  » 0.006, 
which  falls  near  the  middle  of  the  data  band  In  Figure  i|.l. 
Recent  careful  subsonic  measurements  with  very  small  pressure 
transducers  have  yielded  values  of  p/q* as  high  as  0.011, 
suggesting  that  In  many  earlier  measurements  the  high-frequency 
energy  has  been  filtered  out,  due  to  the  finite  transducer  size 
[1083. 


6.1.4  Power  Spectral  Density 

0 

Envelopes  of  normalized  experimental  data  on  the  power 
spectral  density  of  boundary  layer  surface  pressure  fluctua- 
tions are  shown  in  Figure  85.  Separate  envelopes  are  present- 
ed for  subsonic  and  supersonic  data,  based  on  the  summary  of 
experimental  data  given  by  Chaump,  et  al  [106].  The  subsonic 
data  corresponds  to  0.1<N<0.8  and  the  envelope  includes  the  ex- 
tensive flat-plate  data  of  Blake  [109 3 » which  was  not  included 

» 

In  t 106)«  The  supersonic  data  corresponds  to  l.^^Mj<3.^5i  the 
slender-cone  data  for  M ■ 4.0  and  8.0  [106]  falls  In  the 
middle  of  the  envelope. 


SUPERSONIC  DATA 


FORMALIZED  ONE-THIRD  OCTAVE  BAND  POWER  SPECTRAL  DENSITY  OF  PRESSURE 
FLUCTUATIONS; ATTACHED  TURBULENT  BOUNDARY  LAYER  FLOW. 


Both  tho  subsonic  and  the  superset  ''-‘■r  zi<:.  pi  wttencea 
hi'^e  only  In  <^iivt>xupe  form,  because  no  systematic  parametric 
trends  are  evident.  For  Instance,  there  Is  no  definite  sys- 
tematic relationship  between  Mach  number  and  where  the  data 
lie  within  the  envelope. 

The  simple  seml-emplrlcal  relation  of  Houbolt  [110], 

S(tt)U  . 2/v 

p 6»  1 + (w««/U)^ 

Is  plotted  in  Figure  85,  as  is  the  relation  of  Robertson 

[111], 

(66) 

which  la  based  on  supersonic  data.  The  two  equations  may  be 
seen  to  yield  similar  results  at  high  frequencies  and  to  agree 
reasonably  well  with  the  average  of  the  high-frequency  data. 

At  very  low  frequencies,  the  second  equation  gives  a value 
that  exceeds  that  of  the  first  equation  by  a factor  of  it,  at 
low  frequencies,  equation  (65)  represents  the  average  of  the 
subsonic  data  reasonably  well,  but  equation  (66)  agrees  more 
closely  with  the  supersonic  data. 

5.1.5  Convection  Velocity 

The  convection  velocity  Is  a measure  of  the  speed  at  which 
turbulent  pressure  disturbances  move  along  with  the  boundary 
layer  flow.  The  convection  velocity  depends  on  the  frequency 
and  spatial  scale  of  the  disturbance;  low-frequency  disturb- 
ances generally  correspond  to  large  scale,  and  high-frequency 
dlsturba.'^oes,  to  small  scale.  Uirge  scale  disturbances  tend  to 


travel  at  nearly  the  freo-r>tream  speed,  whereas  smal]  scale 
disturbances  (which  occur  primarily  very  near  the  wall)  travel 
more  slowly,  4ue  to  the  velocity  gradient  in  the  boundary  layer. 

Convection  velocity  usually  Is  evaluated  from  the  time 
delay  Indicated  b;;.  the  space-time  correlation  function  ob- 
tained from  pressure  measurements  made  by  sensors  located  at 
known  distances  apart.  Figures  86  shows  an  envelope  of  broad- 
band convection  velocity  data,  as  a function  of  streamwise 
separation  distance.  The  flattening  of  the  curve  at  large 
separation  distances  Is  due  to  smaller,  slower  eddies  decaying 
mere  rapidly  and  having  no  effect  at  sufficiently  large  separ- 
ation distances.  Figure  87  shows  the  variation  of  narrowband 
convection  velocity  with  frequency,  for  various  separation 
distances.  (These  figures  are  based  on  data  by  Chyu  and  Hanly 
[112]  and  Bull  [113],  as  cited  in  [106].) 

The  empirical  formula  given  by  Lowson  [10?],  when  modi- 
fied by  taking  6 » 86»,  gives  the  narrow-band  convection  velo- 
city as 

U^/U  - 0.75  + 0.3*  exp(-0.886«w/U)  - 0. 25  • exp(-0. 156/6* ) (67) 

corresponding  to  which  the  broadband  convection  velocity  is 
given  by 

Ueb/U  » 0.8  -0.25*  exp(-0. 156/6*).  (68) 

These  formulas  are  plotted  in  Figures  86  and  87,  together 
with  corresponding  data.  The  agreements  may  be  seen  to  be 
reasonable  in  general. 
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5.1.6  Broadband  Correlations 

Normalized  broadband  space-time  correlation  functions 
are  shown  In  Figures  88  and  89  for  various  streamwlse  and 
lateral  spatial  separations,  C and  n,  respectively.  The  data 
spans  a Mach  number  range  from  0.3  to  3-^5*  Analytical  for- 
mulas to  fit  this  data  are  not  available.  A time-delayed 
peak  occurs  only  In  the  streamwlse  correlations,  since  turbu- 
lent eddies  convect  in  that  direction  only. 

The  contours  of  broadband  spatial  correlation  at  zero 
time  delay  (t  = 0)  shown  in  Figure  90  reveal  an  elongated 
streamwlse  structure  associated  with  the  greater  length  scales 
In  that  direction.  This  figure  is  based  on  experimental 
results  [ iiH  ] obtained  from  subsonic  wind-tunnel  tests.  Data 
for  a wider  Mach  number  range  are  not  available  at  present; 
neither  are  analytical  expressions  for  the  contour  shapes. 

5.1.7  Narrowband  Correlations 

The  principal  data  on  narrowband  correlations  has  been 
summarized  by  Chaump,  et  al  [106],  A compilation  of  data 
on  the  narrowband  streeimwlse  and  lateral  spatial  cross- 
correlation coefficients  for  zero  time  delay  appears  in 
Figures  91  and  92,  respectively. 

Corresponding  narrowband  space-time  correlation  data  are 
shown  in  Figures  93  and  9^  respectively.  In  these  two 
figures,  the  results  presented  are  for  those  values  of  the  time 
delay  t which  maximize  the  correlation  for  the  given  value  of 
C or  n;  thus,  these  figures  essentially  show  the  decaylng- 
exponentlal  envelope  of  oscillatory  functions  that  appear  much 
like  tho-e  of  Figures  91  and  92. 
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REDUCED  LATERAL  SEPARATION  DISTANCE,  77 /S* 


FIG.  89.  NORMALIZED  LATERAL  BROAD  BAND  SPACE-TIME 
CORRELATION  FOR  ATTACHED  TURBULENT 
FLOW  Cl 06]. 


NORMALIZED  NARROW  BAND  CORRELATION  [R(0,  C»0)/R  {0,0,( 


FIG.  92.  FREQUENCY-DEPENDENT  NARROWBAND  LATERAL  SPATIAL 
CORRELATION  [106],  ' 
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FIG.  94.  LATERIAL  VARIATION  OF  NARROWBAND  SPACE-TIME 

CORRELATIONS  FOR  ATTACHED  TURBULENT  FLOW  [106]. 


Thi  data  of  Figures  93  and  may  be  seen  to  con- 
verge for  large  values  of  the  Strouhal  numbers  shown,  but 
to  differ  considerably  for  small  values  of  the  Strouhal  num- 
bers, l.e.,  for  low  frequencies.  The  asymptotic  low- 
frequency  values  are  plotted  In  Figure  95,  as  functions  of 
non-dlmensionallzed  separation  distance. 

The  space-time  correlation  and  spectral  density  functions 
comprise  a Fourier  transform  pair: 

/OD 

R(C,n,t)  e^“^dt  (69) 

2TrR(C,n,T)  «/  S(5,n,w)  e“^“^dw,  (70) 

.00 

and  S(C,  n,  m)  Is  often  expressed  in  terms  of  the  co-spec tral 
density  (Its  real  part)  and  the  quad-spectral  density,  (its 
imaginary  part ) : 

S(5,n,w)  = - J<^i(?,n,w).  (71) 


The  co-spectral  density  Is  generally  used  in  defining  character- 
istics of  the  fluctuating  pressure  field,  Lowson  [107  ] and 
Robertson  [llTl,  using  data  obtained  by  Bull  [113],  described 
the  cross-correlation  function  for  turbulent  flow,  assuming  that 
the  cross-power  spectral  density  Is  a function  of  only  the  sep- 
aration distances,  by 


♦rU,0*w)  / ♦^(?) 

/ ?^(n) 


A^(C,  «)  cos(u)C/0^) 


OT 


(72) 


A^(n»  m) 


(73) 


>*’‘v 


where  (?)  " [♦r  (C^)  and  (C^)  Is  evaluated 

at  one  measurement  location  and  (C  ) at  another,  a distance 
_ n 2 

C away.  tn)  is  analogously  defined.  The  correlation  co- 
efficients may  be  expressed  as  [II3] 

A^(C,  w)  - expC  O.iew/U^]  exp[-0.03^C/6»J  (7^) 

A^(n,  u)  = exp[-0.72nw/u^]*{0. 3 + 0.7  exp[-0. 5n/6*]}  , (75) 

where  C and  n are  always  taken  as  positive. 

These  coefficients  are  plotted  in  Figures  91-95*  from 
where  one  can  perceive  their  relatively  good  fit  to  the  data. 
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5.2  Pressure  Fluctuations  In  Transitional  Boundary  Layers 

5.2.1  Introduction 

Boundary  layer  transition  Is  the  process  by  which  a laminar 
boundary  layer  changes  Its  structure  to  become  a fully  turbulent 
boundary  layer.  The  transition  process  takes  place  over  a stream 
wise  distance  which,  for  a typical  aircraft,  represents  a small 
fraction  of  the  total  surface  area.  However,  transitional  flow 
may  cover  much  of  the  surface  of  a ballistic  re-entry  vehicle. 
Boundary  layer  transition  is  a sensitive  phenomenon  which  depends 
strongly  on  such  factors  as  Reynolds  number,  surface  roughness 
and  discontinuity,  and  the  pressure  gradients  to  which  the  bound- 
ary layer  has  been  subjected. 

In  a transitional  boundary  layer,  unstable  disturbances  are 
amplified  and  lead  to  the  onset  of  turbulence.  The  turbulence 
occurs  first  in  patches,  so  that  transitional  flow  becomes  a 
changing  mixture  of  turbulent  and  disturbed  laminar  flow  regions. 

Because  of  their  complexity  and  sensitivity,  transitional 
flov’s  are  particularly  difficult  to  characterize.  Furthermore, 
little  reliable  data  is  available  and  the  general  validity  of 
particular  test  results  is  often,  difficult  to  ascertain.  Thus, 
use  of  any  data  or  prediction  scheme  may  be  expected  to  involve 
considerable  uncertainty  for  any  specific  case  for  which  a pre- 
diction is  desired. 

I 

5.2.2  RMS  pressure  fluctuations 

In  many  cases,  particularly  those  involving  supersonic  and 
hypersonic  flow,  the  fluctuatlni,  pressure  levels  in  transitional 
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boundary  layers  have  been  observed  to  exceed,  by  considerable 
amounts  (as  much  as  10  to  20  dB),  corresponding  levels  for  fully 
developed  turbulent  boundary  layers.  Figure  96  shows  an 
envelope  of  available  data*  for  the  overall  rms  fluctuating 
pressure  level  In  transitional  boundary  layer  flows,  as  a function 
of  Mach  number.  Comparison  of  these  levels  with  Fig.  4.1,  which 
pertains  to  fully  developed  turbulent  flows,  shows  the  levels 
associated  with  transitional  flows  to  be  substantially  higher. 


Also  shown  In  Fig.  96  Is  a plot  of  the  equation 

£ ^ 0.006 

q 1 + 0.013M* 


(76) 


which  represents  a modification  of  an  expression  given  by  Chaump 
et  at,  [106]  to  fit  data  in  th<*  Mach  number  range  of  4 to  8. 
Equation  (76),  which  pertains  to  transitional  flow,  gives 
higher  rms  pressure  values  than  does  Eq.  (64a),  which  pertains 
to  turbulent  flow.  Equation  (76)  was  constructed  so  that  the 
values  predicted  by  it  for  low  Mach  numbers  agree  with  the  cor- 
responding turbulent  flow  values.  However -^%>ils  range  Is  out- 
side the  data  base  for  transitional  flows,  and  use  of  the 
equation  may  lead  to  considerable  error. 


5.2.3  Power  spectral  density 

A power  spectrum  representative  of  supersonic  transitional 
flow  is  shown  in  Fig.  76.  This  spectrum  is  described  by 


*The  data  envelope  of  Fig.  96  Includes  the  data  of  [106,  115, 
116],  ell  taken  on  conical  models  in  wind  tunnels. 


I 
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NORMALIZED  POWER  SPECTRAL  OENSiTY  FLUCTUATING  PRESSURE  LEVEL 
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FIG.  96.  RMS  PRESSURE  FLUCTUATION  LEVEL  OF 
TRANSITIONAL  BOUNDARY  LAYER  FLOWS 
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3(m)U  . 2/ff 

p*  1 + (w«^/U)* 


which  relation  was  suggerced  by  Chaump  et  at.  [106].  Equation 
(120)  is  identical  to  Eq.  (65)*  except  for  the  transitional 
flow  length  pareuneter  6^  which  here  replaces  the  turbulent  bound- 
ary layer  displacemen-,  thickness  6*.  Because  6^  > 6*  In  general, 
transitional  spectrum  in  effect  is  shifted  toward'  lower  frequency 
with  respect  to  the  turbulent  boundary  layer  spectrum.  For  Mach 
numbers  between  3.7  and  8.1,  the  relation 


^ » 2.9  (1  + 0.013M®)* 


(78) 


has  been  suggested. 

It  should  be  mentioned  that  various  researchers  have  mea- 
sured  spectra  that  differ  considerably  from  that  given  by  Eq. 
(77);  some  spectra  show  a pronounced  low-frequency  peak. 

Note  also  that  the  data  on  which  Rq.  (77)  is  based  lies 
within  an  envelope  that  is  up  to  6 dB  wide,  depending  on  the 
observation  location  on  the  vehicle. 


5.2.4  Correlation  coefficients  and  convection  velocities 

Reliable  values  for  correlation  coefficients  in  transitional 
flow  are  not  available.  The  previously  discussed  formulae  for 
fully  turbulent  flow  may  be  used,  but  are  likely  to  be  seriously 
in  error  in  acnae  oases. 


the 
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Broadband  convection  velocities  have  been  found  to  be  ap- 
proximately 0.6  of  those  in  turbulent  flow.  This  trend  is 
physically  cdnsistent  with  the  previously  indicated  fact  that 
6^  > 6». 


5,3  Pressure  Fluctuations  in  Separated  Flows 
5.3.1  Introduction 

The  phenomenon  of  flow  separation  occurs  when  a boundary 
layer  detaches  from  the  surface,  resulting  in  a region  of  flow 
reversal  which  is  often  characterized  by  unsteadiness  and  the 
formation  of  large  turbulent  eddies  or  vortices.  Separation 
often  occurs  as  the  result  of  a loss  of  energy  from  the  boundary 
layer,  usually  due  to  an  adverse  pressure  gradient,  or  as  the 
result  of  a viscous  boundary  layer’s  inability  to  turn  a very 
sharp  corner.  Separation  leads  to  large  energy  losses  and 
attendant  drag  increases;  therefore,  designers  of  flight 
vehicles  take  care  to  avoid  its  occurrence.  Nevertheless,  in 
some  situations  separation  is  still  encountered. 

In  subsonic  flow,  separation  is  associated  with  either  an 
adverse  pressure  gradient  (due  to  body  shape)  or  a sharp  edge 
or  corner.  Common  areas  of  subsonic  flow  separation  are  the 
leading  or  trailing  edges  of  wings  (depending  on  the  airfoil 
configuration),  portions  of  flaps,  and  behind  bluff  bodies, 
such  as  protuberances,  struts,  or  landing  gear.  The  separation 
behind  bluff  bodies  is  termed  "base  flow",  and  is  discussed 
in  Sec.  fj. In  supersonic  flow,  the  occurrence  of  flow  separa- 
tion and  the  presence  of  shock  waves  are  often  closely  related, 
because  shock  waves  are  associated  with  body  shape  changes  and 
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also  produce  very  large  adverse  pressure  gradients.  Often  the 
region  of  flow  separation  and  the  associated  shock  wave  Interact 
dynamically  to  produce  a shock  oscillation  phenomenon  which 
produces  particularly  Intense  pressure  fluctuation  levels. 

The  remainder  of  this  section  discusses  data  and  prediction 
methods  for  several  phenomena  associated  with  flow  separation, 
first  for  the  very  Important  case  of  supersonic  flow,  and  later 
for  flow  over  wings.  Because  flow  separation  characteristics 
depend  strongly  on  surface  configurations  and  flow  conditions. 

It  must  be  noted  that  the  sug^^sted  prediction  methods  are 
limited  to  cases  that  correspond  to  the  data  on  which  they  are 
based.  In  some  situations.,  significant  deviations  from  predic- 
tions obtained  by  these  methods  may  be  expected. 


5.3.2  RMS  pressure  leve'^s 

♦ 

The  pressure  fluctuations  in  separated  flow  depend  on  the 
condition  of  flow  separation;  therefore,  several  sets  of  data 
need  to  be  presented.  Figure  98  shows  the  envelope  of  data 
for  separated  flow  induced  by  an  expansion,  such  as  occurs  at 
a slanted  rearward  step  or  at  the  junction  between  a cone  and 
cylinder.  The  data  on  which  the  envelope  Is  based  is  taken 
from  the  work  of  several  researchers,  as  summarized  by  Chaump 
et  al.  [IOC],  and  corresponds  to  expansion  angles  between  about 
15°  and  35°.  An  equation  to  predict  expansion-induced  separated 
flow  pressure  levels,  as  proposed  by  Robertson  [111],  is 


^rms 


0.0i»5 

1 + M* 

e 


(79) 


where  denotes  the  Mach  number  downstream  bf  the  expansion,  and 
q represents  the  upstream  dynamic  pressure.  As  evident  from  Pig. 
98,  Eq,  (79)  agrees  well  In  the  available  data. 
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FIG.  98.  OVERALL  FLUCTUATINfi  PRESSURE  LEVEL  IN 

REGIONS  OF  FLOW  SEPARATION  CAUSED  BY  AN 
EXPANSION  CORNER 
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FIG.  99.  OVERALL  PLUCTUATING  PRESSURE  LEVEL  IN 
REGIONS  OF  FLOW  SEPARATION  i RODUCFD  BY 
A COMPRESSION  CORNER 


Figure  99  shows  a data  envelope  for  separated  flow  up- 
stream of  a compression  corner  or  region  of  rapid  area  Increase 
of  a body.  The  available  data  [112,  117]  is  limited  to  com- 
pression corner  angles  of  ^5®.  According  to  [106 ]I,  an  analytical 
expression  has  not  been  developed,  because  insufficient  information 
is  available  on  local  flow  conditions.  The  assumption  of  a con- 
stant level  gives  acceptable  agreement  for  the  presently  available 
data : 


^rms 


0.022 


(80) 


At  the  point  on  a surface  where  flow  separation  begins,  the 
levels  tend  to  be  particularly  high  because  the  flow  alternates 
between  being  attached  and  separated.  In  the  case  of  supersonic 
flow,  a shock  wave  oscillation  may  be  associated  with  the  point 
of  Initial  separation.  Chaump  et  at,  [106],  on  the  basis  of  data 
from  [ill,  112],  have  suggested  the  following  expression  for  the 
overall  fluctuating  rms  pressure  at  a separation  point : 

Eee?-  . . (81) 

q«  1 + 0.5M* 

The  available  data  and  this  expression  are  indicated  In  Pig.  lOO. 

A 

Note  that  the  available  data  are  very  limited.  One  may 
expect,  for  instance,  that  the  actual  levels  in  the  transonic 
region  are  higher  than  those  predicted  by  Eq.  (8l)..  High 
levels  may  also  occur  In  the  region  of  flow  reattachment,  par- 
ticularly if  a shock  wave  Is  present.  Some  data  points  for 
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shock  reattachment  from  [112  and  118]  are  a3so  Indicated  In 
Pig.  100.  This  data  Is  consistent  with  the  levels  shown  In 
Figs,,  98  and  99  for  regions  of  separated  flow  In  the  same 


5.3.3  Power  spectral  density 

An  expression  for  the  power  spectrum,  normalized  with  respect 
to  dynamic  pressure,  was  suggested  by  Roberton  [ill]: 


_ _ _ 

f 

f 5 
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i 

0.17  f\Jej 

(82) 


where  6 denotes  the  upstream  boundary  layer  thickness,  and  is 
the  velocity  above  the  separation  region.  A comparison  of  this 
expression  with  corresponding  data  appears  in  Pig.  101. 


The  power  spectrum  normalized  to  the  rms  pressure  is 
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(83) 


An  alternative  formula,  based  on  a data  curve  fit,  has  been 
derived  by  Coe  and  Chyu  [119]  and  is  given  by 


S(f)  U 


^ - exp  C-8.09^-1.239x-0.259x*-0.090x*-0.01^x‘*-0,001x*] 
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where 


X « An 


(85) 


POWER  SPECTRUM  FOR  SEPARATED  FLOWS 


This  formula  and  the  associated  data  envelope  are  also  shown  In 
Pig.  101  for  comparison.  It  is  recommended,  however,  that  Eq. 

(82)  be  used  fcr  predictions  of  spectra  for  all  types  of  separa- 
tion, with  the  rms  pressure  determined  from  the  appropriate  case 
of  Sec.  .3.2. 

An  experimental  power  spectrum  [106]  for  shock  wave  oscil- 
lations with  separated  flow  induced  by  two  and  three  dimensional 
protuberances  is  shown  in  Pig.  102.  The  flow  disturbances  were 
provided  by  cylindrical  protuberances  at  approximately  M = 1.5 
and  by  a ^5®  wedge  at  M = 2.0. 

5.3.4  Narrowband  convection  velocities 

Available  data  show  narrowband  convection  velocities  in 
separated  flow  to  Vc»ry  from  about  0.2  times  the  freestream  velo- 
city at  low  wavenumbers  to. nearly  1.0  times  at  high  wavenumbers. 
Due  to  this  wide  variation,  broadband  convection  velocities  have 
considerably  less  meaning  in  separated  flow  than  for  attached 
turbulent  boundary  layer  flow. 

Some  data  for  M ■ 2.0  in  compression-comer  produced  separa- 
tion are  shown  in  Pig.  103.  A representative  curve  for  attached 
flow  is  also  shown  in  the  figure  for  comparison;  the  trend  is 
seen  to  be  quite  different.  A trend  similar  to  that  of  the  data 
of  Pig.  103  also  appears  in  the  data  of  Pig.  10^1',  where  frequency 
is  normalized  with  respect  to  the  distance  from  the  separation 
point  and  the  flow  velocity  outside  the  separation  region.  The 
behavior  in  the  low  Strouhal  number  region  is  dominated  by  the 
slow  reolrculatory  flow  within  the  separation  region;  the  behavior 
at  hlgd^er  Strouhal  numbers  shows  the  effect  of  the  outer  Higher 
speed  separated  boundary  layer  flow. 
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FIG.  104.  NARROW  SAND  CONVECTION  VELOCITY  VERSUS 
STROUHAL  NUMBER  BASED  ON  DISTANCE  FROM 

SEPARATION  POINT  [120]  ■ 
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S.3.5  Broadband  correlations 


Normalized  broadband  space-time  correlation  functions  for 
separated  flow  are  shown  In  Fig.  105  for  several  streamwldo 
distances.  These  should  be  coinpai'ed  with  the  corresponding 
functions  presented  earlier  for  turbulent  boundary  layer  flow. 

No  data  for  correlations  in  the  lateral  direction  are  available 
for  separated  flow.  No  expressions  have  been  developed  to  fit 
this  data. 

5.3.6  Narrowband  correlations 

The  principal  data  on  narrowband  correlations  have  been 
summarized  by  Chaump  et  at.  [106].  Data  for  the  streamwise 
spatial  cross-correlation  coefficients  for  zero  time  delay  are 
shown  in  Pig.  106.  Coherence  function  data  are  presented  in 
Pig.  107,  which  compares  functions  for  attached  and  separated 
flow.  Corresponding  analytical  expressions  for  the  cross-spectra 
coefficients,  as  defined  In  Sec.  5.1.7,  are  [112j; 

exp  [-0.75n]  for  f6/Ug  < 

\ ‘I 

exp  c-0.75n  (f6/6xlO"»U^)**’3  for  f60_ '>  6x10“’* 

^ (86) 

for  f«/Ug  < 6x10“ ® 
for  6x10“ * < f6/U^  < 6x10”* 
for  f6/U^  > 6x10”*  . 

(87) 

Figure  108  shows  the  degree  of  correlation  that  exists 
under  a detached  shock  wave  and  the  surrounding  attached  and 
separated  flow  region.  Weak  correlation  is  observed,  except  for 
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FIG.  lOS.  NORMALIZED  BROAD  BAND  STREAHUISE  SPACE  TIME  ^ 
tORRtLATION  FUNCTIONS  FOR  SEPARATED  FLOW  [112] 
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ARISC«  OF  COHERENCE  FUNCTIONS  FOR  ATTACHED 
SEPARATED  FLOW  [112] 
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FIG.  108.  CORRELATION  OF  PRESSURE  FLUCTUATIONS  BETWEEN 
SHOCK  WAVE  AND  ADJACENT  ATTACHED  AND 
SEPARATED  FLOW  REGIONS  [111] 
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the  region  in  close  proximity  to  the  shock.  Not  surprisingly, 
the  correlation  between  the  shock  region  and  the  upstream  attached 
flow  region  l-s  the  weakest. 

5.3.7  Pressure  spectra  for  the  flow  separation  on  wings  in 
subsonic  flow 

A wide  variety  of  flow  separation  phenomena  occur  on  wings 
in  subsonic  flow.  The  presence  of  flow  separation  on  wings  is 
often  associated  with  the  occurrence  of  flow  phenomena  wh^-ch 
lead  to  stalling  (loss  of  lift)  of  the  wing  if  the  angle  of 
attack  is  Increased  sufficiently.  A comprehensive  review  of 
Information  on  flow  separation  and  stall  on  airfoils  has  been 
made  by  Heller  et  al,  [121]. 

On  straight  wings  and  airfoils,  regions  of  flow  separation 
can  develop  in  two  different  ways  as  the  angle  of  attack  is 
Increased.  For  thicker  sections,  the  flow  separation  begins  near 
the  trailing  edge  and  moves  forward  as  the  angle  of  attack 
increases.  For  wings  with  thinner  airfoil  sections  and  sharp 
leading  edges,  a region  of  flow  separation  occurs  at  or  near 
the  leading  edge  and  reattaches  farther  back  along  the  wing 
chord.  As  the  angle  of  attack  increases,  the  reattachment 
point  moves  farther  back;  stall  occurs  when  the  flow  fails  to 
reattach.  Some  sections  may  also  exhibit  tooth  leading  and 
trailing  edge  separation  phenomena  simultaneously.  The  highest 
fluctuating  pressure  levels  occur  in  the  reattachment  region 
of  leading  edge  flow  separation.  Maximum  spectral  levels,  in 
this  case,  tend  to  be  roughly  10  dB  higher  than  the  naxlmum 
levels  associated  with  trailing  edge  flow  separation. 

In  an  experiment  using  an  airfoil,  Heller  et  at.  [121] 
were  able  to  collapse  the  pfeasure  spectra  at  the  reattachment 
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point  measured  under  different  conditions  into  a relatively 
narrow  envelope,  as  shovm  in  Pig.  109*  Peak  spectral  levels 
away  from  the  reattachment  point  were  as  much  as  20  dB  lower. 

Separation  regions  on  swept  wings  can  also  be  characterized 
according  to  whether  they  occur  In  the  leading  or  trailing  edge 
region.  Again,  leading  edge  separation  phenomeiTa  produce  the 
highest  fluctuating  pressure  levels.  On  swept  wings,  the  leading 
edge  separation  flow  has  a very  three-dimensional  character.  An 
Intense  vortex  forms  In  the  separation  x-egion,  becoming  larger 
near  the  wing  tip,  Spanwlse  flow  occurs  in  the  vortex,  so  that 
fluid  entering  the  vortex  region  spirals  outward  toward  the  wing 
tip.  The  presence  of  the  vortex  tends  to  have  the  effect  of 
bringing  some  orderly  structure  to  the  flow  separation  region. 

This  is  most  clearly  Illustrated  in  the  case  of  highly  swept 
sharp-edged  delta  wings.  In  this  case,  although  the  flow  separa- 
tion region  may  cover  much  of  the  upper  surface,  the  vortex  is 
so  Intense  and  well  ordered  that  the  flow  in  and  ben  ath  the 
vortex  is  essentially  steady.  However,  for  moderately  swept 
wings  (say  a sweep  angle  of  less  than  50**),  the  highest  fluctuating 
pressure  levels  occur  in  the  region  from  directly  beneath  the 
vortex  back  to  the  reattachraent  line.  Heller  et  at.  [121]  were 
able  to  collapse  fluctuating  pressure  data  taken  in  this  region 
into  a reasonably  narrow  envelope.  The  results,  which  are  based 
on  measurements  made  on  a 45^  swept-wlng  model,  Are  shown  in 
Fig.  110.  The  fact  that  spectrum  peak  occurs  at  higher  frequency 
(similarly,  nondlmenslonal  peak)  in  Pig.  110  than  in  Pig.  99 
can  be  attributed  to  the  ordering  effect  of  the  volrtex  which  tends 
to  diminish  the  quantity  of  large-scale,  low-frequency  disturbances. 
The  spectra  of  Figs.  109  and  liO  can  be  considered  to  provide 
bounds  on  the  effects  of  sweep  between  0*^  and  45°. 
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FIG.  no.  NORMALIZED  FLUCTUATING  PRESSURE  SPECTRUM  IN  mTEX- 
FLOW  REGIME  (on  45®  - SWEPT  WING)  [121] 
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5.4  Surface  Pressure  Fluctuations  In  Base  Flow  Regions 
5.4.1  Introduction 

The  region  of  flow  separation  behind  a blunt  termination  of 
a body  is  referred  to  as  base  flow.  A typical  base  flow  region 
is  characterized  by  flow  recirculation  and  turbulent  vortex 
shedding.  Base  flows  are  found  behind  blunt  bodies,  dive  brakes, 
spoilers,  and  the  blunted  rear  of  conical  re-entry  vehicles. 

The  data  and  prediction  formulae  reported  in  these  sections  are 
the  result  of  measurements  on  circular  bases  behind  otherwise 
streamlined  shapes. 


5.4.2  RMS  pressure  fluctuations 

Subsonic  data  for  base  pressure  fluctuations  were  obtained 
by  Eldred  [122]  for  the  blunt  base  of  a cylinder  with  a Bti'-eam- 
llned  nose.  High  Mach  number  measurements  at  the  base  of  a cone 
were  reported  by  Chaump  et  al.  [106].  A formula  relating  the 
rms  pressure  to  the  base  static  pressure  and  to  the  Mach 
number  Mj^  of  the  adjacent  flow  Just  outside  the  wake  was  devcloiJtti 
by  Houbolt  [123]  and  later  modified  by  Chaump  [106]. 


Figure  111  presents  the  rms  pressure  level  as  a function 
of  Mach  number.  Eldred  observed  that  levels  wei*e  higher  away 
from  the  centerline  of  the  cylinder  in  his  experiment;  his 
results  (for  M < 0.4),  which  may  be  summarized  by 


1.5  X 10“^ 

1 X 


at  65X  radius 


(88) 


at  the  centerline 


are  plotted  in  Fig.  Ill  on  the  basis  of  the  assumption  that 


q„  . - I PV»  = I YP^MJ  . (89) 

The  high  Mach  number  data  of  Cl06]  are  also  indicated  in  Fig. 
111.  It  is  evident  that  there  Is  a large  Mach  number  interval 
for  which  no  data  is  available;  predictions  for  this  interval 
cannot  be  made  with  high  confidence. 

The  prediction  formula  suggested  in  [123]  and  modified  in 
[106]  is 

■»rm«  . “.OlWg 

’’b  1 + O.OI)M“  ’ 

which  is  also  plotted  in  Pig.  HI* 

For  very  low  Mach  number,  Eq.  (90)  reduces  to 

Prins  * Ob  • 

which  is  in  general  agreement  with  Eldred's  results.  At  high 
Mach  numbers,  Eq.  (90)  approaches 

"rm*  ‘ Pb  • 

implying  that  the  rme  pressure  level  here  depends  only  on  the 
base  static  pressure. 

It  should  be  noted  that  the  values  of  pj^  and  must 
generally  be  determined  from  a knowledge  of  the  aerodynamics 


FIG.  111.  RMS  BASE  PRESSURE  FLUCTUATION  LEVELS 


of  the  particular  configuration.  The  assumption  that  Is 

appropriate  for  Eldred’s  experiment,  because  the  flow  is  subsonic 
and  the  base  region  is  at  the  end  of  a streamlined  body  with 
nearly  parallel  mean  flow  upstream  of  the  base.  For  more  blunt 
bodies  in  subsonic  flow,  and  all  bodies  in  supersonic  and  hyper- 
sonic flow,  this  assumption  could  lead  to  serious  errors. 

S.4.3  Power  spectral  density 

A formula  for  the  power  spectral  density,  suggested  in  [106] 
on  the  basis  of  [123],  Is 


S(f)Vjj 


(91) 


This  equation  is  plotted  in  Pig.  112,  together  with  the  low 
Mach  number  data  envelope  of  [106],  as  reported  in  [123],  and 
the  high  Ma-ch  number  data  envelopes  of  [106].  The  agreement  is 
seen  to  be  relatively  good,  considering  the  large  Mach  number 
differences  between  the  data  sets. 


5.4.4  Correlation  functions 

No  data  or  validated  procedures  are  available  that  are  suf- 
ficiently general  to  warrant  inclusions  here.  However,  a formu- 
lation for  the  cross  correlation  R(t)  applicable  to  the  base  of 
a re-entry  vehicle  was  developed  by  Houbolt  [123]  and  is  dis- 
cussed in  [1C6]. 


5.5  List  of  Symbols  for  Sec.  V* 


Symbol 

Definition 

Unitst 

M 

Mach  number 

R(5,n,t) 

space-time  correlation 
function 

IbVff* 

(NVmM 

S(w);S(f) 

power  spectral  density 

Ib^sec/ft** 

(N^s/m**) 

U 

free-stream  velocity 

ft/sec 

(m/s ) 

«c 

convection  velocity 

ft/sec 

(m/s) 

P 

rms  fluctuating  pressure 

Ib/ft* 

(N/m*) 

q 

dynamic  pressure 

Ib/ft* 

(N/m^ ) 

u(y) 

local  flow  speed 

ft/sec 

(m./s) 

y 

distance  from  solid 
boundary 

ft 

(m) 

boundary  layer  thickness 

ft 

(m) 

«* 

boundary  layer  displace- 
ment thickness 

ft 

(ns) 

«t 

transitional  flow  length 
perometer 

ft 

(m) 

C,n 

separation  distances  In 
axial  and  lateral  direc- 
tions 

ft 

(m) 

p 

air  density 

slug/ft* 

(kg/m*) 

T 

delay  time 

sec 

(s) 

♦(C,n.w) 

spectral  density  function 

Ib^sec/ft"* 

(N^s/m** ) 

(1) 

radian  frequency 

sec”* 

(s-») 

•Specially  defined  symbols,  such  as  those  for  empirically  deter- 
mined constants,  which  are  used  only  once  in  the  text,  are  not 
Included  here.  Such  symbols  are  defined  in  the  text  where  they 
occur. 

The  units  given  here  are  typical  ones.  SI  units  are  given  In 
parentheses  where  appropriate.  Note  that  some  empirical  predlc 
tlon  metho4a  require  the  use  of  specific  units. 
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SECTION  VI 


CAVITY  NOISE 

6.1  Introduction 

V 

High-speed  flow  over  cavities  or  cut-outs  in  aircraft 
surfaces  frequently  produces  an  Intense  acoustic  field  wltVi 
single-frequency  (tonal)  pressure  fluctuations.  The  fluctu- 
ating pressures  can  be  high  enough  to  jeopardize  the  integrity 
of  nearby  structural  components  or  sensitive  instruments,  and 
can  subject  the  areas  adjacent  to  the  cavity  mouth  to  Intense 
radiated’  sound  fields. 

A 

Plow-Induced  pressure  oscillations  in  cavities  have  been 
studied  extensively,  particularly  for  the  basic  rectangular 
cavity  configuration  shown  in  Figure  113«  This  configuration  is 
of  considerable  practical  Interest  itself,  but  also  represents 
geometries  that  do  not  deviate  too  greatly  from  the  basic  rec- 
tangular shape.  Presently  available  information  pei'mits  one 
to  make  reasonable  estimates  of  the  frequency,  levels  and  modal 
behavior  for  shallow  rectangular  cavities  (i.e.,  for  length-^ 
to-depth  ratios  greater  than  2)  in  flows  that  have  boundary 
layer  thicknesses  that  are  much  smaller  than  the  cavity  depth. 
The  rectangular  cavities  that  have  been  studied  to  date  tend 
to  respond  two-dimenslonaily ; the  cavity  width  seems  to  be  of 
little  importance.  Pressure  variations  do  occur  in  the  depth 
direction,  but  these  variations  usually  are  small. 

The  oscillation  process  occurs  because  of  two  factors: 

(1)  interaction  between  the  pressure  modes  within  the  cavity 
and  the  shear  layer  over  the  cavity  mputh;  and  (2>  interaction 
of  this  shear  layer  with  the  cavity  trailing  edge.  If  the 
shear  layer  that  impinges  on  the  cavity  trailing  edge  is  de- 
flected slightly,  then  there  results  a net  mass  addition  or 
removal  at  the  rear  of  the  cavity;  this  variable  ma.ss 


addition  has  an  effect  similar  to  that  of  replacine  the  rear 
bulkhead  with  a moving  piston,  which  excites  pressure  modes 
in  the  cavity.  Under  some  conditions,  the  shear  layer  motion 
and  the  trailing  edge  mass  addition  reinforce  each  other, 
setting  the  system  into  intense  oscillations.  The  amplitude 
is  controlled  by  damping  and  nonlinearity  effects,  which  often 
cause  the  shear  layer  to  roll  up  into  discrete  vortices.  Typi- 
cally, oscillations  begin  to  appear  at  low  Mach  numbers  and 
become  most  Intense  slightly  above  M = 1.0;  at  very  high  Mach 
numbers,  say  M > 3-0,  the  oscillations  are  often  found  to 
decrease . 


6.2  Prediction  of  Tonal  Frequencies 

The  oscillation  frequencies  of  shallow  cavities  (say 
L/D  >2.0)  can  be  predicted  with  reasonable  accuracy  by  the 
modified  Rossiter  Equation  [12^1], 


f = ^ 
n L 


(n  - 0.25) 


M 


, n = 1,2 ,3.  • . 


(1  + 1^  MM 


V ' 


+ 1.75 


(92) 


where  f^  denotes  the  frequency  of  the  nth  mode,  U represents 
the  free  stream  speed  M the  Mach  number,  L the  cavity  length, 
and  Y = 1.^  is  the  ratio  of  specific  heats  of  air.  This  for- 
mula is  derived  from  a seml-einpirlcal  form  originally  suggested 
by  Rossiter  [12‘j],  modified  to  account  for  the  nearly  adiabatic 
temperature  recovery  in  the  cavity. 


In  Figure  11  refjults  calculated  from  Equation  (92)  are 
compared  to  the  best  available  data.  One  may  note  that  the 
predictions  best  match  the  data  lor  M > 0,5-  For  M < 0.5>  the 
data  and  Equation  (92)  agree  loss  well,  and  thore  is  some 
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evidence  oT  a systematic  L/D  dependej-kce ; however,  one  should 
note  that  the  low  Mach  numbe:*  date  is  all  based  on  tests  of 
rather  small-scale  models. 


For  deep  cavities  that  respond  primarily  in  depth  modes, 
the  oscillation  frequencies  (at  least  for  low  Mach  numbers) 
may  be  expected  to  correspond  to  the  resonance  frequencies  of 
a pipe  that  is  open  on  one  end  and  closed  on  the  other,  namely 


f « 


2ri-l  c 
IT  D 


n = 1,?,3. . • 


(93) 


where  c denotes  the  speed  of  sound  iii  air.  At  higher  Mach 
numbers,  the  shear  layer  over  the  cavity  mouth  may  contribute 
an  effective  stiffening  effect  at  the  open  end,  leading  to 
deviations  from  Equation  (93)* 


Frequency  prediction  is  most  uncertain  for  cavities  that 
are  neither  clearly  shallow  nor*  deep,  say  for  1.0  < L/D  < 2.0, 
foi'  which  the  modal  response  is  difficult  to  anticipate  and 
for  which  occurance  of  ”sl  Ilow”  or  "deep"  cavity  behavior  may 
depend  on  Mach  number  (with  increasing  Mach  number  tending  to 
make  cavities  effectively  more  shallow)  [126]. 


6.3  Prediction  of  Mode  Shapes  and  Level  In  Cavities 

Two  approaches  are  reviewed  here:  (1)  that  of  Shaw  and 

Smith  [12?],  whose  curve  is  fitted  to  /lir  Force  flight  test  data 
and  (2)  that  of  Holler  and  Bliss  f 126‘J,  which  is  based  on  wind 
tunnel  tests.  Both  of  these  involve  the  modal  frequencies  as 
detennined  from  Equation  (92). 


6.3.1  The  Method  of  Shaw  and  Smith 


Accoi’ding  to  this  method,  the  peak  sound  pressure  in 
the  one-third  octave  bands  containing  the  first  three  cavity 
oscillation  frequencies  are  obtained  from  the  following  rela- 
tions : 


?0  log  &^/q)  = 9.0  -3.3  (L/D)  + 20  log  (-MH2M-0.7) 

(9^la) 

20  log  (?"/q)  = 20  log  (^/q)  - 2(L/D)*  + 26(L/D)  - 86 

1 2 

(94b) 

120  log  (l"yq)  - l],for  L/D  < 4.5 
20  log  & /q)»  for  L/D  > 4.5 

(94c) 

Here  represents  the  maximum  value  of  the  fluctuating  pres- 
sures associated  with  the  nth  mode  (which  maximum  pressure 
occurs  near  the  rear  bulkhead)  and  q denotes  the  free-stream 
dynamic  pressure. 

The  variation  of  the  one-third  octave  band  pressure  ampli- 
tude with  distance  from  the  cavity  leading  edge  (i.e.,  the  pres- 
sure mode  shape)  Is  given  by 

20  log  CP„(x)/q]  - 20  log  C^^/q]  - 10  log  [l.O  + (0.33  L/D 

- 0.60)(l-x/L)  - I cos(o^x/L) I j • , n ■ 1,2,3 

(95) 


where  0^ 

* 3.5 

a 

- 6.3 

2 

a 

» 10.0 

6.3.2  The  Results  of  Heller  and  Bliss 

In  a series  of  wind  tunnel  tests,  Heller  and  Bliss  [126] 
measured  pressures  at  the  front  and  rear  of  a variable  depth 
cavity  and  carried  out  corresponding  mode  shape  surveys.  Some 
of  their  results  are  reproduced  here,  both  as  the  basis  for  an 
alternative  prediction  scheme  and  to  provide  an  indication  of 
the  variation  between  predictions. 

Figures  II5  and  II6  show  one-third  octave  band  pressure 
levels  measured  in  the  cavity  front  and  rear  bulkhead  regions, 
for  three  modes  and  three  length-to-depth  ratios,  as  a function 
of  Mach  number.  Note  that  different  modes  are  dominant  in 
different  parts  of  the  Mach  number  range.  Also  shown  in  these 
figures  is  some  of  the  Air  Force  flight  test  data  on  which  the 
previously  discussed  method  is  based,  as  well  as  data  points 
from  several  other^tests.  The  sometimes  considerable  spread 
in  the  data  suggests  that  factors  other  than  those  current? y 
taken  into  account  in  prediction  schemes  and  scaling  methods 
(e.g.,  Reynolds  number,  cavity  width,  upstream  boundary  layer 
thickness)  may  also  be  important. 

P’igures  II7  and  118  show  measured  mode  shapes  at  Mach  num- 
bers of  0.8  and  1.0,  respectively,  for  length-to-depth  ratios 
of  2.3,  ^.0  and  5«1.  In  Figure  117 » Air  Force  flight  test  data 
points  are  again  shown  for  comparison.  By  taking  the  most  rear 
ward  value  shown  to  be  3 dB  above  the  tralling-edge  levels 
given  in  Figure  llG,  one  obtains  modal  pressure  distributions 
that  are  within  the  experimental  error.  In  general,  however, 
it  may  be  more  convenient  to  obtain  mode  shapes  using  the  meth- 
od of  Shaw  and  Smith. 


• 

LE  PC  8'x  6' 
TUNNEL 

o 

MIT  NAV  SUP 
SON  TUNNEL 

o 

CM 

• 

o 

FLIGHT  TEST 

-40 

■ 

SMALL  SCALE 

MODEL  TEST 

-80 

I I I I 


t 2 
Mco 


FIG.  11 S.  MACH  NUMBER  DEPENDENCES  OF  RESONANT  MODE  LEVELS  IN 
UAUINfi-Em  AREAS  OF  CAVITIES  [126^. 


^^1 


msi 


MODE  3 MODE'2  MODE  1 


SHAPES  FO 
.8  (conti 


6.4  Prediction  of  Broadband  Spectra  In  Cavities 

According  to  Shaw  and  Smith  [127],  one  may  determine  the 
one-third  octave  band  spectrum  from  Figure  II9,  with  the  pres- 
sure Pjj  corresponding  to  the  peak  of  this  spectrum  given  by 

20  log  CP.  (x/L)/q]  » 20  log  ("P  /q)  + 3^3(L/D)  - 28 
o i 

+.3(1-L/D)  (1-x/L)  (96) 

One  may  note  that  Equation  (96)  provides  for  dependence  on 
length-<’o-depth  ratio  and  on  distance  from  the  cavity  leading 

edge,  and  that  P /q  can  be  determined  from  Equation  (9i|a). 

2 

The  narrow  region  into  which  the  broadband  spectra  ob- 
tained by  Heller  and  Bliss  collapsed  for  fairly  wide  ranges  of 
Mach  numbers  and  length-to-depth  ratios  is  illustrated  in 
Figure  120.  Also  shown  in  this  figure  is  the  result  obtained 
by  substituting  Equation  (9^a)  into  (96),  which  then  becomes 
independent  of  L/D  for  x/L  » 1.0  (i.e.,  for  the  trailing  edge 
region).  The  agreement  may  be  seen  to  be  quite  reasonable. 

6.5  Consecutive  CAvities 

Closely  spaced  cavities  arranged  one  behind  the  other 
sometimes  occur  in  aircraft  designs.  In  such  arrangements , one 
may  expect  to  be  able  to  treat  the  upstream  cavity  as  being  un- 
affected by  the  presence  of  the  downstream  cavity,  but  one  may 
expect  the  downstream  cavity  to  be  affected  by  the  upstream  one. 
The .downstream  cavity  is  likely  to  be  subjected  bo  a much  thicker 
incident  boundary  layer  due  to  the  upstream  cavity,  and  this 
incident  flow  will  contain  strong  periodic  disturbances  due  to 
the  oscillation  of  tne  upstream  cavity.  However,  since  thick- 
ening the  incident  boundary  layer  tends  to  reduce  pressure  os- 
cillations, while  the  presence  of  upstream  disturbances  at  the 
proper  frequencies  tend  to  increase  these  oscillations,  these 
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FIG.  120.  ENVELOPE  OF  ALL  NONDIMENSIONAL  BROADBAND  SPECTRA  IN  1/3-OCTAVE 

BANDS  AT  CAVITY  TRAILING  EDGE,  AS  FOUND  BY  HELLER  AND  BLISS  [126] 
COMPARED  WITH  THE  RESULTS  OF  SHAW  AND  SMITH  [127]. 


two  contributions  may  in  effect  cancel  each  other.  If  the 
cavities  are  very  closely  spaced,  the  downstream  cavity  la 
likely  to  "lock-in"  on  the  upstream  one,  so  that  the  two  will 
oscillate  almost  exactly  in  phase. 

Although  consecutive  cavities  have  not  been  studied  ex- 
tensively, some  corresponding  leading  and  traJlirg  edge  spec- 
tra were  measured  [1P6]  in  consecutive  cavities  with  L/D-3- 
The  coi’responding  spectra  in  the  front  and  rear  cavities  were 
found  to  be  very  nearly  identical,  particulai’ly  as  regards 
the  broadband  portions.  Figure  1?1  shows  the  observed  differ- 
ences between  the  front  and  rear  cavities  levels  in  the  tran- 
sonic and  supersonic  Mach  number  range.  The  data  indicates 
that  (1)  the  differences  between  the  pressures  in  consecutive 
cavities  may  not  be  and  (2)  there  is  a general  trend 

towards  higher  pressure  levels  in  the  rear  cavity  with  increas- 
ing Mach  number. 

6.6  Effect  of  Internal  Stores 

The  presently  available  information  is  not  sufficient  for 
the  systematic  prediction  of  the  effect  of  inserting  a given 
store  into  a given  cavity,  particularly  since  an  almost  infinite 
variety  of  store  and  cavity  configurations  is  possible.  This 
section  therefore  only  reviews  the  available  results,  in  order 
to  indicate  the  nature  of  the  effect  and  to  provide  very  general 
guidelines . 

The  insertion  of  an  Internal  stoi’e  into  a cayity  has  often 
been  observed  to  produce  a slight  reduction  in  the  oscillatory 
pressure  levels.  [l^i'l,  128,  1291.  The  reduction  is  more  dra- 
matic if  the  store(s)  interfere  directly  with  the  cavity  shear 
layer  over  a substantial  fraction  of  the  cavity  length  and  span. 
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Apparently,  this  Interference  tends  to  reduce  the  shear  layer 
motion  that  is  essential  to  the  oscillation  process.  Flpure 
122  summarizes  the  results  of  Heller,  Holmes,  and  Covert  Cl 30] 
for  one  and  two  cylindrical  oclve  nose-cone  stores  in  a simple 
rectangular  cavity.  Tests  were  performed  with  the  storo(s) 
attached  to  the  movable  floor  of  the  model.  The  distance  be- 
tween the  top  of  a store  and  the  cavity  mouth  varied  as  a 
function  of  L/D;  at  L/D*7«0  the  outer  surface  of  a store  was 
just  in  line  with  the  cavity  mouth.  Notice  that  in  some  cases 
the  presence  of  a store  actually  caused  a slight  increase  In 
fluctuating  pressure  level. 

6.7  Noise  Radiation  To  Outside  of  Cavities 

Relatively  little  quantitative  information  is  available 
regsirdlng  tne  noi-ie  radiation  from  cavities.  Pressure  levels 
on  (!',ternal  surfaces  immediately  upstream  and  downstream  of 
cavities  have  been  measured  by  several  investigators  T 
129 » 131]-  Their  results  indicate  that  the  levels  mcasui-od  on 
adjacent  ups v ream  and  downstream  surfaces,  at  distances  not  ex- 
ceedJr>c  the  cavity  depth  dimension,  ai*e  very  no.;i-ly  equal  to 
tho.se  jT'-asured  in  the  cavity  at  the  front  and  rear  bulkh'  Is. 
(Note  that  disturbances  can  travel  upcti’can:  threurh  the  sub- 
sonic portion  of  the  boundary  layer,  even  in  supciTonl c I'low, 
and  intense  upstream  dlsturbanceo  were  fo\».nd  to  persist  at 
supersonic  speeds.)  Pressure  measurements  on  surfaces  to  the 
sides  of  cavities  appear  not  to  have  been  made,  but  levels 
commensurate  with  adjacent  points  in  the  c:»vity  cari  be  expected. 
Thus,  proD-nurcs  on  surfaces  around  an  o;;c i 1 Tat  1 ng  cavity  may  be 
expected  to  be  dominated  by  the  effects  of  that  cavity,  aad  the 
pressure  levels  may  be  expected  to  bo  very  nearly  equal  to  those 
within  the  cavity,  for  a distance  of  at  least  ouo  cavity  depth 
dimension  around  the  cavity  mouth. 
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Schlleren  pnotcgraphs  of  the  external  radiation  pattern 
from  cavities  have  been  made  by  several  researchers,  e.g., 
[l?5]>  [132  J,  At  supersonic  epeed-'i,  the  radiation  Is  ex- 
tremely directional,  wit  the  most  Intense  region  falling  in 
a "beam"  between  the  Macn  cenes  at  the  leading  and  trailing 
edges.  The  most  Intense  portion  of  the  beam  is  likely  to  be 
in  the  direction  normal  to  the  cavity  mouth,  where  it  is  not 
likely  to  affect  aircraft  surfaces.  However,  any  portions 
of  the  aircraft  surface  which  intercept  the  region  between  the 
Mach  cones  may  be  expected  to  be  subjected  to  portions  of  the 
intense  radiated  field,  quite  possibly  for  substantial  dis- 
tances away  from  the  cavity.  Unfortunately,  no  quantitative 
data  is  presently  available  to  allow  estimation  of  this  poten- 
tially important  nearfleld  noise  source.  At  high  subsonic 
speeds  the  radiation  still  appears  directional,  but  is  no 
longer  confined  to  a narrow  beam;  the  most  Intense  radiation 
appears  tc  be  in  the  forward  arc. 

At  low  subsonic  speeds,  a cavity  radiates  sound  nearly 
like  a simple  monopole  [1333.  The  source  of  the  monopole  is 
presuiiiably  the  vnasn  addition  and  removal  process  at  the  cavity 
trailing  edge.  The  mean  square  pressure  p'  radiated  by  this 
monopole  may  be  estimated  from 

«.] ' 

where  q denotes  the  dynamic  pressure,  S » S',  rouhr*! 

number  giver,  by  Figure  11'*  or  Equation  (*^?)»  W repr«miVs  the 
cavity  width  r Is  the  distance  from  the  cavity  trailing 
edge  to  the  obaervatlors  point.  Equation  (9T)  generally  ovi:r> 
estidiutes  tho  presarre  levels;  it  Is  no  2ppi:cubl<t  for  :i>0.3 
or  for  modes  alcove  the  third. 


0,03*^  qS 
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6.8  Suppression  Devices 


Several  devices  which  can  be  used  to  reduce  cavity  oscil- 
lation levels  substantially  were  investigated  by  Heller  and 
Bliss  [126].  As  a general  rule,  all  suppression  devices  are 
more  effective  at  higher  values  of  L/D,  quite  possibly  be- 
cause the  levels  to  be  suppressed  are  generally  less  intense 
in  this  case. 

Suppression  of  pressure  oscillation  in  cavities  was 
achieved  by  placing  small  spoilers  or  vortex  generators  Just 
upstream  of  the  cavity  leading  edge,  which  spoilers  extend 
through  the  upstream  boundary  layer  and  are  set  at  high 
angles  of  attack  so  as  to  disturb  the  incident  flow  signifi- 
cantly and  to  thicken  the  shear  layer  over  the  cavity  mouth. 

In  a Mach  number  range  of  0.8  to  1.5  and  for  length-to-depth 
ratios  between  2.3  and  5*1>  the  dominant  second-mode  tone 
was  always  reduced  by  at  least  10  dB  and  was  often  suppressed 
entirely.  The  first  mode  tone  was  always  suppressed  by  at 
least  5 dB  and  often  entirely.  The  effect  of  the  spoilers  on 
the  higher  mode  tones  was  less  systematic j in  some  cases  they 
produced  substantial  reductions,  in  others,  none. 

Another  effective  suppresaion  means  consisted  of  cutting 
back  the  upper  portion  of  the  cavity  trailing  edge  at  a 
degree  slant.  This  ir«odificatlon  has  the  effect  of  stabiliiiing 
the  shear  layer  impingement  process  at  the  cavity  trailing 
edge.  In  order  to  achieve  this  stabilization,  the  height  of 
the  slanted  region  must  exceed  the  shear  layer  thickness  at 
the  rear  of  the  cavity,  so  that  typically  the  slant  height 
nust  exceed  G.2L.  Vhen  tested  in  the  same  Mach  number  range 
88  the  aforementioned  spoilers,  tl»e  slanted  trailing  edge  re- 
sulted In  reductions  of  the  dominant  second  mode  level  by  at 


least  15  dB,  and  In  reductions  of  5 to  20  dB  in  the  first  mode 
levels.  Some  suppression  of  the  third  mode  levels  was  usually 
obtained  as  well.  A reduction  of  the  broadband  spectrum  level 
of  at  least  5 dB  was  also  achieved,  except  at  high  frequencies 
(typically  beyond  the  third  mode  level). 

Leading  edge  spoilers.  In  combination  with  a trailing 
\ edge  slant,  provided  more  effective  suppression  than  each  of 

I the  two  modifications  used  separately,  but  the  effects  of  these 

I I 

> two  devices  are  not  directly  additive;  often  the  first  three 

I 

' modes  were  suppressed  entirely.  This  combination  was  also 

i 

found  effective  in  suppressing  oscillations  in  consecutive 
cavities . 


Another  suppression  device,  a slotted  trailing  edge  cowl, 
was  found  to  be  quite  effective,  but  its  effectiveness  appears 
to  be  quite  sensitive  to  its  shape  and  location. 


6.9  Illustrative  Calculation  and  Comparison 

6.9.1  Problem  Statement 

Compute  the  trailing  edge  spectrum  and  the  shape  of  the 
second  mode  of  a cavity  10  ft.  long  with  L/D  * 5.0  for  M ■ 1.0 
at  sea  level.  Estimate  the  effect  of  Inserting  two  cylindrical 
stores.  Estimate  the  minimum  reduction  In  levels  obtained  by 
using  simple  oscillation  suppression  devices. 

6.9.2  Results 

At  sea  level  the  density  Is  p ■ 
the  sound  speed  Is  c = 1120  ft/sec. 
ft/sec  and  q = 1^19  Ib/ft^. 

The  resonance  frequencies  found  from  Equation  (92)  are; 

fj  = 31.5  Hz  fg  - 73.6  Hz  f^  » 115.7  Hz 

The  peak  third-octave  band  sound  pressure  levels,  as  determined 
from  Equation  (9^)  are; 

20  log  —17.96 

20  log  (^i/q)  —23.96 

20  log  (^3/q)  —17.96 

The  leading  and  trailing  edge  levelr,  20  log  CPj^(x)/q],  found 
by  use  of  Equation  (95)  with  x/h  ■ 0 and  1,  respectively,  are 
listed  below,  together  with  corresponding  wind  tunnel  results 
from  Figures  II6  and  117;  with  the  latter  shown  In  parentheses. 


n 

Leading 

Mbs. 

1 

- 2B.5 

(-27) 

-18.0 

(-18) 

2 

- 3^.5 

(-33) 

-2H.6 

(-21) 

3 

28.5 

(-35) 

-20.0 

(-23) 

2.38  X 10"^  Ib/ft^  and 
Thus,  at  M = 1.0,  U » 1120 


A.', 
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Broadband  levels  predicted  from  Figure  120  are  shown  In 
Figure  123*  Agreement  between  the  two  methods  of  prediction 
is  seen  to  be^ relatively  good. 

Figure  12*<  shows  the  shape  of  the  second  mode  computed 
from  Equation  (95) » together  with  the  corresponding  experi- 
mental mode  shape  from  Figure  118.  (Since  Figure  Il8  does  not 
show  absolute  levels,  its  level  was  determined  by  extrapolating 
the  shape  to  the  trailing  edge  point,  where  the  level  is  knovm 
from  Figure  II7.)  It  is  apparent  that  the  two  mode  shapes,  al- 
though qualitatively  similar,  show  appreciable  quantitative 
differences  at  some  locations  In  the  cavity. 

The  range  of  effect  of  inserting  two  cylindrical  stores  in 
the  cavity  can  be  estimated  from  Figure  122.  One  would  be  con- 
servative in  assuming  that  there  is  no  effect  at  this  Mach  num- 
ber. However,  depending  on  the  configuration.  Figures  5. 10  sug- 
gests that  a reduction  in  peak  level  of  as  much  as  10  dB  might 
actually  occur. 

Figure  123  shows  that  the  first  three  modal  peaks  extend 
a little  moi’e  than  10  dB  above  the  broadband  spectrum.  The  addi- 
tion of  upstream  spoilers  will  reduce  these  peaks  by  at  least 
5 dB,  and  very  likely  much  more.  A trailing  edge  slant,  or  a 
slant  and  spoiler  combination,  can  be  expected  to  reduce  all 
three  peaks  to  the  broadband ' level. 
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FIG.  124.  ILLUSTRATIVE  COMPARISON  OF  CALCULATED  AND  MEASURED 
SHAPE  OF  SECOND  NODE  OF  A CAVITY. 


6.10  List  of  Symbols  for  Sec.  VI* 


Symbol 

Def  1 n1  tion 

Un1ts+ 

D 

cavity  depth 

ft 

(m) 

L 

cavity  length 

ft 

(m) 

M 

Mach  number 

— 

peak  sound  pressure  in  ntA 
mode 

lb /ft® 

(N/in*) 

S 

Strouhal  number 

— 

U 

free-stream  velocity 

ft/sec 

(m/s) 

w 

cavity  width 

ft 

(m) 

c 

speed  of  sound  in  air 

ft/sec 

(m/s) 

^n 

frequency  of  nth  mode 

Hz 

n 

mode  number 

— 

P* 

mean-square  pressure 

lb /ft* 

(N/m*) 

q 

dynamic  pressure 

Ib/ft* 

(N/m*) 

Y 

ratio  of  specific  heats  of  air 

— 

^Specially  defined  symbols,  such  as  those  for  empirically  deter- 
mined constants,  which  arc  used  only  once  in  t'le  text,  are  not 
included  here.  Such  syihbols  are  defined  in  tre  text  where  they 
occur. 

fThe  units  siver  here  are  typical  ones.  SI  units  are  given  in 
parentheses  vhere  appropriate.  Note  that  sem?  empirical  predic- 
tion methods  require  the  use  of  specific  units. 
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SECTION  VII 


BLAST  PRESSURES  GENERATED  BY  ARMAMENT 
7.1  Overview  of  Literature 

The  key  conclusions  one  may  draw  from  the  available  liter- 
ature (Refs.  [13^1]  “ [1^^])  are: 

(a)  In  order  to  account  for  the  different  basic  blast- 
generation  mechanisms  and  the  correspondingly  different  blast 
fields,  one  may  consider  weapons  in  three  g;roups : 

Closed-breech  weapons 

Open-breech  weapons 

Rockets . 

(b)  Closed-breech  weapons  have  been  studied  widely;  corres- 
ponding scaling  laws  and  blast-field  prediction  techniques  are 
well  developed. 

(c)  Open-breech  weapons  also  have  been  investigated  tc  a 
considerable  extent,  and  much  data  is  available.*  The  back-blast 
through  the  nozzle  at  the  breech  tends  to  be  more  severe  than 
the  muzzle  blast. 

(d)  Although  steady-state  noise  fields  from  rockets  ex- 
hausts can  be  predicted  reliably,  attempts  to  scale  and  pj edict 
blast  fields  from  rockets  generally  ha.>e  been  unsuccessful. 

(e)  For  weapons  fired  from  high-speed  aircraft,  the  for- 
ward motion  of  the  aircraft  results  in  Increasing  the  effective 
blast  energy  forward  of  the  blast  source  and  in  decreasing  that 
aft  of  the  source. 
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7.2  Scaling  of  Blast  Fields  From  Stationary  Explosives 


Scaling  laws  for  blast  fields  about  weapons  have  their 
origin  in  the  scaling  laws  used  for  explosives.  Therefore  It 
Is  helpful  to  review  the  scaling  laws  for  explosives  that  are 
stationary;  the  theory  must  be  modified  considerably  to  ac- 
count for  motion  of  the  explosive.  (However,  moving  explo- 
sions best  describe  the  blast  fields  generated  by  weapons.) 

In  1915>  Hopklnson  established  a pair  of  scaling  laws  relating 
peak  overpressure  and  positive  impulse  to  charge  weight  (or 
charge  energy)  and  standoff  distance.  These  laws  may  be  ex- 
pressed In  the  following  forms: 


AP 

e 

Hi 

(l/w^/3) 

a 

Ha 

(L/W^/-) 

or 

AP 

B 

Hi 

(l/e^''3) 

8B 

Ha 

(L/E^'^^) 

or 

AP 

9 

Hi 

(L/d) 

I/d 

m 

Ha 

(L/d) 

where  IIi 

and  II2  represent  functions 

(98) 


(99) 


(100) 


AP  “ free- field  peak  overpressure  (psl) 

1 >•  positive  Impulse  (area  under  positive  overpressure 
tine  history)  (psl-msec) 

L ■ standoff  distance  (distance  from  explosion  center 
to  observer)  (ft)  « 

W charge  weight  (proportional  to  E and  to  d'^)  (lb) 

B • energy  release  of  explosion  ■ (1,55*10*  ft*lb/lb)* 
W for  TNT  (ft-lb) 

d ■ equivalent  spherical  diameter  of  charge  (ft) 


miT  Is  generally  taken  as  the  standard  against  which  the  "yield" 
of  other  explosive  materials  is  evaluated.  Figure  125  shows 
curves  of  AP  and  as  functions  of  for  TNT  explosions 
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OVERPRESSURE,  AP  (psi) 

SCALED  IMPULSE.  ({psi- msec)/lb 


IV  I IV  IV  iv 

SCALED  STANDOFF  DISTANCE.  L/Vl/'*  (fl/lb’'^) 

FIS.  125.  EMPIRICAL  RELATIONSHIPS  BETWEEN  OVERPRESSURE.  IMPULSE, 
CHARGE  HEIGHT  AND  STANDOFF  DISTANCE  FOR  GROUND  SURFACE 
FIRING  OF  TNT  [134j. 
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occurring  at  the  surface  of  the  ground.  These  curves  were 
developed  as  the  result  of  statistically  averaging  a large 
quantity  of  blast  measurements. 

The  above  scaling  laws,  when  generalized  to  stationary  ex- 
plosions occurring  at  any  altitude  above  the  ground,  where  re- 
duced ambient  air  pressure  effects  must  be  Included,  may  be 
written  as 

AP/Po  = n,  [L/(E/P^)^/3] 

I/(E/Po)^/3  » H2  [L/(E/Pq)^/33 

where  P Is  the  ambient  atmospheric  pressure  (psi).  All  of  the 
above  scaling  laws  correspond  to  free-field  blasts,  for  which 
the  direct  blast  wave  dominates  and  there  are  no  reflections 
from  ground  or  other  surfaces. 

7.3  Scaling  of  Free-Fleld  Blasts  From  Stationary  Weapons 

Scaling  laws  for  free  blast  fields  around  stationary 
weapons  are  developed  in  a manner  similar  to  those  of  Section 
6.2  above,  except  that  the  parameters  are  modified  to  account 
for  the  more  complicated  blast  fields  generated  by  weapons. 

The  scaling  laws  for  weapons  may  be  expressed  in  the  form 
[135] 

AP*C^*b/E  « niCX/C,  R/C] 

IC^a^/E  = naCX/C,  R/C]  (102) 

where  Iliandlla  again  represent  functions  and 

C ■ caliber  ■ inside  diameter  of  barrel 

b « barrel  length 

a^“  speed  of  sound  in  air 

S ■ energy  release  from  gun  blast 

X **  axial  distance  from  muzzle  to  observer  ■ L cos  0 
R ••  lateral  distance  from  muzzle  to  observer  •>  h sin  6 
L " standoff  distance  from  muzzle  to  observer 
9 * angle  between  barrel  axis  and  muzzlo»observer  line 
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7.4  Free-Field  Blasts  From  Stationary  Closed-Breech  Weapons 

Figure  126  shows  a summary  of  side-on  overpressure  data 
points  measured  in  the  muzzle  plane  (x=0)  at  various  lateral 
distances  for  fourteen  different  Army  weapons.  All  data  per- 
tain to  ground  level,  and  thus  do  not  contain  the  effects  of 
altitude,  and  all  pertain  to  short  standoff  distances  of  less 
than  400  calibers.  Scatter  about  the  regression  line  is  seen 
to  be  relatively  small,  considering  the  wide  range  of  caliber 
sizes,  energy  levels,  barrel  lengths,  and  standoff  distances 
included.  Thus,  the  regression  equation  given  in  this  figure 
provides  a fairly  accurate  method  for  predicting  side-on  over- 
pressure levels  for  a wide  range  of  closed-bi’eech  weapons. 

Figures  127  summarizes  Impulse  data  in  the  same  manner 
that  Figure  126  summarizes  overpressure  data.,  again  showing 
tight  clustering  about  the  regression  line.  Although  both  of 
those  figures  only  indicate  variations  in  the  muzzle  plane, 
similar  scaling  appears  to  hold  also  for  variations  in  the 
direction  of  the  barrel  axis,  so  that  one  can  develop  curves 
of  constant  reduced  (seeled)  overpressure  and  of  constant  re- 
duced impulse,  as  indicated  in  Figures  128  and  129. 

7.5  Scaling  Laws  for  Blast  Fields  From  Moving  Explosions 

Thornhill's  theoretical  solution  [136]  for  the  blast  field 
generated  by  a moving  explosion  may  be  expressed  [137]  as 


(10^1) 


SCALED  OVERPRESSURE  AP  * C*  • b /E 


4 10  20  40  100  200  400 


REDUCED  f.ATcRAL  DISTANCE  IN  MUZZLE  PLANE,  R/C 


FIG.  126.  SCALED  SIDE-ON  OVERPRESSURE  FOR  CLOSED-BREECH 
WEAPONS  [135]. 
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FIG.  129.  CURVES  OF  CONSTANT  SCALED  IMPULSE  IC^^^  b^^^/E 
FOR  CLOSED-BREECH  WEAPONS  [1353- 


where 


V 

a. 

K 

n 

E 

E, 

e 

f 


velocity  of  moving  charge 

speed  of  sound  In  the  ambient  atmosphere 

constant  associated  with  explosive  material 
0.157  for  TNT 

explosive  energy  per  unit  mass  of  explosive  material 
effective  stationary  energy  of  moving  explosive 
energy  of  explosive 

azimuth  angle  measured  from  the  direction  of  motion 
parameter  evaluated  from  Equation  (103) 


If,  as  usually  is  the  case,  the  quantities  V,  a^,  K,n  are 
known,  one  may  use  Equation  (103)  to  solve  for  the  parameter 
f.  Then  one  can  determine  the  effective  stationary  energy  E 
corresponding  to  the  moving  charge  from  Equation  (10^1). 

The  effect  of  the  velocity  V Is  to  alter  the  directivity 
pattern  of  the  blast  pressure  (increasing  the  pressure  forward 
of  the  moving  explosive  and  decreasing  It  aft  of  the  explosive) 
and  to  increase  the  net  energy  release  from  the  explosive. 


7.6  Theoretical  Prediction  of  Blast  Field  From  Closed-Breech 
Weapon 


Smith  C137]  ascribes  the  blast  field  about  a stationary 
closed-breech  weapon  to  the  gas  efflux  from  the  barrel  subse- 
^ quent  to  emergence  of  the  projectile.  The  gas  efflux  velocity 
V is  given  by 

V ••  -XT  a + U (105) 

*e  y*l  8 Y+1  s 

where 

Y * ratio  of  specific  heats  for  the  propellant  gas 
a.  *■  speed  of  sound  in  the  compressed  propellant  gas 

* speed  of  projectile  as  it  emerges  from  the  barrel 

8 


(106) 


The  flow  energy  density  n Is  given  by 

K K - < 

" " T \ Trrrrr 


where  a is  the  speed  of  sound  In  the  propellant  gas  at  ambient 

" 2 

conditions.  According  to  Smith,  the  extra  energy  term  V /n  in 

Equation  (103)  and  (10^)  is  not  needed,  because  the  gas  efflux 

velocity  is  obtained  at  the  expense  of  the  internal  energy  of 

the  propellant  gas.  Thus,  the  momentum  equation.  Equation  (103). 

for  gas  efflux  reduces  to 

f(l-0.20f^)  = Ve  a^/Kn  (107) 

and  the  energy  equation.  Equation  (10^1),  reduces  to 

-f-  = f«cos6  + (1-f^-sin^e)^^^  = [P(f,e)]^  (108) 

o 


The  results  of  applying  these  equations  to  a 7.62  mm 
closed-breech  rifle  [137]  are  shown  in  Figure  130,  together 
with  corresponding  experimental  data.  Prom  these  results,  the 
values  of  f and  K foi’  an  equivalent  stationary  explosion  were 
found  to  be  f » 0.70  and  K = O.I85. 

Smith’s  gas  efflux  theory  appears  to  provide  a good  method 
for  theoretically  predicting  the  disti'ibutlon  of  the  blast  field 
about  a closed-breech  weapon.  For  this  purpose,  the  above  value 
of  K might  be  x*easonable  for  different  weapons.  Interestingly, 
this  value  of  K is  within  20%  of  the  value  for  a TNT  explosion. 
Smith  also  uses  the  above  equations,  along  with  experimental 
data,  to  chow  that  the  near field  blast  of  a weapon  depends  on 
gas  efflux,  whereas  the  farfiold  blast  is  determined  only  by 
the  total  energy  release. 
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COMPARISON  BETWEEN  THEORY  AND  EXPERIMENT 
FOR  7.62  mm  RIFLE  [137]. 
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FU.  131.  EFFECT  OF  FORWARD  SPEED  ON  THEORETICAL 
6UN  BLAST  PRESSURE  DISTRIBUTION  [|37]. 


Smith  [137]  also  extends  the  above  theory  to  the  predic- 
tion of  the  modified  blast  field  about  a weapon  having  a for- 
ward velocity^,  '.he  general  effect  of  the  forward  velocity  can 
be  seen  in  Figure  13I,  which  shows  the  Increased  blast  energy 
forward  of  the  gun  and  the  reduced  blast  energy  aft  of  the  gun. 

7.7  Prediction  of  Free-Field  Blasts  From  Recollless  Rifles 

A recollless  rifle  has  a nozzle  at  the  breech,  so  as  to 
allow  a fraction  of  the  propellant  gas  to  oxcape  rearward,  in 
order  to  reduce  the  transient  shock  load  on  the  rifle  support. 
Thus,  a recoilless  rifle  may  be  designated  as  an  open-breech 
weapon . 

In  a recollless  rifle,  a blas1^  field  is  generated  rearward 
of  the  breech,  as  well  as  forward  of  the  muzzle.  Because  of 
the  highly  directional  nature  of  the  sound  waves  emanating  f.^'orn 
these  two  blast  fields,  overpressure  and  impulse  values  forward 
of  the  muzzle  are  dominated  by  the  blast  field  generated  at  the 
muzzle,  whereas  the  blast  field  aft  of  the  breech  is  dominated 
by  that  associated  with  the  breech  [I38].  Thus,  for  purposes 
of  evaluating  measured  data,  scaling  of  the  bast  fields,  and 
blast  field  predictions,  these  two  blast  fields  can  be  treated 
independently . 

7 . 7 . I Seal i ng  Laws 

For  closed-breech  weapons,  the  scaling  law  for  overpressure 
LEquation  (102)]  employs  the  ratio  AP/(E/C^b).  The  quantity 

P 

(E/C  b)  has  the  units  of  pressure;  is  described  in  Reference 
[I38],  for  recoilless  weapons  this  quantity  should  be  replacec 
by  the  chamber  pressure  P for  best  scaling  results,  because  a 
given  recolllass  rifle  can  be  operated  with  different  nozzles 
that  result  in  different  chamber  pressures  and  dlffere  ; blast 


field  characteristics.  The  scaling  laws  that  yield  the  best 
results  [138]  are  of  the  form 

AP/P^  - [x/C,  R/C] 

IC/CP^E^]^/3  - Rg  [x/C,  R/C]  (109) 

7.7.2  Blast  Field  Data 

Reference  [138]  contains  an  extensive  summary  of  data, 
collected  from  eight  different  sources,  for  thirteen  different 
rifle  models. 

Figures  132  and  133  summarize  blast  overpressure  data  aft 
of  the  breech,  for  observation  positions  along  lines  at  two 
angles  from  the  weapon  axis;  Reference  [138]  also  contains  simi- 
lar data  for  other  angles.  Figure  131;  shows  all  of  the  mean 
lines  of  these  data  plots  and  indicates  how  overpressure  de- 
creases with  increasing  angle  from  the  axis. 

Figures  135  and  136  summarize  blast  overpressure  data 
at  three  angles  ahead  of  the  muzzle.  The  mean  lines  of  these 
figures  and  of  other  similar  ones  of  Reference  [138]  essentially 
overlap,  indicating  absence  of  a directivity  effect  (as  for  a 
spherical  shock)  ahead  of  the  muzzle. 

In  addition  to  the  overpressure  data  discussed  above. 
Reference  [I38]  also  presents  similar  data  on  the  impulse;  the 
scatter  here  is  somewhat  greater,  but  the  general  conclusions 
are  the  same. 

Analysis  [138]  of  these  data  correlations  leads  to  the  em- 
pirical prediction  equations  given  below. 
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Forward  of  Muzzle; 


10*P 


11.8 


p 

c 


Aft  of  B**eech: 
10 'P 

TT 

10  I. 


3.35 


08 


(110) 


exp  (5.82  - 0.0<<280)‘fe 


P 

c 


exp  (6.30  ” 0.03620) 


(hT- 


003790  - O.QOil) 
(111) 


Note  that  the  pressure  field  forward  of  the  nozzle  is  indepen- 
dent of  the  anelo  0,  whereas  the  field  aft  of  the  breech  exhi 
bits  an  angle-dependence. 


No  data  is  available  on  the  effect  of  forward  motion  on 
the  blast  field. 
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S^m. 

Weapon 

Sym. 

Weapon 

X 

57  mm  M18A  1 

o 

106  mm  T 1 VOEl 

p 

90  mm  T219  (PAT) 

□ 

106  mm  M40A  1 

o 

75  mm  T21 

A 

105  mm  M27 

105  mm  T19 

0 

57  mm  M18 

0 

75  mm  T21 

o 

105  mm  M27 

57  mm  T66E6 

(modified) 

SCALED  AXIAL  DISTANCE  AFT  OF  BREECH 

X/C 

FIG.  132.  BLAST  OVERPRESSURES  OF  RECOILLESS  RIFLES,  ALONG  THE 
AXIS,  AFT  OF  THE  BREECH  [138]. 
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Sym. 


W capon 


90  mm  T219  (PAT) 
75  mm  T21 
105  mm  T19 
75  mm  T21 
57  mm  T66E6 


Sym. 

W capon 

o 

106  mm  T170E1 

A 

105  mm  M27 

O 

57  mrn  Ml 8 

o 

105  mm  M2  7 

(modified) 

o 

^ qO 

^^0 


0.04 


40U  200  100  40  20  10 

SCALED  STANDOFF  DISTANCE  AFT  OF  BREECH 

L/C 

FIG,  133.  BLAST  OVERPRESSURES  OF  RECOIUESS  RIFLES,  AFT  OF 
BREECH.  AT  30*  FROM  AXIS  [138]. 


> 
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10  20  40  ^0  100 

SCALED  STAHDOFF  DISTANCE,  L/C 


FIG.  135 


OVERPRESSURES  OF  RECOILLESS  RIFLE  AHEAD  OF 
NUZZLE  AT  S*  AND  30**  FROK  AXIS  [138]. 


SCALED  OVERPRESSURE,  AP  lOVPe 


Sym.  W capon 


X 57  mrn  M18A  1 
[P  90  mm  T219  |PAT) 
^ 105  mm  T19 

O 106  mm  T170E1 
O 57  mm  M18A  1 
O 106  mm  M40A1 
O 90  mrn  M67 


FIG.  135.  OVERPRESSURES  OF  RECOlUESS  RIFLE,  TO  SIDE  OF  NUZZLE 


7.8  Li'it  of  Symbols  for  Sec.  VII 


Symbol 

Def i ni ti on 

Units'^ 

C 

caliber 

ft 

(m) 

£ 

energy  release 

ft  lb 

(Nm) 

I 

impulse 

lb  sec/ft ^ 

(Ns/m* ) 

L 

standoff  distance 

ft 

(m) 

Po 

ambient  atmospheric  pressure 

Ib/ft® 

( N/m  * ) 

R 

lateral  distance  from  muzzle 
to  observer 

ft 

(m) 

U 

s 

muzzle  velocity  of  pro- 
jectile 

ft/sec 

(m/s ) 

V 

velocity  of  moving  charge 

ft/sec 

(m/s ) 

Ve 

gas  efflux  velocity 

ft/sec 

(m/s) 

W 

charge  weight 

lb 

(N) 

speed  of  sound  in  air 

ft/sec 

(m/s) 

speed  of  sound  in  propellant 
gas 

ft/sec 

(m/s) 

b 

barrel  length 

ft 

(m) 

d 

equivalent  charge  diameter 

ft 

(m) 

X 

axial  distance  from  muzzle 
to  observer 

ft 

(m) 

Y 

ratio  of  specific  heats  of 
propellant  gas 

— 

n 

explosive  energy  per 
unit  mass 

ft  Ib/siug 

(Nm/kg) 

e 

angle  between  barrel  axis 
and  source-to-observer  line 

deg 

AP 

peak  overpressure 

Ib/ft^ 

( N/ni  * ) 

*Speclally  defined  symbols,  such  as  those  for  empirically  deter- 
mined constants,  which  are  used  only  once  in  the  text,  are  not 
Included  here.  Such  symbols  are  defined  in  the  text  wh  re  they 
occur. 

tThe  units  given  here  are  typical  ones,  SI  units  are  given  in 
parenthesos  where  appropriate.  Note  that  .^omc  empirical  predic- 
tion methods  require  the  use  of  specific  units. 
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